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INTRODUCTION 


This  report  includes  abstracts  and  bibliographic  lists  on  major 
contractual  subjects  that  were  completed  in  May,  1972.  The  major 
topics  are:  laser  technology,  effects  of  strong  explosions,  geosciences, 
and  particle  beams.  A  report  on  biocybernetic  studies,  selected  as  the 
optional  topic  for  May,  will  be  issued  under  separate  cover.  The 
abstracted  material  includes  some  selections  prior  to  1972  that  have 
not  otherwise  been  reported. 

To  avoid  duplication  in  reporting,  only  laser  entries  concerning 
high-power  effects  have  been  included,  since  all  current  laser  material 
will  appear  routinely  in  the  quarterly  bibliographies. 

An  index  identifying  source  abbreviations  and  an  author  index  to  the 
abstracts  are  appended. 
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1  •  Laser  Technology 


A.  Abstracts 


Ashmarin,  I.  I.,  Yu.  A.  Bykovskiy,  V.  A. 
Gridin,  V,  F.  Yelesin,  A.  I.  Larkin,  and 
I.  P.  Sipaylo.  Shock  waves  generated  by  the 
action  of  laser  radiation  on  transparent  bodies. 
IN:  Sbornik.  Kvantovaya  elektronika.  Moskva, 
Izd-vo  Sovetskoye  radio,  no.  6,  1971,  126-128/ 


laser  rad?af  lul  exp*rioments  ls  described  on  interaction  of  powerful 

r riter'a  f  u"  yP®  8  glaas’  with  the  obiect  of  determining  the 

criteria  for  shock  wave  generation  at  the  focal  point.  A  Q-switched  ruby 

was  used  generating  50  Mw  pulses,  focused  to  a  0.  2  mm  spot  both  internally 
wa8°u"«d  tg;7  8UrfaCe'  A  plus  high-speed  framing  method  " 

the  tra^sitio  fr  rmi"!  WaV<!  propaSali°n  velocity.  Results  confirmed  that 
he  transition  from  a  longitudinal  sonic  wave  to  a  shock  wave  generally 

occurred  near  the  damage  threshold.  Wave  velocity  was  observed  to  be  a 

asnsee„"ta  FiUg“eimt,e',thy’  “  T'"  °f  dela’r  Ume  in  P“>“  >PP»cation, 

in  Fig.  1.  I„  theory  the  sonic  wave  converts  to  a  shock  wave  when 
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Fig.  1.  Wave  velocity  vs.  delay  (squares) 
and  laser  pulse  intensity  (circles). 


pressure  p  in  the  focal  region  attains  a  value  of  p 
density  and  CQ  is  the  sonic  velocity  in  the  glass. 


QC0  ,  where  pQ  is  glass 
A  sample  calculation  based 
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on  actual  parameters  was  made  on  this  assumption,  showing  that  shock 
wave  velocity  should  be  1.3  times  sonic;  this  agreed  closely  with  measured 
results.  With  the  beam  focused  on  the  surface  a  dual  spherical  wave  was 
typically  generated,  as  shown  in  Fig.  2.  Their  characteristics  were  essentially 
the  same  as  for  the  internally  focused  pulse. 


Fig.  2.  Dual  wave  generation  from 
surface  irradiation. 


Cherkun,  Yu.  P.  ,  and  I.  N.  Konopel'ko. 
Current  status  of  experimental  and  theoretical 
studies  in  the  field  of  combined  heat  exchange. 
I-FZh,  v.  22,  no.  4,  1972,  757-758. 


A  brief  summary  is  given  of  several  papers  on  the  subject 
of  radiation  -  induced  heat  exchange,  presented  at  a  session  of  a  general 
conference  on  mass  and  heat  transfer  held  on  Sept.  15-17,  1971  at  the 
Thermophysics  Institute,  Siberian  Branch  AN  SSSR.  The  general  theme  of 
the  session  was  phenomena  related  to  combination  heat  transfer  mechanisms, 
namely  radiative  -  conductive  and  radiative  -  convective. 


S*  S*  Kutateladze  and  D.  I.  Avaliani  reported  on  experimental 
findings  of  laser  beam  attenuation  in  liquid  media  with  pulsating  turbulence; 
this  work  was  judged  important  from  a  heat  transfer  viewpoint  as  well  as  in 
regard  to  beam  propagation  in  a  turbulent  medium.  A  related  paper  by 
.  .  Kostylev  and  V.  Ya.  Belostotskaya  dealt  with  the  relation  between 

conductive  and  radiative  heat  transfer  in  optically  thin  layers  of  finely-dispersed 
media;  they  show  that  in  this  case  it  is  incorrect  to  use  the  hypothesis  of 
additive  fluxes. 


A  report  by  N.  A.  Rubtsova  and  A.  E.  Berte  covered  the 
t  ermal  state  of  pure  metal  surfaces  exposed  to  a  constant  thermal  flux. 
Results  showed  anomalies  in  surface  temperature  distribution  (the  Jacques 
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effect)  on  metals  including  Al,  Fe,  Sn  and  others;  several  possible  models 
for  this  phenomenon  are  advanced. 

V.  M.  Yeroshenko  discussed  variants  of  ablation  techniques 
in  which  blowoff  of  dispersed  media  is  used  as  a  protection  against  radiative 
flux;  this  refers  to  processes  at  levels  of  100  atm  and  temperatures  in  the 
20,000°C  range.  Several  papers  by  other  authors  dealt  with  heat  exchange 
studies  in  various  types  of  combustion  reactions. 

Special  mention  is  giver,  to  a  paper  by  A.  I.  Leont'yev  and 
A.  M.  Pavlyuchenko  which  is  cited  as  giving  the  most  complete  picture 
to  date  on  temperature  distribution  in  a  turbulent  boundary  layer  and  laminar 
sublayer  in  a  thermal  radiation  field.  As  with  Kostylev  et  al,  these  authors 
also  demonstrate  the  fallacy  of  the  additive-flux  hypothesis  under  their  test 
conditions. 


In  closing,  S.  Kutateladze  stressed  the  need  for  further  experi¬ 
mental  data  on  radiative  hnat  exchange  resulting  from  interaction  of  e-m 
waves  with  materials. 


Epshteyn,  E.  M.  Thermal  instability  of  a 
semiconductor  in  a  laser  beam.  IVUZ  Radiufiz, 
no.  1,  1972,  33-37. 

A  theoretical  study  is  presented  on  a  possible  mechanism  of 
laser-induced  destruction  of  a  semiconductor.  The  case  is  limited  to  photon 
energies  less  than  the  width  of  the  forbidden  zone,  so  that  the  main  exchange 
mechanism  is  absorption  by  free  carriers.  In  the  general  case  a  cumulative 
heating  and  increase  in  free  carriers  will  occur,  resulting  in  a  stationary 
temperature  field.  However,  above  some  laser  threshold  level  the  stationary 
field  breaks  down  and  a  rapid  thermal  reaction  sets  in  which  is  analogous  to 
breakdown  in  a  d-c  field,  or  the  thermal  explosion  from  an  exothermic 
chemical  reaction.  To  examine  the  governing  phenomena  the  author  assumes 
the  simplest  case  of  a  cylindrical  semiconductor  target,  axially  excited  by 
a  uniform-intensity  laser  beam  such  that  \«  p  «  R,  where  A  is  laser 
wavelength,  and  p  and  R  are  beam  and  cylinder  radii.  Expressions  are  derived 
for  the  limit  conditions  of  the  stationary  field  solution,  and  the  corresponding 
threshold  power  of  the  laser  is  determined.  Finally,  it  is  shown  that  the 
elapsed  time  required  to  reach  threshold  under  practical  conditions  will  lie 
in  the  1  millisec  -  1  sec  range  -  i.  e.  demanding  in  effect  a  c-w  laser  regime. 
Therefore  the  described  mechanism  cannot  account  for  breakdown  observed, 
for  example,  in  transparent  dielectrics  under  nanosecond  pulse  exposure. 


-3- 


Frolov,  V.  V.  Temperature  fields  in  multi¬ 
layer  translucent  coatings  under  conditions  of 
pulsed  radiation  heating.  I-FZh,  v,  22,  no.  4, 

1972,  755-756.  (Annotation  of  deposited  paper.) 

A  summary  is  given  of  an  analysis  of  heat  transfer  in  a  multi¬ 
layer  coating  system  under  beamed  irradiation.  The  problem  is  treated 
as  nonstationary  and  one-dimensional,  with  distribution  of  absorption  taken 
into  account.  Reflection  and  concentrated  absorption  at  layer  interfaces 
are  treated;  the  absorption  process  is  considered  to  be  linear. 

As  an  example  the  temperature  distribution  is  calculated  in 
a  two-layer  heat  shield  in  which  the  heat  load  is  assumed  to  be  applied  in 
triangular  pulses  of  a  given  energy,  at  durations  ranging  from  10  to  60 
seconds.  Owing  to  changing  boundary  conditions,  it  is  shown  that  the  system 
must  be  solved  as  a  class  of  piecewise-discontinuous  functions.  On  this 
basis  graphical  solutions  are  obtained  showing  the  time  characteristic  of 
heat  absorption  as  a  function  of  absorption  coefficient  and  pulse  duration. 
Instantaneous  temperature  distributions  are  also  presented  for  the  cases  of 
absolutely  opaque  as  well  as  semitransparent  coating  systems. 

The  analysis  shows  that  failure  to  account  properly  for  internal 
heat  absorption  can  lead  to  significant  errors  in  theoretical  temperature 
fields  and  level  of  total  absorption;  furthermore,  the  error  increases  with 
reduction  in  pulse  width. 


Bessarab,  Ya.  Ya.  ,  Yu.  B.  Tkach,  V.  P. 

Zeydlits,  N.  P.  Gadetskiy,  and  V.  V,  Dyatlova. 

Study  of  collective  processes  in  a  plasma  using 
light  of  stimulated  emission.  ZhETF,  v.  62, 
no.  2,  1972,  569-572. 

The  use  of  coherent  light  emission  to  study  the  collective 
processes  in  a  plc.3ma-beam  discharge  is  described.  It  is  shown  that  an 
effective  method  of  studying  these  processes  is  the  direct  observation  of  the 
form  of  oscillation  in  the  intensity  of  spontaneous-emission  spectral  lines, 
with  a  subsequent  application  of  Fourier  analysis.  Study  of  the  oscillation 
spectrum  was  performed  through  investigation  of  the  fluctuation  intensity  of 
cohorent0emission  in  the  Ar  II  transitions  4p2Dgy2  — ►  4s2P3/2  at  a  wavelength 
of  4880  A,  and  4p^D^/2  *"  at  a  wavelength  of  5145  R.  The  maximum 

frequency  excited  in  the  plasma,  which  depends  upon  modulation  of  the 
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emission  intensity,  is  determined  by  the  relationship  wmaxT^  where  t  is 
the  lifetime  at  the  upper  level  of  the  observed  transition.  In  this  study 
0)  max  ~  100-120  MHz.  In  the  case  of  coherent  emission,  the  spectrum 
of  the  oscillation  being  investigated  has  a  lower  limit  given  by 

_2t9E*c(i~R) 

<0m<B  *  (Nt-Nt)h«>ltd' 

where  €Q  is  the  dielectric  permeability  of  the  active  medium,  E  is  the  field 
intensity  of  the  light  wave  in  the  resonator,  c  is  the  speed  of  light,  R  is  the 
reflection  coefficient  of  the  resonator  mirrors,  d  is  the  distance  between 
mirrors, “H  is  Planck's  constant,  (o ^  is  the  frequency  of  the  observed  transition 
and  Nj  and  N2  the  corresponding  population  densities  of  the  upper  and  lower 
laser  levels,  respectively.  Fig.  1  compares  obtained  spectra  of  coherent 
and  spontaneous  emission. 


S,  rel.  units 


Fig.  1.  Realization  and  type  of  spectrum 
obtained,  a  -  by  modu’  *tion  of  coherent 
emission,  b  -  by  modulation  of  spontaneous 
emission. 
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Sultanov,  M.  A,  Destruction  of  transparent 
dielectrics  under  the  action  of  free-running 
neodymium  and  ruby  lasers .  Mekhanika 
polimerov,  no.  2,  1972,  359-360. 

The  author  notes  that  a  variety  of  theories  exist  on  the  basic 
destruction  mechanism  of  laser  interaction  with  transparent  solids;  among 
these  are  thermal  explosion,  acoustic  phonons  from  stimulated  Brillouin 
scattering,  and  photolysis  of  the  target  material.  Experiments  are  described 
which  were  designed  to  clarify  the  predominant  mechanism  of  laser  damage 
in  several  dielectrics,  including  polymethyl-methacrylate,  polystyrene  and 
several  types  of  glass.  Both  ruby  and  neodymium  crystal  lasers  were  used, 
in  free-running  regimes  at  energies  of  3  to  15  j.  Characteristic  platelets  of 
cracks  at  angles  up  to  90°  from  the  beam  axis  were  observed  for  both  laser 
types;  however  in  the  Nd  target  specimens  the  damage  area  was  more 
extended  along  the  beam  axis,  whereas  the  ruby  damage  was  localized  near 
the  beam  focus.  (Repeated  reference  is  made  to  damage  photos,  which  are 
unfortunately  omitted  from  the  text.  v 

The  test  results  lead  Sultanov  to  the  following  concept  of  the 
damage  process.  On  focusing  a  sufficiently  powerful  beam  in  the  transparent 
dielectric,  a  local  region  is  heated  to  a  liquid  state  in  which  further  energy 
is  absorbed  by  inhomogeneities,  and  crack  propagation  begins.  A  plasma  is 
generated  which  expands  at  high  velocity;  thus  the  effect  is  that  of  a  powerful, 
instantaneous  point  explosion.  The  process  is  a  hydrodynamic  one,  including 
shock  wave  generation  which  causes  further  destruction  both  before  and 
behind  the  beam  focal  point.  Examples  of  remote  damage  regimes  in  PMMA 
specimens  are  cited  in  added  support  of  the  theoretical  model. 
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2.  Effects  of  Strong  Explosions 


A.  Abstracts 


Aslanov,  S.  K.  and  I.  P.  Kopeyka.  Generalization 
of  triple  discontinuity  configurations.  IN:  Sb. 

Matematika  i  mekhanika  tezisy  dokladov  4-y 
Kazakhstan  mezh,  -vuz.  nauch.  konf.  po  mat  i 
mekh.  Ch.  2.  Alma-Ata,  1971,  171-173. 

A  theoretical  generalization  of  triple  discontinuity  configurations 
was  undertaken  to  aid  in  calculating  spin  detonation.  The  generalization 
consists  of  introducing  into  tie  shock  system  an  additional  self-similar 
rarefaction  wave  following  a  Chapman- Jouguet  self-sustaining  detonation. 
Three  possible  variants  of  the  generalization  are  analyzed.  Representing 
different  contact  angles  between  the  reflected  wave  and  the  incident  flow,  the 
variants  are  extended  to  the  case  when  the  Mach  wave  is  a  self-sustaining 
detonation  front  and  the  reflected  front  is  a  supercompressed  detonation.  The 
calculated  pressure  values  near  the  triple  point  agree  well  with  data  from 
experimental  oscilloscope  traces. 


Balakin,  V.  B.  and  V.  V.  Bulanov.  A _ 

numerical  solution  of  the  problem  of  interaction 
of  a  shock  wave  with  a  cylindrical  body  in 
supersonic  flow.  PFZh,  v.  32,  no.  6,  1971, 


1033-1039. 


A  two-dimensional  finite  difference  scheme  of  the  second  order 
accuracy  is  introduced  to  compute  the  pr  ocess  of  collision  of  a  nonstationary 
shock  wave  with  a  stationary  flow  around  the  nose-end  of  a  cylinder  moving 
at  a  supersonic  speed.  The  initial  set  of  equations  of  two-dimensional  flow 
in  divergent  form  is  presented  in  vector  notation 


Of 

01 


V-  /••(/)  :•  v-fif/'M-  ---"W  0 

Ox  Oy  tj 
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and  allowance  is  made  for  through  computation  of  shock  waves.  Eq.  (1)  is 
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then  approximated  in  two  steps  by  a  finite  difference  equation  which 
contains  artificial  viscosity  terms  in  the  x  and  y  directions  measured 
along  and  from  the  flow  axis,  respectively.  The  necessary  and  sufficient 
stability  criteria  of  the  finite  difference  scheme 


were  established  by  a  Fourier  method.  A  =  dF/df  and  B  =  dG/df  in  (2) 
are  Jacobians  of  the  F  and  G  vectors  in  a  random  point  of  the  flow,A  x  and 
Ay  are  the  difference  grid  pitches  along  x  and  y,  A  t  is  the  time  pitch,  and 
E  is  the  matrix  unit.  By  use  of  a  through  computation  method,  a  common 
numerical  solution  is  outlined  for  two  corollary  problems:  the  initial 
supersonic  streamline  flow  around  the  body,  and  the  oncoming  shock  wave 
with  establishment  of  a  new  regime  of  supersonic  streamline  flow.  The 
initial  flow  parameters  are  those  of  the  first  shock  wave.  The  parameters  of 
the  second  shock  wave  are  calculated  by  means  of  Hugoniot  formulas.  Numer 
ical  solutions  of  the  problem  are  illustrated  graphically. 


Limarev,  A.  Ye.  and  A.  D.  Chernyshov. 

Propagation  of  shock  waves  in  a  reinforced 
elasto- plastic  medium.  PMM,  no.  6,  1971, 

1083-1088. 

Propagation  and  structure  of  shock  waves  in  a  translationally 
-reinforced  elasto-plastic  medium  is  described  by  the  closed  set  of  discon¬ 
tinuity  equations 


[3„*  ]  —  2pj  \eu*  ]  =  2  l«pl  (2p1me*i/  -  2p03*^  —  PuP* \e„*  ]  + 
4*  Po  +  •  '• 
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[3iy]  Vy=  —  p+GlPj),  Ittj,  /]  =  0>iV/,  [yj  =  —  C(0t  (4) 
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where  or-,  a  ^  are  stress  components,  e^j  are  strain  components,  X^, 
p»  /^1»  /*2  are  elas1-ic  coefficients,  v  j  is  normal  to  discontinuity  surface, 

G  is  the  propagation  velocity  of  a  nonremovable  discontinuity  surface  in  a 
continuous  medium,  Uj,  Vj  are  the  displacement  and  velocity  of  a  particle, 
and  the  asterisk  denotes  a  tensor  deviator  component.  Equations  1-4  were 
derived  assuming  a  linear  dependence  between  two  discontinuous  functions 
within  the  shock  wave  transition  layer.  This  linear  approximation  made 
possible  a  common  solution  to  the  problems  of  propagation  and  structure  of 
shock  waves  in  a  plastic  noncompressible  material. 

Using  (1)  and  (4),  an  expression  is  derived  for  the  characteristic 
of  irreversible  plastic  deformation  of  a  continuous  medium 

2  (q>]  =  />„  (n0(o,A  ~ 

I^i  =  P+G3  —  —  2nlt  Ft  =  p +G2  —  Xj  —  3/3 Hj  (1  n2/p0) 

where  SP*  is  the  stress  component  of  a  plastic  element  p  of  the  rheological 
model  of  the  material,  the  index  n  and  the  sign  +  indicate  direction  of  the 
normal  v  ^  and  the  value  p  ahead  of  the  shock  wavefront,  respectively.  The 
velocity  p+G2  of  discontinuity  surface  propagation  is  determined  by  the  roots 
of  the  equation  of  the  4th  degree  in  p+G2,  which  was  derived  from  the 
equations  1-5.  The  root  p+G2  - of  the  cited  equation  corresponds  to  a 
transverse  neutral  wave  without  discontinuity  of  plastic  deformation  and 
satisfies  the  condition  £  J  =  0.  This  particular  case  is  illustrated  by  an 
example  of  a  rectangular  elasto-plastic  plate  subjected  to  either  a  tensile  or 
a  shear  stress.  Two  other  roots  of  the  4th  degree  equation  indicated  the 
possibility  of  generating  two  plastic  shock  waves  in  a  three-dimensional  re¬ 
inforced  elasto-plastic  medium.  One  of  the  plastic  wave  j  propagates  in  a 
medium  under  spherically  symmetric  stress,  the  other  in  a  medium  where 
S*?*+  =  0. 


Pogorelov,  V.  I.  The  impact  of  an  incalculable 
supersonic  jet  on  a  plane.  IFZh,  v.  21,  no.  5, 

1971,  941-942. 

Shock  wave  separation  occurs  in  front  of  a  plane  positioned 
normally  to  the  nozzle  axis  when  an  incalculable  supersonic  jet  impinges  on 
the  plane.  The  separated  wave  is  convex  toward  the  plane  and  interacts  with 
the  hanging  shock  of  the  jet  if  the  nozzle-to-plane  spacing  is  such  that  the 
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separated  wave  interferes  with  the  expansion  of  a  free  jet.  A  flow  measure¬ 
ment  equation  is  established  for  gas  flow  ahead  of  and  behind  the  separated 
shock  wave  with  allowance  for  a  relationship  in  the  shock  wave  and  on 
assumption  that  the  density  and  the  velocity  component  parallel  to  the  plane 
are  constant  behind  the  separated  wave.  This  flow  equation  is  used  to 
derive  a  differential  equation  which  defines  the  geometry  of  the  generatrix 
of  the  separated  wave.  The  boundary  condition  of  this  equation  is  the  derivative 
of  the  shock  wave  equation  for  a  triple  configuration  of  shock  waves.  The 
coordinate  of  the  triple  configuration  is  determined  by  a  simple  algebraic 
correlation  using  the  numerical  value  of  the  derivative.  Gasdynamic  parameters 
of  subsonic  flow  behind  the  shock  wave  can  be  also  determined,  if  the  relative 
position  versus  the  plane  and  the  shape  of  the  shock  wave  are  known.  The 
calculated  coordinates  of  the  triple  point  and  geometry  of  the  separated  shock 
wave  were  found  to  be  in  good  agreement  with  experimental  data  obtained 
from  flow  shadow  photographs. 


Kozachenko,  L.  S.  and  B.  D.  Khristoforov. 

Parameters  of  a  shock  wave  from  explosion 
at  the  bottom  of  a  water  vessel.  FGiV,  no.  1, 

1971,  127-135. 

Measurements  are  described  of  the  relative  pressure  profile 
p/pi  (t),  duration  T,  and  specific  impulse  I  of  an  explosive  shock  wave  in 
water  vessels  with  three  different  bottom  types.  The  measurements  were 
prompted  by  the  inadequacy  of  known  formulas,  which  describe  a  wave  field 
in  water  without  allowance  for  significant  nonlinear  effects.  The  shock 
waves  were  generated  by  the  underwater  detonation  of  a  spherical  explosive 
charge.  The  shock  wave  parameters  were  measured  with  a  minimum  +  10% 
accuracy  by  pressure  sensors.  Data  were  recorded  by  PID-9  devices  with 
a  2  Hz  -  300  kHz  range.  The  effect  of  the  free  surface  was  eliminated  under 
the  experimental  conditions.  Three  series  of  experiments  were  carried  out 
in  separate  vessels  with  air- water- saturated  sand  (I),  polystyrene  (II),  or 
concrete  bottoms  (III).  The  three  bottom  types  differed  from  each  other  by 
the  nature  of  their  interaction  with  a  shock  wave  in  water.  This  interaction 
is  defined  by  the  relationship  of  sound  velocity  a  in  water  to  the  velocities  Cj 
of  longitudinal  and  C2  of  transverse  waves  in  the  ground.  The  inequalities 
C2<Ci  <  a,  Cj>a>  C2,  and  Cj  >  C2  >  a  characterize  the  I,  II,  and  III  bottom 
types.  Depth  charges  of  100  kg  cast  TNT,  lg  compacted  TEN,  and  0.2  kg 
compacted  TNT  were  detonated  in  3,  1,  and  0.  5  m  deep  basins  with  I,  II,  and  III 
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bottoms,  respectively,  at  relative  distances  H  =  1-11,  1-12,  and  1-4  charge 
radii  from  the  bottom.  The  pressure  sensors  were  located  at  a  depth  h 
and  at  relative  distances  R  =  30-120  for  vessel  I  and  60-120  for  vessel  III 
from  the  blast  center,  and  on  the  periphery  of  a  hemisphere  of  0.45  m  radius 
for  vessel  II  with  the  center  at  a  point  below  the  charge. 

Typical  oscilloscope  traces  of  shock  waves  in  water  (Fig  1  and 
2)  and  the  experimental  t(h)  and  t(  «  )  plots  (Fig.  3  and  4)  are  shown.  The 


Fig.  1.  Relative  pressure  p/pj  versus 
reduced  duration  x^sec/m)  of  explosive 
shock  waves  in  a  water  basin  with  a  sand 
bottom. 


Fig.  2.  Relative  pressure  p/pj  versus 
reduced  duration  t(msec/m)  of  explosive 
shock  waves  near  the  polystyrene  bottom. 
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Fig.  3.  Duration  t(msec/ra)  of  shock 
wave  versus  distance  h  from  the  pickup 
to  the  sand  bottom. 


Fig.  4.  Reduced  duration  T  (msec/m) 
of  a  shock  wave  versus  angle  «  for  an 
explosion  at  a  polystyrene  bottom. 

nearly  parabolic  shape  of  the  p(i)  curves  (Fig.  1)  suggests  a  nonlinear 
interaction  with  rarefaction  waves  propagating  from  the  bottom.  The  latter 
decrease  the  frontal  pressure,  as  shown  by  the  plots  of  p(h).  In  the  case 
of  the  I  type  bottom,  pmax  of  the  reflected  compression  wave,  which  is 
observed  only  at  large  angles  of  incidence  o  of  the  forward  wave  to  the  bottom, 
attains  30%  of  p  in  the  leading  shock  wave.  A  bottom  wave  behind  and  a  leading 
shock  wave  ahead  of  the  forward  wavefront  were  observed  in  the  case  of  the 
II  type  bottom  (Fig,  2).  The  leading  wave  is  generated  only  at  an  angle  •  as  a  t 
which  is  determined  by  the  intersection  of  two  wavefronts.  The  pmax  may  ° 
exceed  by  40%  the  explosion  pressure  in  an  infinite  liquid  volume  due  to  the 


superposition  of  the  leading  and  forward  wave  pressures  (Fig,  2,  d,e).  The 
pmax  near  *  =  *0*20-30°  increases  sharply  because  of  the  interaction 
between  the  forward  and  leading  waves,  as  shown  on  the  p(«)  plot.  Interaction 
of  the  forward  wave  with  the  bottom- reflected  rarefaction  wave  manifests 
itself  in  a  decrease  of  t  with  a  decrease  in  •  (Fig.  4).  The  shaded  area  in 
Fig.  4  <?  signates  the  near-bottom  zone,  where  total  r  of  the  wave  field  is  much 
longer  than  t  of  the  forward  wave. 

Tabulated  experimental  data  (Table  1)  show  that  the  parameters 


h 

•J*- 

— 

Ic%U 

TP» 

MeeXiM 

l  P'  )s 

90,  H=.\ 

R=9Q, 

,  //=--  3 

0 

0.28 

I  0,45 

0,384 

0 

0,38 

0.40 

1 

0,28 

1  0,64 

0.3S4 

0,123 

0,70 

0.55 

2 

0,78 

1  0,75 

0,376 

0.25 

0.78 

0,64 

3 

0,73  | 

1  0.78 

0,376 

0,37 

0,99 

0,72 

R~  C)0,  W=1 

/?  — 120,  Hrr-.] 

0 

0..19J 

0.31 

3.5 

0 

0,33 

0,30 

4 

0.52 

0,58 

2.1 

0  49 

0,41 

0,51 

8 

0,57 

0,71 

0,93 

0,87 

0,46 

0,62 

Table  1.  Parameters  p/pj  and  t  of  explosive 
shock  wave  in  a  water  basin  with  concrete 
bottom.  Footnote:  The  indexes  e  and  p  indicate 
the  experimental  data  and  the  data  calculated 
from  the  formulas  given  in  the  text,  respectively. 

of  the  forward  wave  for  the  III  type  bottom  are  greatly  influenced  by  the 
bottom- reflected  rarefaction  wave.  The  pmax  of  the  forward  wave  decreases, 
when  the  sensor  is  positioned  closer  to  the  bottom,  and  approaches  the  pmax 
value  calculated  using  an  empirical  formula.  However,  the  experimental 
t  disagrees  with  the  T  calculated  from  an  empirical  formula,  because  of  the 
effect  of  a  strong  bottom  wave. 

The  experimental  data  show  that,  at  great  distances  from  the 
blast  center,  the  effect  on  shock  wave  parameters  in  water  of  all  three  bottom 
types  is  analogous  to  that  of  a  free  surface.  The  empirical  formulas  proposed 
can  be  used  to  calculate  shock  wave  parameters.  The  actual  t  of  the  shock 
wave  measured  near  the  bottom  is  usually  longer  than  the  calculated  x  ,  because 
of  bottom  wave  interference.  Pressure  in  the  lateral  and  leading  waves  is 
generally  not  over  20-30%  oi  that  in  the  forward  wave.  The  maximum  pressure 
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near  the  intersection  of  the  leading  and  forward  waves  is  30-40% 
higher  than  the  explosion  pressure  measured  in  an  infinite  liquid  volume. 
This  paper  is  an  extension  of  a  similar  work  by  the  authors  reported  in 
Effects  of  Strong  Explosions,  no.  2,  p.  111. 


Dmitriyevskiy,  V.  A.  ,  V.  I.  Fedulov,  and 
V.  F.  Nikolayeva.  Study  of  the  properties 
of  SFa  and  UFa  in  a  shock  tube.  ZhETF, 
v.  61,  no.  4,  1971,  1427-1433. 

Tests  are  described  whose  object  was  to  study  thermodynamic 
properties  of  selected  gases  under  shock  heating.  The  test  gases  were 
SF ^  and  UF^,,  since  multiatomic  molecules  of  this  type  could  be  expected 
to  exhibit  some  anon  .lous  behavior  during  relaxation.  An  electrically 
discharged  diaphragm-type  shock  tube  was  used,  with  He  as  the  driver  gas. 
Trace  CC>2  and  other  heavy  gases  were  added  on  occasion  to  reduce  the 
velocity  of  the  wave  front.  One  anomaly  observed  was  a  drop  in  thermal 
capacity  of  the  tested  gases  with  increase  in  temperature.  Results  show 
that  both  SF ^  and  UF^  have  at  least  two  discrete  vibrational  relaxation  modes 
under  the  cited  test  conditions. 


Poluboyarinov,  A.  K.  Propagation  of  a 
shock  front  reflected  by  a  blunt  body. 

MZhiG,  no.  2,  1971,  70-77. 

Shook  wave  front  interaction  with  a  blunt  sphere  or  cylinder  is 
analyzed  on  the  assumption  that  the  front  plane  is  normal  to  the  symmetry 
axis  or  plane  of  the  body  and,  in  the  case  of  the  cylinder,  is  parallel  to  the 
body  generatrix.  An  exact  differential  equation  is  derived  to  describe 
nonstationary  separation  of  the  symmetry  center  of  the  shock  reflected  by  a 
static  body.  The  equation  describes  dimensionless  axial  shock  separation  e 
as  a  function  of  the  independent  variable  M  of  the  reflected  front.  The  equation 
also  contains  two  additional  unknown  quantities:  the  derivative  {  dv/  d  x)j  of 
gas  velocity  immediately  behind  the  shock,  and  the  radius  r  of  the  reflected 
shock  front  curvature  along  the  symmetry  axis.  An  exact  formula  was 
derived  for  initial  (  d  v/  d  x)ls!c  as  a  function  of  r*  and  M*.  Approximate 


formulas  define  r*  and  (dr/dt)*  in  the  ?uitial  separation  phase  of  the  front 
from  the  body.  A  fairly  simple  approximate  method  is  introduced  to 
determine  nonstationary  e  for  a  sphere,  on  the  basis  of  the  derived  formulas 
and  the  precise  equation  of  £.  The  method  consists  of  integrating  the 
equation  d  £  /dM  =  ,  where  <Pis  a  function  of  r  and  (dv/dxjj.  A  numerical 

example  shows  that  the  method  correctly  describes  the  initial  separation 
phase,  and  also  its  final  phase  if  a  correction  to  the  ^function  is  introduced. 


Yeliseyev,  Yu.  B.  and  A.  Ya.  Cherkez.  The 
effect  of  an  increase  in  stagnation  temperature 
in  a  streamline  gas  flow  past  a  deep  cavity. 

MZhiG,  no.  3,  1971,  8-18. 

An  earlier  observed  anomalous  increase  in  temperature  of 
the  separated  gas  flow  within  a  cavity  above  the  stagnation  value  T0  is 
discussed,  on  the  basis  of  a  theoretical  model  of  compression  and  decom¬ 
pression  waves  in  a  cylindrical  tube.  Piston  oscillatory  movements  at  the 
open  end  of  the  tobe  model  simulate  conditions  at  the  streamline  flow- separated 
gas  interface.  Characteristics  of  the  nonstationary  self- oscillating  shock 
wave  process  at  the  interface  are  formulated  and  a  numerical  example  is 
given.  A  preliminary  conclusion  was  drawn  that  the  observed  phenomenon 
is  reproducible  with  models  of  different  dimensions  but  is  influenced  by 
increased  boundary  layer  thickness.  Aerodynamic  experiments  with  different 
models  at  M<  1  and  M  =  1.6  and  2.0  are  described.  In  agreement  with 
theoretical  calculations,  the  stationary  air  temperature  at  the  base  of  a 
model  increased  by  a  minimum  of  100°K  above  T0.  The  experimental  data 
confirmed  the  hypothesis  that  the  shock  wave  thermal  effect  and  heat  release 
follow  energy  dissipation  in  shock  waves.  The  intensity  of  the  observed 
effect  depends  greatly  on  the  shape  of  the  leading  edge  and  the  ?ngle  formed 
by  the  model  exit  section  plane  with  the  axis.  The  experiments  suggest  the 
possibility  of  reproducing  the  effects  in  the  surface  cavity  of  an  object  flying 
in  the  atmosphere.  Methods  for  minimizing  the  undesirable  thermal  effect 
are  given. 


Pustovalov,  V.  K.  Self- similar  gas  motion 
behind  a  shock  wave  front  sustained  by  radiation. 
DAN  BelSSR,  v.  15,  no.  12,  1971,  1079-1081. 


Adiabatic  expansion  of  a  cold  perfect,  gas  in  vacuum,  initiated 
by  a  strong  shock  wave,  is  analyzed  with  allowance  for  radiant  energy  q 
supplied  to  (q  >  0)  or  withdrawn  from  (q<0)  the  expanding  gas.  Shock  wave 
propagation  through  the  gas  is  sustained  by  radiation  which  determines  the 
motion  pattern  of  the  shock  front.  The  radiation  flux  density  is  defined  by 
q  =  qQt*  ,  where  i  is  an  integer.  This  is  a  self- similar  gasdynamic  problem 
of  the  first  kind,  since  there  are  two  independent  dimensional  parameters 
qQ  and  pQ.  The  problem  is  described  by  the  set  of  equations 
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where  P\,  pj,  vj  are  the  unknown  dimensionless  density,  pressure,  and 
velocity  of  the  gas,  f  is  a  self- similar  variable,  Y  is  the  adiabatic  exponent, 
=  1  +  f/3  is  the  self- similarity  index,  k  =  // 3,  and  m  =  2  /  /3.  The 
boundary  conditions  at  the  shock  front  and  on  the  vacuum  side  complete  the 
formulation  of  the  problem. 


Solution  of  (1)  with  the  boundary  conditions  is  given  as  follows: 
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7  — 5y 

8(2  v)s(V  —  1) 


(3) 


where  a  is  one  of  the  dimensionless  constants  which  define  the  boundary- 
condition  at  the  shock  front.  Solution  (2)  for  a  =  0  and  Y=  7/5  describes 
shock  wave  propagation  in  the  absence  of  radiation  (q  =  0).  Solution  (2) 
for  a>0,  i.  e. ,  q>0,  describes  propagation  of  a  shock  wave  with  radiant 
energy  absorption.  A  numerical  solution  of  (1)  with  boundary  conditions 
for  a<0,  i.  e.  ,  q<0,  describes  propagation  of  a  shock  wave  front  with 
self-radiation. 


Dremin,  A.  N,  ,  K.  K.  Shvedov,  and  V.  S. 

Klapovskiy.  Effect  of  charge  shape  on  detonation 
parameters  of  PZhV-20  ammonite.  IN: 

Detonatsiya  vzryvchatykh  veshchestv  i  bezopasnost' 
vzryvnykh  rabot.  Moskva,  Izd-vo  Nedra,  1967, 

129-135. 

This  work  presents  some  experimental  materials  concerning 
the  study  of  the  detonation  parameters  of  spherical,  cylindrical  and  flat 
charges.  The  experiments  were  performed  with  charges  of  various  diameters 
and  thicknesses.  The  explosive  used  was  PZhV-20  ammonite,  with  a 
density  of  1.0  g/crn^  in  all  experiments.  The  experimental  data  showed 
that  with  identical  characteristic  dimensions  of  the  charge,  the  shortest 
wavelength  and  sharpest  drop  was  observed  in  spherical  charges.  The  time 
and  wavelength  in  flat  charges  is  approximately  30%  greater  than  in 
cylindrical  charges.  Thus,  a  flat  charge  with  a  thickness  equal  to  the 
diameter  of  a  cylindrical  charge  or  with  equal  weight  will  be  more  effective 
from  the  point  of  view  of  generation  of  shock  waves  of  higher  intensity  and 
longer  duration  in  the  surrounding  medium. 


Afanasenkov,  A.  N.  ,  I.  M.  Voskoboynikov,  and 
A.  Ya,  Apin.  Transmission  of  detonation  through 
an  air  gap.  IN:  Detonatsiya  vzryvchatykh  veshchestv 
i  bezopasnost'  vzryvnykh  rabot.  Moskva.  Izd-vo 
Nedra,  1967,  101-110. 

The  purpose  of  this  work  was  to  produce  a  simple  criterion 
for  explosive  sensitivity  and  to  study  the  possible  mechanism  of  transmission 
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of  detonation  through  an  air  gap  on  the  basis  of  available  data.  Experiments 
were  performed  using  passive  charges  20  mm  in  diameter  and  70  mm 
long  of  hexogene  of  various  densities,  and  active  charges  of  hexogene  and 
TG  50/50  of  various  weights  and  shapes.  The  experiments  determined  the 
significance  of  explosion  products  and  the  air  shock  wave,  and  produced  a 
quantitative  characteristic  of  the  sensitivity  of  explosives  to  detonation 
transmission  through  an  air  gap.  The  dependence  of  the  distance  of  trans¬ 
mission  of  detonation  on  explosive  charge  density  was  explained. 


Khristoforov,  B.  D.  Small-scale  electromagnet 
for  explosion  investigations.  PTE,  no.  1,  1972, 

251. 

An  electromagnet  has  been  built  for  laboratory  measurements 
of  the  mass  and  wave  velocities  of  shock  and  detonation  waves.  The  field  of 
intensity  across  a  200  x  200  x  200  mm3  gap  is  ~  400  Oe;  in  addition  there 
is  a  permanent  magnetic  field  which  extends  over  a  50  x  50  x  40  mm3  region 
near  the  gap  center.  The  70  mm  diameter  magnetic  circuit  is  made  of  900  x 
900  mm^  Arrmco  iron.  The  corners  of  the  device  are  strengthened  by  welding 
and  special  straps  to  increase  resistance  to  explosive  loads.  Dimensions  of 
the  magnet  gap  and  the  permanent  magnetic  field  are  controlled  by  specially 
designed  pole  shoes  made  of  a  200  mm  dia  type  St.  3  steel.  The  field 
intensity  across  the  gap  can  be  increased  in  inverse  ratio  to  the  end-area 
of  a  tapered  attachment  to  the  pole  shoe.  Two  paralleled  magnetic  coils, 
each  composed  of  2000  turns  of  1.  8  mm  copper  wire,  are  used  for  superposed 
magnetization.  The  24  volt  current  supply  is  supplied  by  a  VSA6  rectifier. 
The  magnetic  circuit  is  installed  in  a  channel  welded  from  channel  bars;  the 
channel  itself  is  welded  to  a  table  made  of  angle  steel  to  increase  the  rigidity 
of  the  assembly.  Over  the  past  three  years,  the  electromagnet  has  withstood 
several  hundred  detonations  of  explosive  charges  up  to  0.5  kg  without  sub¬ 
stantial  damage.  A  photograph  of  the  electromagnet  assembly  is  included. 


Simonov,  I.  V.  Diffraction  of  a  strong  shock 
wave  on  a  weakly-defined  wedge.  ZhPMTF, 
no.  6,  1971,  107-114. 

A  theoretical  study  is  presented  on  diffraction  of  a  strong 
two-dimensional  steady  shock  wave  on  a  wedge  whose  properties  (initial 
density,  behavior  under  impact  load)  differ  only  slightly  from  those 
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of  the  medium  in  which  the  shock  wave  propagates.  A  linear  approximation 
of  diffraction  -  initiated  small  perturbations,  disregarding  differences  in 
strength  and  sound  velocity  c  in  the  medium  and  wedge,  is  justified  under 
conditions  of  high  pressure  (~10^  atm)  and  temperature  behind  the  strong 
shock  wavefront.  The  validity  of  the  approximation  is  supported  by  the  fact 
that  the  effect  of  strength  characteristics  diminishes  with  increased  shock 
wave  intensity.  In  this  approximation,  the  unknown  perturbation  velocity 
components  u1  and  w'  and  perturbation  pressure  p'  satisfy  the  ordinary 
linear  equations  for  the  two-dimensional  flow  of  an  ideal  compressible  fluid. 
The  equations  are  given  in  the  form 


(i) 


where  u,  w,  and  p  are  dimensionless  variables.  Conditions  in  the  pertur¬ 
bation  region  of  the  shock  wavefront  are  described  by 


P  —  A  (J  —  yf '),  u  =  Bp,  w  =  —Mf  at  x  -k 


(2) 


in  the  medium,  and 
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in  the  wedge.  Vo  and  Uo  in  (3)  are  specific  volume  and  mass  velocity, 
respectively,  behind  the  shock  wavefront,  and  V  is  the  specific  volume  of 
the  medium  ahead  of  the  front.  The  problem  was  reduced  to  determining  the 
function  p(x, y),  which  is  continuous  along  the  discontinuity  boundary  line. 
Two  possible  diffraction  wave  configurations  are  examined.  At  wedge  angles 
al»  a2  <  a  *»  t*le  I!1™*  anglei  a  triple  configuration  (Fig.  1)  forms  near  the 
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intersections  L  and  F  of  the  shock  wave  with  the  edges  of  a  wedge.  The 
nonstationary  perturbation  from  the  shifted  wedge  center  O  propagates 
within  the  region  hounded  by  the  ABCD  arc  of  the  Mach  circle,  the  attached 
reflected  acoustic  wavefront;,  and  the  AD  portion  of  the  shock  wave.  At 
0  1,2  >  a  *>  there  is  an  irregular  refraction  near  L  and  F  and  only  one  region 
of  non  stationary  flow,  which  includes  the  shock  wave  sections  outside  the 
wedge  (Fig.  2). 


A 


Fig.  2.  Diffraction  pattern  at  aj,  a  ^  >  a  *. 
The  a*  value  in  the  system  of  x'y'  coordinates  is  determined  by 


a*  =  arc  sin  ( fc0  /  +  A0*)  (*i  =  Vi  -  kl) 

where  kQ  =  D0/C0,  DQ  =  shock  wave  velocity. 


(4) 


In  any  case,  the  problem  of  determining  p(x, y)  is  transferred 
into  the  plane  Z  =  xj  +  iyj  (Fig.  3)  and  the  luneis  represented  conformally 


Fig.  3.  Region  of  nonstationary  flow 
in  the  Z  complex  plane. 


on  the  upper  half-plane  £  by  means  of  the  transform 


t  =  l  +  ir\  (5) 

The  analytical  function  F*(£)  of  £  half-plane  satisfies  the  condition 

j-hdP- 


where 


a^  +  b^  =  d 


a  =  £/l-£a,  b  =  (K I  <  1) 

a  =  0,  b=B-Bi  (|5|>1) 

d  =  (B-  Bt)  (p,6  (s  -  l,)  -  p26  a  -  S,))  (a,  ,2  <  a*) 

d  =  (^  +  ^  A)  (6  (6  -  l3)  -  6  (s  -  £,))  >  a,) 

£i  =  Ib  =  V1  (cosec  9j  —  A-  ctg  OJ,  =  £c  =  «r*  (cosec  02  —  A  clg  02) 
£s  =  II—  —  x_1clgai,  £4  =  If  =  x_1  ctg  x2 


(6) 


(7)  and 


8  (0)  is  a  Dirac  function. 


The  formulas  (6)  and  (7)  represent  the  inhomogeneous  Hilbert 
problem  with  continuous  parameters.  The  problem  is  solved  by  determining 
the  function  having  a  zero  point  not  lower  than  second  order  at  infinity.  The 
solution  must  satisfy  the  condition 

0<A*(l.-/)<xV  (1  +/)  (8) 


which  is  equivalent  to  the  condition  of  stability  of  a  two-dimensional  stationary 
shock  wave  in  a  homogeneous  medium.  The  boundary  4>+(£)  of  the  analytical 
function 


<S>+(£)  = 


Vi-f- 


(9) 


satisfies  the  condition  (8),  since  it  does  not  exhibit  singularities  in  the 
upper  half- plane.  The  function  F+(£)/  (£)  must  be  regular  at  infinity. 

It  follows  that 


F*(£)  =  O'  (£)(TMC)  +  Col  (10) 
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The  constant  CQ  in  (10)  is  determined  from  the  condition  of  smooth 
conjugation  of  the  wavefront. 

Pressure  is  determined  from 

p(g,T|)=  Im  J  (2jjj-+  id-Hr)dt  + pA  (I2) 

where  PA  =  Pj  for  a  j  ^  2  <  a  *  and  pA  =  0  f or  a  j  2  >  a  *. 

Pressure,  as  well  as  f',  u,  and  w,  exhibit  a  logarithmic 
singularity  in  L  and  F.  A  formula  is  also  derived  for  p  distribution  p(y) 
along  the  curved  section  AD  of  the  shock  wave.  The  functions  u  and  w  are 
determined  from  (1)  and  expressed  in  terms  of  p  in  a  closed  system  of 
equations. 


It  is  noted  that  the  results  obtained  may  contribute  to  the 
development  of  a  theory  of  irregular  refraction. 


Salamandra,  G.  D.  and  N,  M.  Ventsel'. 
Measurement  of  gas  velocity  with  the  use  oi 
the  contact  discontinuity  surface.  ZhTF, 
no.  11,  1971,  2463-2465. 


,  A  method  is  described  for  gas  velocity  measurements  behind 

a  detonation  wave  front.  The  contact  dis  continuity  surface  formed  in  the 
gas  is  used  as  the  point  of  reference  for  gas  motion;  the  contact  discontinuity 
surface  is  assumed  to  be  formed  by  the  interaction  between  shock  and 
detonation  waves.  The  surface  is  visualized  by  the  Toepler  technique.  For 
several  experiments,  the  numerical  value  of  gas  velocity  measured  16%  lower 
an  the  theoretical  value.  Shock  tube  experiments  show  that  the  contact 
is  continuity  surface  can  give  a  reasonably  accurate  value  of  shock  wave 
velocity,  and  can  furthermore  serve  as  a  reference  for  evaluating  other 
velocity  measurement  references,  e.  g.  the  heat  inhomogeneity  technique. 


Zaydel1,  R.  M.  Shock  wave  passage 
through  a  curved  boundary  interface 
between  two  media.  MZhiG,  nw.  1,  1972, 

111-121. 

Shock  wave  interaction  with  a  slightly  curved  boundary  surface 
and  shock  wavu-induced  surface  motion  are  analyzed  in  an  x-y  coordinate 
system.  The  shock  wave  theoretically  propagates  with  velocity  D  through 
a  light  medium  (X  >  0)  toward  the  boundary  surface  X  =  0  and  interacts  with 
it  during  the  time  t  =  2  a0/D0  (Fig.  1).  During  the  interval  t,  the  A  point 


Fig.  1.  Shock  wave  interaction  with 
the  free- stream  interface  X  =  0  between 
two  media. 

shifts  to  A'  and  the  points  between  A  and  A'  are  distributed  at  a  time  r  along 
the  boundary  6 

«  (//.  }  )  (jj 

where  kaQ  S  1  is  a  small  perturbation.  The  free- stream  boundary  is  set  in 
motion  after  T  at  a  velocity  U  relative  to  an  uncompressed  heavy  medium 
(X<0).  Vj  and  V2  are  the  relative  velocities  of  the  reflected  and  transmitted 
shock  waves.  A  set  of  gas  dynamic  equations  with  boundary  and  initial 


conditions  is  formulated  for  small  perturbations.  Introducing  polar 
hyperbolic  variables 

y,  —  r,  ch  0„  x,  =  r,  sh  8„  r,  =  Vy.J  —  x,\  th-0,  —x,/y,  (2 ) 

and  using  a  Laplace  transform,  the  author  obtains  a  solution  to  the  equations 
in  the  general  form 

w."(q„  9*)  =*’.(?. +  0.)  f  0.(ar.  —  8.)  (s  =  1,  2)  (3) 

where  F  and$are  arbitrary  functions,  ps  =  sinhgs,  s  =  1  in  the  0<X  <Vjt 
region,  and  s  =  2  in  the  -V^t  <0  region.  Integration  of  the  original  gas 
dynamic  equation  in  polar  hyperbolic  coordinates,  using  (3),  three 

finite  difference  equations  which  have  a  unique  solution  in  the  form  of  expansion 
in  an  exponential  series.  In  the  Ap — ►O  approximation,  where  A.p  is  the 
shock  wave  intensity,  the  solution  is 

<I>t(2«)  =  ~'l2Vt°e-\  Ft(qi)  =  xkvl°e~^ 

(4) 


u(t)  =  ka0Uty(T)  (t  =  kcxt) 

|  t 

Jl|5.(i)d*  +  ^(T)  +  tJ),(T)  (5) 

0 

where 

V(p)  =  p-i'T,(p)  +^(p)  +  'P,(p) 

M',(p)  =  (l-  n)[p/r  +  l/J?]-\  ^Fa(p)  (M  —  l/«)  [fr-h  r/JF?]- 

_  ; _  (6) 

Y,(p)  =  (l-a)[l  +  pfl/r]-‘,  r  =  ypJ  +  l,  i?  =  Valp’  +  l 


is  the  transform  of  ^(p),  cj  is  the  sound  velocity  in  the  s  =  1  region. yU  -  p-y/p^, 
and  a=  ci/c2»  Expression  (6)  is  simplified  in  the  cases  of  a  =  1  (ci  =  C2), 

Pl  =P2  {/*■  =  1).  or  the  acoustic  impedance  pj,  cj  =  P  2CZ  (  °/a=  1). 

Integration  of  (5)  gives  the  asymptotic  value 


00)  =  const  =  — — —  =  — — — 
1+P  P:  +  P» 


(7) 


i.  e.  'J'(t)  depends  on  thepj/p2  ratio  only.  A  generalization  of  the  conclusion 
led  to  an  equation  of  motion  for  the  free- stream  boundary 


#a 
dt * 


P»  —  P» 
P»  +  P. 


kga 


(8) 


where  tne  acceleration  g  =  NU,  and  N  is  the  number  of  incident  shock  waves, 
each  imparting  AU  acceleration  to  the  boundary  interface.  This  equation 
is  identical  to  the  equation  of  Rayleigh- Taylor  for  gravitational  instability 
of  noncompressible  fluids.  It  is  concluded  that  the  effect  of  compressibility 
is  insignificant  when  the  wavelength  of  boundary  perturbations  is  sufficiently 
small. 
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Bezhanov,  K.  A.  Diffraction  of  a  shock 
wave  on  a  wedge  moving  at  supersonic 
velocity.  Prikladnaya  matematika  i 
mekhanika,  no.  4,  1969,  631-637. 

A  theoretical  study  is  made  of  the  diffraction  of  an  arbitrarily 
strcng  shock  wave  on  the  surface  of  a  wedge,  moving  with  supersonic  velocity. 
Using  the  notations  of  Fig.  1,  the  author  assumes  that  the  incident  angle  8  of 


Fig.  1.  Flow  model. 


the  shock  wave  is  approximately  the  same  as  the  wedge  angle  a.  In  this 
notation  UQ  defines  absolute  shock  wave  velocity,  while  U  is  its  velocity 
relative  to  adjacent  flow.  With  the  use  of  some  simplifying  assumptions 
the  analysis  is  made  with  a  linearized  equation  for  planar  gas  flow,  using 
dimensionless  parameters.  The  work  is  an  extension  of  some  cited  earlier 
treatments  of  a  similar  problem  in  which  diffraction  of  the  shock  wave  was 
assumed  not  to  take  place. 


Borovoy,  V.  Ya.  and  M.  V.  Ryzhkova. 

Heat  transfer  on  a  semi-conic  surface  at 
high  hypersonic  velocities.  MZhiG,  no.  4, 

1969,  137-142. 

Formation  of  narrow  bands  of  intensive  heat  transfer  on  the 
flat  surface  near  the  edges  of  a  streamlined  spherically-blunted  semi-cone 
(Fig.  1)  was  studied  experimentally  at  a  Reynolds  number  Re  in  the  1.  5  x  10^ 
-  2.  4  x  106  range.  At  a  sufficiently  high  Re,  the  boundary  layer  separates 
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Fig.  1.  Spherically  blunted  semi-cone. 

from  the  flat  surface  and  two  separate  vortex  formations  appear  near  the 
edges  as  boundary  lines  along  the  surface.  Most  of  the  heating  occurs 
along  the  boundary  lines  where  heat  transfer  coefficient  h°  may  be  several 
times  higher  than  h°  on  the  remaining  part  of  the  flat  surface.  The 
experiments  show  that  Fe  is  the  most  important  factor  in  formation  of  the 
boundary  lines  of  intensive  heat  transfer.  At  Re  =  2.  5  x  104  -  10  x  104, 
hmax  increafes  nearly  proportionally  to  Rf},  while  at  Rq  >  10  x  104  the  rate 
of  h^ax  increase  is  somewhat  slower.  The  strong  effect  of  Re  confirmed 
the  assumption  that  the  laminar  flow  changes  to  turbulent  at  the  boundary  of 
the  separation  zone.  Boundary  layer  separation  occurs  because  of  an  over¬ 
flow  of  gas  from  the  conic  to  the  flat  surface.  This  was  confirmed  by  an 
experiment  with  a  barrier.  Experiments  with  a  reduced  thickness  boundary 
layer  and  increased  conic  surface  roughness  did  not  contradict  the  proposed 
flow  pattern. 


Myshenkov,  V.  I.  Subsonic  and  transonic 
viscous  gas  flow  in  the  wake  of  a  two-dimensional 
body.  MZhiG,  no.  2,  1970,  73-79. 

Viscous  gas  flow  in  the  wake  of  a  finite  thickness  plate  is 
analyzed  at  subsonic  and  transonic  velocities,  Prandtl  number  =  0.71,  and 
Reynolds  numbers  R  in  the  1-1,  000  range.  The  problem  is  formulated  by 
Navier-Stokes  equations  for  given  flow  parameters  at  infinity  and  at  a 
finite  distance  from  the  bottom  of  the  plate.  Equations  of  state  and  motion 
are  solved  by  approximation  of  a  finite-difference  scheme  in  two  steps. 

A  numerical  solution  which  satisfies  the  condition  of  a  steady  flow  is  presented 
graphically  for  different  R  and  M  values.  The  solution  shows  that  at  RS  1,7 
gas  flows  without  separation  from  the  bottom,  and  at  R>1.  7  a  region  of  reverse 
flow  with  separation  forms  at  the  stagnation  point.  The  reverse  flow  region 
expands  with  increase  in  R  and  the  separation  point  shifts  upwards  in  the 
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direction  of  the  angular  point  of  the  plate.  A  deformation  of  the  reverse 
flow  region  with  a  decrease  in  relative  velocity  of  reverse  current  is 
observed  at  R  =  100  at  increased  M  (from  0.288  to  1.15).  The  flow  patterns 
are  in  agreement  with  earlier  theoretical  and  experimental  data.  Pressure 
variation  patterns  at  the  bottom  of  the  plate  are  also  shown  in  the  flow 
plane  and  the  plane  normal  to  flow. 


Stulov,  V,  P.  ,  and  L.  I.  Turchak.  Non¬ 
equilibrium  chemical  reactions  in  a  shock 
layer  in  a  streamline  flow  of  a  carbon  dioxide 
-nitrogen-argon  mixture  around  a  sphere. 

MZhiG,  no.  5,  1969,  147-150. 

A  theoretical  evaluation  of  hypersonic  (M  =  16-20)  flow  of 
C02~N2“Ar  mixtures  around  a  sphere  was  made  using  a  known  physico-chemical 
model  of  a  shock-heated  gas  mixture.  The  model  presumes  the  existence 
of  nine  components  in  the  shock  layer,  with  allowance  for  five  molecular 
dissociation  and  five  atom  exchange  reactions  which  are  tabulated.  Dissociation 
rate  constants  were  calculated  for  collisions  with  particles  of  different 
catalytic  activity.  The  presence  of  O2  was  disregarded  in  the  establishment  of 
the  model,  and  molecular  vibrations  at  T  4000°K  were  considered  to  be  at 
equilibrium.  It  was  shown  that  the  flow  changes  significantly  with  a  change 
in  relative  concentrations  Yoo  (CO2)  and  Y*,  (N2)  in  the  shock  layer.  Shock 
wave  separation  increases  by~50%  when  YM  (CO2)  is  decreased  from  0.95 
to  0.  05,  and  thermal  energy  absorbed  by  the  gas  simultaneously  increases 
because  of  a  decrease  in  T  and  an  increase  in  the  fraction  of  atomic  components. 
The  effect  of  Ar  addition  is  analogous  to  that  of  N2.  Most  of  the  heat  is 
absorbed  by  GO2  dissociation  because  of  its  higher  rate. 


Shifrin,  E.  G.  Formation  of  a  "hanging" 
compression  shock  in  the  streamline  flow 
past  a  profile  with  a  broken  generator.  PMM, 
no.  6,  1970,  1159-1167. 

Uniform  supersonic  streamline  flow  past  a  profile  with  a 
convex  salient  point  is  analyzed  in  the  case  of  an  oncoming  stream  with  a 
low  supersonic  velocity,  when  entropy  changes  in  the  shock  wave  can  be 
disregarded.  The  region  behind  the  shock  wave  is  projected  on  the  hodograph 
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Fig.  1.  Projection  of  the  flow  behind 
shock  wave  on  hodograph  plane  T]  -  f}  . 


plane  q  -fj  (Fig.  1)  where  =  (k  +  1)1 /3(  \  _  l)f  /3  is  the  inclination  of  the 
velocity  vector  to  the  symmetry  axis,  K  is  the  velocity  factor,  and  k  is 
the  adiabatic  exponent.  A  streamline  flow  with  attached  shock  wave  past 
a  profile  formed  by  straight  line  intercepts  is  considered  in  a  first  approxi¬ 
mation.  The  equations  of  the  shock  polar  /9  and  intersecting  characteristics, 
e.  g.  ,  aja2,  are  given  in  a  transonic  approximation.  Proof  is  offered  that 
a  single  "hanging"  shock  is  necessarily  generated  behind  the  shock  wave. 

The  "hanging"  shock  propagates  downstream  from  its  starting  or  end  point 
which  is  located  at  an  infinite  distance  from  the  profile.  The  shock  at  0 
is  convex  in  relation  to  the  flow  behind  it.  In  addition  to  the  tail  shock,  there 
is  an  internal  shock,  at  fiQ<  0,  which  intersects  the  tail  shock  and  forms 
a  lambda- shaped  shock. 


Apshteyn,  E.  Z.  ,  L.  G.  Yefimova,  and 
G.  A.  Tirsky,  Intensive  destruction  of  a 
vitreous  body  by  Cshockj  radiation.  MZhiG, 
no.  2,  1971,  131-134. 

Destruction  of  a  vitreous  body  by  radiation  from  a  shock  layer 
in  hypersonic  gas  flow  is  analyzed  assuming  an  optically  thin  boundary  layer 
at  high  injection  rates  and  sufficiently  intense  vaporization.  Presence  of  a 
liquid  film  near  the  stagnation  point  on  the  body  surface  is  presumed,  to 
simplify  calculations  of  destruction  parameters.  Application  of  an  approximation 
of  radiative  heat  conduction  in  the  body  to  the  liquid  film  led  to  a  set  of 
transcendental  equations  which  describe  ablation,  friction,  and  surface 
temperature  as  functions  of  shock  layer  parameters.  Iterative  solution  of  the 


equations  is  illustrated  by  diagrams  of  the  destruction  rate  and  gasification 
r  of  a  heat  insulating  material  with  thermophysical  characteristics  of  tused 
quartz  under  assumed  flight  conditions.  Plots  of  Tversus  body  radius 
Rg  show  that  Tdecreases  with  a  decrease  in  Rg  due  to  liquid  phase  ablation. 

At  great  Rg  values,  coating  destruction  proceeds  by  vaporization  only.  The 
contribution  of  liquid  phase  ablation  to  the  destruction  rate  becomes  significant 
at  small  Rg,  because  of  the  effect  of  the  radiative  heat  conductance  parameter 
\  H  (Fig.  1).  Since  it  increases  greatly  with  increase  in  atmospheric  reentry 


Fig.  1.  Destruction  rate  of  heat 
insulating  coating.  Curves  1,2,  3,  4 
correspond  to  Ag  =  0,  104,  5  x  104, 
105,  respecti’  .y. 


velocity,  radiative  heat  flux  is  thus  a  decisive  factor  in  calculation  of 
destruction.  A  previous  study  on  this  problem  was  reported  by  Apshteyn 
in  Effects  of  Strong  Explosions,  no.  1,  1971,  p.  22. 


Antonets,  A.  V.  Computation  of  three 
dimensional  BuperBonic  flow  past  blunt  bodies 
with  generator  breaks  and  allowance  for 
equilibrium  and  frozen  states  of  the  gas  in  the 
shock  layer.  MZhiG,  no.  2,  1970,  178-181. 

An  algorithm  is  developed  for  computation  of  the  steady  super¬ 
sonic  flow  of  an  inviscid  gas  past  a  blunt  body  with  multiple  breads  of  a 
generator.  The  algorithm  is  based  on  the  principle  of  separation  (behind 
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each  of  an  arbitrary  number  of  breaks)  of  the  regions  of  smooth  solutions 
between  the  body  and  the  character istic  shock  wave  surfaces,  at  a  positive 
(expansion)  or  negative  (compression)  break.  Finite-difference  equations 
in  the  nodal  points  of  a  net  are  substituted  for  differential  equations  in  the 
regions  of  smooth  solution.  The  unknown  velocity  components,  pressure  p, 
and  enthalpy  i  are  determined  from  the  equations  by  the  method  of  successive 
calculations.  The  algorithm  accounts  for  the  effect  of  nonequilibrium 
physicochemical  reactions  by  computing  the  flow  for  either  the  entire  shock 
layer  in  a  nonequilibrium  (frozen)  state  or  for  only  a  partially  frozen  layer 
in  the  area  adjacent  to  the  body  surface,  the  other  part  being  occupied  by  a 
gas  at  equilibrium.  The  effect  of  the  gas  equation  of  state  in  the  shock  layer 
is  illustrated  by  p(x)  plots  for  the  flow  past  a  spherically  blunted  cone  with  an 
apex  half  angle  0  £  10°.  The  plots  show  that  at  9  =  10°  p  in  some  points  on 
the  surface  in  a  flow  at  equilibrium  can  be  twiae  that  in  the  frozen  flow.  The 
computation  method  may  be  useful  in  evaluating  actual  flight  effects  at  altitudes 
above  60  km  and  at  velocities  of  6,  500  m/sec  or  greater. 


Antonets,  A.  V.  ,  A.  V.  Krasil'nikov,  and 
V.  I.  Lagutin.  Experimental  determination  of 
the  position  of  pressure  center  in  hypersonic 
gas  flow  around  a  blunted  cone  at  different 
angles  of  attack.  MZhiG,  no.  2,  1971,  142-143. 

Following  from  the  previous  paper,  the  relationship  between 
the  position  of  a  pressure  center  and  an  angle  of  attack  or  was  determined 
in  a  hypersonic  wind  tunnel  at  =  6  and  15  and  Reynolds  number  Re  S£  10^ 
using  a  spherically  blunted  test  cone.  The  cone  had  a  sphere  radius  R, 
relative  conic  length  x/R  =  5,  and  apex  half-angle  9  =  11°.  The  highly  accurate 
method  of  free  aerodynamic  balance  was  applied  while  varying  the  distance 
xQ  of  the  rotational  axis  from  the  cone  nose.  The  value  for  a  balanced  a 
coincides  with  xQ.  Plotted  experimental  Cj  data  (Fig.  1)  agree  within  a  2% 


Fig.  1.  Pressure  center  versus  angle 
of  attack:  curves  1  and  2  -  calculated  for 
Mao  =  6  and  15,  respectively;  small  circles 
-  experimental  data  for  Mco  =  6;  points  - 
experimental  data  for  M*  =  15. 
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maximum  discrepancy  with  theoretical  data  calculated  by  the  method  of 
Antonets  cited  in  the  previous  paper  [MZhiG,  no.  2,  1970,  178],  but 
disagree  with  the  (  a  )  data  calculated  by  the  Newton  formula  (dotted  line). 
Results  suggest  that  physical  phenomena  on  the  leeward  side  of  a  cone  at 
a  5  6  do  not  significantly  affect  aerodynamic  characteristics. 


Voronkin,  V.  G.  Nonequilibrium  viscous 

flow  of  a  multicomponent  gas  near  the  stagnation 

point  of  a  blunt  body.  MZhiG,  no.  2,  1971,  144-147. 

A  theoretical  study  is  made  of  the  effect  of  shock  layer  variable 
parameters  on  heat  transfer,  in  the  case  of  nonequilibrium-dissociated 
air  at  a  barrier  of  finite  catalytic  activity  near  the  stagnation  point  of  a 
spherically  blunted  body.  Heat  transfer  in  the  nonequilibrium  outer  flow  is 
compared  to  that  in  an  equilibrium-dissociated  gas  flow.  Numerical 
calculation  of  the  nonequilibrium  flow  was  based  on  the  concept  of  a  compressed 
viscous  shock  layer.  Only  free-stream  Reynolds  numbers  R  above  1,500 
were  considered.  Under  these  assumptions,  the  calculation  was  reduced 
to  a  solution  of  the  boundary  layer  equations  with  boundary  conditions  corre¬ 
sponding  to  those  immediately  behind  the  leading  shock  wave.  Temperature 
profiles,  concentrations  of  oxygen  and  nitrogen  atoms  in  the  shock  layer, 
and  the  relative  heat  flux  Qw  were  calculated  at  a  60  km  height,  7.  4  km/sec 
free-stream  velocity,  700°K  wall  temperature  and  1-100  cm  radius  of 
curvature  R.  It  was  shown  that  Qw  in  the  nonequilibrium  shock  layer  at 
small  R  can  be  significantly  increased  (by  a  factor  of  3)  in  relation  to  the 
theoretical  Qw  in  an  asymptotically  thin  nonequilibrium  boundary  layer  with 
equilibrium  parameters  at  its  outer  boundary.  At  increased  R  the  vorticity 
effect  decreases  and  the  effects  of  nonequilibrium  chemical  reactions  become 
predominant. 


Lunev,  V.  V.  ,  and  N.  Ye.  Khramov.  Flow 
past  a  blunt  body  near  the  stagnation  point  in 
a  divergent  hypersonic  stream.  MZhiG,  no.  3, 

1970,  102-105. 

The  inhomogeneous  hypersonic  outer  stream  near  the  stagnation 
point  past  a  blunt  smooth  body  is  analyzed  at  a  small  ratio  k  of  densities  ahead 
of  and  behind  a  straight  shock  wave.  Thickness  8  of  the  shock  layer  at  the 
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forward  stagnation  point  is  formulated  by 


S  _  1  _2_ 

kR *  2  U  k 

_6 _ 1 

a/?*-  i  +  y5iE 


v  =  0 


(v  =  l) 


(1) 


where  R*  is  the  effective  curvature  radius  of  the  body  and  v  =  0  and  v  =  1 
refer  to  two-dimensional  and  axisymmetric  flows.  The  formulas  in  (1) 
were  derived  from  equations  of  motion  assuming  gas  density  to  be  constant. 

It  follows  from  (1)  that  8,  and  hence  the  profile  of  the  shock  layer,  depend 
on  k  and  R*  only.  This  dependence  is  illustrated  by  computed  8(k)  plots 
for  an  axisymmetric,  strongly  underexpanded  conic  jet  around  a  sphere. 

The  initial  flow  parameters  are  tabulated.  The  data  all  lie  on  the  same  curve 
which  coincides  with  the  theoretical  curve  calculated  for  v  =  1  and  R*/R 
boundary  values.  Plots  of  pressure  distribution  p/p^  along  a  barrier  versus 
x/R  ratio  diverge  significantly  for  different  jets,  but  the  same  plots  of  p/p^ 
versus  x/R*  form  a  relatively  narrow  band  which  is  close  to  the  plot  for 
a  sphere  in  a  homogeneous  stream  and  the  Newtonian  plot.  The  effect  of 
stream  divergence  is  thus  accounted  for  by  introduction  of  R*.  The  effects 
of  nonuniform  density  and  velocity  of  the  oncoming  stream  are  taken  into 
account  in  an  analogous  manner. 


Arutyunyan,  G.  M.  Climb  of  a  shock  wave  on 
a  wedge  moving  at  supersonic  velocity.  IN: 

Volny  v  neuprugikh  sredakh  (Waves  in  nonelastic 
media).  Kishinev,  AN  Mold  SSR,  1970,  7-12. 

The  flow  generated  at  the  upper  surface  of  a  moving  wedge  by 
an  oncoming  plane  shock  wave  was  analyzed.  Conditions  are  formulated  for 
the  flow  according  to  the  pattern  shown  in  Fig.  1.  It  was  assumed  that  intensity 
of  the  oncoming  shock  wave  CE  was  equal  to  that  of  the  attached  compression 
shock  OA,  CE  is  normal  to  OB,  and  the  interacting  shocks  CD  and  DA  are 
weak  and  reflected  at  a  $ /2  angle  by  a  rigid  wall  with  resulting  refraction, 
but  without  formation  of  a  third  shock  wave.  It  is  shown  that,  under  these 
assumptions,  the  flow  past  the  wedge  upper  surface  OB  can  be  calculated  using 
linearization  theory,  if  a  point,  defined  by  the  Mach  number  M  of  the  oncoming 
flow  and  the  angle  8  of  the  wedge,  lies  within  certain  boundaries  in  the  M-8 
plane.  All  the  conditions  were  verified  within  an  area  bound  by  the  curves 
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*1  (M»  8  )  -  0,  (M,  8  )»  f3  (M,  8  )»  and  the  abscissa.  If  the  point  defined 


Fig.  1.  Flow  diagram  on  the  upper 
surface  of  a  wedge. 

by  the  wedge  parameters  M  and  5  lies  within  this  area,  the  flow  will  be 
uniform.  Pressure  p3  within  the  area  is  then  given  by 

'£t-hJL£*  L_  f 

p.  vpf^r  (i) 

where  pQ  is  the  pressure  generated  by  interacting  shocks  and  Y  is  the 
adiabatic  exponent. 


Mamadaliyev,  N.  A. ,  and  Kh.  A.  Rakhmatufin. 

Streamline  fluid- solid,  two-phase  flow  past 
a  slender  profile.  IN:  Volny  v  neuprugikh 
sredakh  (Waves  in  nonelastic  media).  Kishinev, 

AN  Mold  SSR,  1970,  146-152. 

A  mathematical  model  of  a  fluid- solid  supersonic  flow  past  a 
slender  profile,  e.  g.  ,  a  wedge,  is  introduced  for  use  in  calculations  of  jet 
parameters  from  supersonic  nozzles,  flight  mechanics  in  a  gas-dust  medium, 
etc.  The  model  makes  allowance  for  post- shock  reflection  of  particles  from 
the  profile  surface.  For  a  specific  volume  <1,  the  effects  of  particles  on  the 


gas  flow  and  friction  between  particles  in  the  region  of  a  three-velocity 
flow  can  be  disregarded.  Solution  of  linearized  equations  of  motion  for 
particles  impinging  on  and  reflected  from  a  wedge  with  a  given  apex  angle 
/3  0  thus  leads  to  the  conclusion  that  the  boundary  between  the  three -velocity 
and  the  two-velocity  flows  is  a  straight  line  with  slope  fi  =  2  fi  o*  The 
boundary  condition  at  y  =  0  is  also  a  straight  line,  and  formulates  the  law 
of  mirror  reflection.  For  the  case  of  a  two-phase  medium  with  a  significant 
particle  concentration,  and  disregarding  the  effect  of  the  reflected  stream, 
a  solution  to  the  inverse  problem  is  presented  in  the  form  of  an  equation  of 
the  profile  surface  for  a  given  flow  separation  boundary  with  a  constant 
inclination  /8  .  In  contrast  to  the  first  analyzed  case,  the  profile  surface  is 
weakly  concave  in  relation  to  the  stream,  and  the  three-velocity  region  II 
narrows  with  the  increase  in  Xj,  as  shown  in  Fig.  1. 


Fig.  1.  1  -  separation  line,  profile  surface; 
2  -  at  a  small  and,  3  -  at  a  significant 
particle  concentration. 


Lunev,  V.  V.  Nonviscous  flow  near  the 
critical  point  of  a  blunt  body  having  a  variable 
density  along  its  axis.  MZhiG,  no.  6,  1971, 

56-61. 

The  flow  in  a  thin  compressed  layer  around  a  smooth  blunt 
body  is  described,  in  the  framework  of  hypersonic  theory.  A  set  of 
conservation  equations  is  derived  using  an  x,  y  curvilinear  coordinate  system 
with  the  stagnation  point  as  origin.  A  localized  closed  solution  near  the 
stagnation  punt  is  possible  when  the  shape  of  the  shock  wave  approaches 
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that  of  the  body  and  the  ratio  k  of  deneitiee  ahead  and  behind  the  shock 
wave  tends  towards  zero.  A  solution  in  the  vicinity  of  the  axis  is  sought 
in  the  form  of  an  expansion  in  xK  powers:  ® 

w==  ~va(y)  u  —  Kxiii (y)  r  =  «  +  ... 

P  =  My)  +  ...,  Ap  —  Po(y)  +  ...  (1) 

where  K  =  R-l  is  the  body  curvature,  uV  and  vV  are  velocity  components 

°"g  X  T  -y  !,XeS  •  °f  the  lncident  flow’  r  is  the  distance  from  the  symmetry 
”  '  ”  P,18  .en8lty*  By  substitution  of  variables,  conservation  equations 

are  derived,  which  are  similar  to  adiabatic  flow  equations  for  a  perfect  eas 
with  nearly  constant  p.  Derivation  of  the  Busemann  formula  from  (1)  g 

entire  fr  t  universf^Y  of  pressure  distribution  can  be  extended  to  the 
ntire  front  surface  of  a  smooth  body.  The  Busemann  formula  also  led  to 

srsswsrsr  rs:.a-  - — * 


<p<consl4-'i-v)^A  (2) 

h'Y  T  solution  in  similar  variables  f.  t.  including  t  8.  depends  only  on 
e  body  values  of  p  distribution  and  the  axial  pressure  gradient  (at  f  -  0) 
and  is  the  same  as  the  solution  in  the  case  of  a  constant  p  =  k  “  alone  thi’ 
thickness  8  of  a  shock  layer.  The  solution  of  (2)  for  an  axisymmetrical  flow 

of  th  d‘°  b;  lTly  depenieat  «  k°-  «1>“«  1  narrow  nr'will  Lgion 

of  the  order  t-k..  .  The  solid  line,  in  Fig.  1  are  plotted  for  the  condition 


Fig.  1.  Axial  velocity  profiles  versus 
initial  coordinate  y  £  ;  the  numbers  1-5 
at  the  curves  correspond  to  those  in  Fig.  2. 
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kQ  =  0.  1.  The  tg  =  »|  g  /k  =  0.  80  value  calculated  from  the  plots  in 

Fig.  2  in  effect  coincides  with  the  8  /kR  —  0.78  value,  which  was  obtained 

earlier  for  a  streamline  flow  of  dissociated  air  at  an  equilibrium  within 

k  =  0.05  -  0.20  range.  But  the  true  thickness  8 /kR  calculated  for  curves 

1-5  depends  significantly  on  p  distribution.  The  solution  for  a  two-dimensional 


0  0.5  t.O 


Fig.  2.  Density  profiles  versus  initial 
coordinate  y/8  for  an  axisymmetric  flow. 


flow  (v  =  0)  is  seen  to  be  always  dependent  on  kQ;  but  the  dependence  on 
kQ  =  0.05  -  0.2  is  weak  as  shown  in  Fig.  3,  where  the  plots  1-4  correspond 
to  k0  =  0.2,  0.05,  0.01,  and  0.001.  A  certain  degree  of  universality  of  the 
f  (  t)  /  T)  g  )  function  is  consequently  also  to  be  expected  in  the  case  of  two 
-dimensional  flow. 


Fig.  3.  Axial  velocity  profiles  versus 

similar  coordinates  f,  t,  for  a  two-dimensional 

flow. 


According  to  the  second  flow  pattern,  P  distribution  near  the 
wall  exhibits  a  singularity.  A  significant  variation  occurs  in  the  p  profile 
which  is  bound  to  affect  the  flow,  particularly  in  the  near-wall  layer.  Even 
in  this  case,  however,  the  f  (  *1  /  *1  g  )  profile  is  expected  to  remain  universal, 
because  of  the  weak  dependence  of  the  solution  on  the  parameter  kQ,  i.e. , 
the  gas  p  near  the  wall. 


It  is  concluded  that  in  all  except  the  near-wall  sublayer 
case,  the  solution  should  be  the  same  as  for  a  constant  p  equal  to  the  p 
value  behind  the  shock  wave.  Calculation  of  the  t  tructure  of  an  inviscid 
shock  layer  can  thus  be  simplified,  given  the  universality  of  the  axial 
flow  profile  in  variables. 


Karasev,  A.  B.  ,  and  T.  V.  Kondranin.  The 
effect  of  ablation  products  on  heat  transfer" 
during  destruction  of  graphite  in  radiative 
air  plasma.  MZhiG,  no.  1,  1971,  23-31. 

A  steady  hypersonic  flow  of  viscous,  heat-conductive,  radiating 
air  around  the  leading  stagnation  point  of  a  graphite  sphere  is  analyzed  in 
the  region  between  the  separated  shock  wave  and  the  body  surface.  The  flow 
is  described  by  a  set  of  equations  with  self- similar  variables  and  boundary 
conditions  behind  the  shock  wave  and  on  the  surface.  Heat  transfer  in  the 
radiant  air  was  calculated  by  solving  the  equations  by  the  method  of  successive 
approximations.  A  binary  model  was  adopted  for  diffusion  of  ablation  products. 
Chemical  reactions  were  considered  to  be  frozen  in  the  boundary  layer  and 
equilibrated  at  the  boundary  layer- shock  wave  interface.  Gar,  state  on  the 
ablation  surface  was  determined  from  the  chemical  equilibrium  condition, 
assuming  equality  of  vapor  pressure  with  stagnation  pressure.  Radiative 
and  convective  heat  transfer  were  calculated  in  the  presence  or  absence  of 
ablation.  Analysis  of  plotted  data  shows  that  the  spectral  range  of  radiant 
energy  transfer  is  narrow  in  the  streamline  flow  around  an  impervious  surface 
at  high  atmospheric  reentry  velocities.  In  the  absence  of  ablation,  the 
total  radiant  heat  flux  changes  insignificantly  by  transfer  through  a  cold 
boundary  layer.  Thus  the  presence  at  the  boundary  layer  of  ablation  products 
capable  of  absorbing  or  emitting  radiation  means  that  heat  flow  to  the  surface 
will  increase  if  the  vapor  product  is  opaque  in  the  visible  range,  but  will 
decrease  if  opacity  extends  to  \  <  0.  115  micron.  Actually,  the  shielding 
effect  at  A<0. 115/*  is  to  a  large  degree  cancelled  by  the  radiative  effect  at 
A  >  0. 115/*.  A  previous  paper  by  these  authors  on  graphite  ablation  was 
published  in  Effects  of  Strong  Explosions,  no.  1,  1971,  p.  41. 


-42- 


Karasev,  A.  B. ,  and  A.  N.  Lyakh.  Study 
of  radiative  and  convective  heat  transfer 
by  an  emissive  mixture  of  carbon  dioxide  and 
nitrogen  flowing  past  the  stagnation  point. 

MZhiG,  no.  2,  1971,  39-47. 

Hypersonic  flow  characteristics  of  a  viscous  thermally 
conductive  mixture  (90%  C02,  10%  N2)  around  a  spherically  blunted  body 
were  calculated  with  allowance  for  radiation  energy  transfer  and  reabsorption. 
Self-similar  variables  and  boundary  condition  equations  describe  the  gas 
flow  in  a  compressed  thin  shock  layer,  assuming  local  thermodynamic 
equilibrium  through  the  entire  shock  layer.  A  solution  was  obtained  by 
iteration  and  successive  computation  for  flights  of  the  spheres  at  9,000  - 
12,  000°K  and  1-10  atm  pressure  behind  the  shock  wave.  The  spheres 
were  of  1  m  rauius  and  had  a  3000°K  wall  temperature.  Analysis  of  resulting 
data  revealed  an  insignificant  radiation  effect  on  convective  heat  transfer. 

The  relative  contribution  of  different  spectral  ranges  to  the  radiative  heat 
flux  on  the  wall  was  found  to  be  approximately  43,  40,  20,  and  15%  for  the 
0.  128  -  0.33^  (CO  band),  0.  33  -  0.  66/Hy'CN  bands),  0.80  -  1.15^  and 
0.  45  -  0. 80/<  ranges,  respectively.  The  contribution  of  the  UV  spectral 
range  decreases  while  that  of  the  visible  range  increases  with  an  increase 
in  pressure.  Radiation  heating  is  five  times  greater  than  convection  heating 
at  10  atm  and  12,OOOOK.  Convective  and  radiative  heat  fluxes  in  the  COrN2 
mixture  exceed  the  corresponding  values  in  air  by  20  and  40%. 


Neyland,  V.  Ya.  Upstream  propagation  of 
perturbations  from  interaction  of  a  hypersonic 
flow  with  a  boundary  layer.  MZhiG,  no.  4,  1970, 
40-49. 


Hypersonic  viscous  gas  flow  past  a  slender  body  is  analyzed  for 
weak,  moderate  and  strong  free  interactions.  In  the  case  of  a  weak  hyper¬ 
sonic  interaction,  perturbations  generated  in  the  vicinity  of  a  singular  point, 

(e.  g.  ,  bottom  cut  or  angular  point),  are  seen  to  propagate  upstream  a  short 
distance  only.  Locally  nonviscous  flow  regions  may  also  form  with  high 
pressure  gradients.  Boundary  layer  separation  into  sub-regions  of  different 
sizes  disappears  in  the  presence  of  a  moderate  or  strong  hypersonic  interaction. 
The  perturbations  in  this  case  propagate  up  to  the  leading  edge  of  the  body  in 
the  process  of  free  interaction  along  the  entire  boundary  layer.  Proof  was 
obtained  that  localized  flow  regions  of  very  high  pressure  gradients  are 
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not  generated  at  the  interaction  parameters  X  £  1.  The  possibility  exists, 
however,  of  a  separation  zone  formation  on  the  order  of  magnitude  of  one 
body  length.  The  hypersonic  boundary  layer  at  X2  1  is  described  in  a  first 
approximation  by  a  set  of  equations  derived  from  Navier-Stokes  equations 
with  boundary  conditions.  Boundary  conditions  at  the  aft  end  of  a  body  must 
be  satisfied  in  addition  to  those  near  the  leading  edge.  A  complete  set  of 
boundary  conditions  was  established  in  the  form  of  a  nonsingular  solution. 
The  conclusion  was  drawn  that  the  flow  at  Mj-*®  and  X  2  1  is  described 
by  the  same  boundary  layer  equations  along  the  entire  separation  zone, 
except  for  small  regions  such  as  those  near  the  angular  points. 


Bogashchenko,  I.  A.  ,  A.  V.  Gurevich,  R.  A. 

Salimov,  and  Yu.  I.  Eydel'man.  Flow  of  a 
rarefied  plasma  around  a  body.  ZhETF, 
v.  59,  no.  5,  1970,  1540-1555. 

Experimental  data  are  compared  with  the  theory  on  the  structure 
of  a  perturbation  zone  in  the  vicinity  of  a  disc  in  a  rarefied  plasma  stream. 
Quasi-neutral  potassium  plasma  was  produced  and  its  streamline  flow  around 
a  0.25  cm  radius  disc  was  studied  in  a  Q-device  (Fig.  1).  Plasma  parameter 
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Fig.  1.  Experimental  Q-device: 

1  -  tungsten  ionizer,  2  -  anode,  3  -  disc  no.  1, 

4  -  disc  no.  2,5-  probe. 

measurements  were  made  using  Langmuir  probes  in  both  the  axial  and  radial 
directions.  Density  and  velocity  of  the  plasma,  temperature  T  of  the 
ionizer,  surface  potential  and  position  of  the  disc,  and  magnetic  field  H 
were  varied.  The  experimental  conditions  in  the  Q-device  effectively 
reproduced  those  in  the  ionosphere  at  ~ 200-1,  000  km  height.  In  the  near 
zone  at  a  distance  Z  ^  1  cm  from  the  disc,  the  self-consistent  electric 
field  effect  on  plasma  distribution  was  predominant.  This  effect  was  mani¬ 
fested  in  an  increased  ion  concentration  near  the  axis.  The  body  trace  had 
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an  oscillatory  structure  in  the  far  zone.  Oscillation  damping  is  determined 
by  the  ionic  temperature  T^j]  which  tends  to  equalize  with  T- .  as  the  result 
of  collisions.  A  change  of  sign  of  the  relative  body  potential  strongly 
affects  the  structure  of  the  perturbation  zone  between  two  discs.  The  experi¬ 
mental  data  agreed  well  with  the  theory  of  neutral  approximation,  except 
in  the  far  zone  at  increased  H  values. 


Lifshits,  Yu.  B.  and  E.  G.  Shifrin.  On  the 
problem  of  transonic  flow  over  a  convex  corner. 

MZhiG,  no.  2,  1971,  67-69. 

Mathematical  analysis  of  two-dimensional  potential  flow  in  the 
vicinity  of  the  sonic  point  O  shows  that  a  similar  solution,  satisfying  the 
conditions  of  Vaglio-Laurin's  solution  [j.  Fluid  Mech.  ,  I960,  v.  9,  no.  l], 
is  possible  for  equations  of  transonic  flow  over  a  convex  corner.  A  solution 
to  the  ^  (u,  v)  hodograph  plane  is  obtained  by  analogy  with  that  of  Tricomi 
equations  which  become  equal  to  zero  at  the  v  =  2/3  u^^  characteristic  of 
the  first  family  of  flow  lines.  A  solution  C  =  |  up/v^  to  the  ^=0  equation 
was  found  in  the  third  quadrant  of  the  u,  v  plane.  It  is  concluded  that  the 
profile  in  the  subsonic  region  must  be  convex.  A  one- sheeted  plot  of  the 
solution  in  the  hodograph  plane  indicates  that  the  sonic  line  is  convey,  in  relation 
to  the  region  of  subsonic  velocities.  The  surface  velocity  vector  at  a  transition 
point  from  the  subsonic  region  is  given  by  u  =  (K^C)1^  <*>2/3  f  where  K  is 
the  curvature  of  profile  in  the  O  point.  The  similar  solution  derived  satisfies 
a  boundary  condition  in  the  vicinity  of  the  salient  point,  as  well  as  the  condition 
of  finite  curvature  in  the  subsonic  region.  The  solution  also  describes  the 
principal  member  of  asymptotic  expansion  in  axisymmetric  and  rotational  flow. 


Filatov,  Ye.  I.  The  optimum  shape  of  a  lifting 
body  at  hypersonic  speeds.  MZhiG,  no.  1,  1972, 

82-86. 

Optimization  of  the  aerodynamic  characteristics  of  a  horizontally 
flying  hypersonic  body  is  attempted,  omitting  some  previous  simplifying 
assumptions.  The  body  is  given  a  rectilinear  coordinate  system  in  which 
directions  of  the  x  and  z  axes  are  opposite  to  the  flight  direction  and  vertically 
downward,  respectively.  It  is  assumed  that  xz  is  the  symmetry  plane;  that 
z  ■■  my  and  z  *  f{x, y)  are  the  upper  and  lower  boundary  planes  of  the  semi-body 
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y>0,  and  that  the  rectilinear  leading  edge  is  formed  by  intersection  of  the 
z  =  my  and  x  =  ky  planes.  It  is  further  assumed  that  the  bottom  section  is 
in  the  x  =  kt  =  !\  plane,  where  t  is  the  semi- span;  the  function  f(x,  y)  can 
be  approximated  by  a  polynomial  of  a  low  order  in  any  x  =  const,  y  =  const 
section,  and  that  the  thickness  of  the  leading  edge  is  negligibly  small  in 
relation  to  the  characteristic  dimension  of  the  optimum  body  in  the  z  direction, 
i.  e.  at  x  =  ky  and  f(x,  y)  =  mz.  The  optimum  body  shape  problem  is  then 
formulated  in  term??  of  a  variational  method.  Solution  of  the  problem  is 
reduced  to  finding,  within  the  class  £  =  £  (f ,  rj )  of  surfaces  satisfying  the 
condition  £  (0,  T) )  =  0,  a  surface  which  would  satisfy  the  maximum  value  of  the 
functional 

K—  (Yi  +  Yt)  (X.  +  Ai)-'  (1) 

where  Xj,  X£  are  the  drag  factors  relative  to  the  double  velocity  head  on 
the  semi-body  and  its  leading  edge,  respectively,  and  Yj,  Y2  are  the 
corresponding  relative  lift  forces.  The  volume  of  the  semi-body  must  be 
V  =  const,  and  the  m,  !,  and  k  parameters  are  assumed  to  vary  with  surface 
geometry.  The  new  coordinates 


6  —  (x—ky)  1 1, 


P 


IF' 


T ]  =  y/l,  t—(z  —  my)H 

^^3,-pk  +  m 


(2) 


are  introduced  in  the  solution.  Using  the  Ritz  method,  the  unknown  optimum 
surface  is  described  by  the  equation 


C -g{6 /**/'+ fl„(ri-0.25)  + 
•f  flu(ri  —  O.i)  -f 

+  (£  —  k  /  2)  +  02i  (£?}  — 

—  ft/30)+...} 


(3) 


The  problem  thus  formulated  contains  two  implicit  parameters:  the  coefficient 
of  friction  Cf  and  the  relative  radius  R  of  leading  edge  curvature.  Numerical 
calculation  of  the  unknown  coefficients  a^j,  a^*  etc*  in  (3)  was  made  by 
the  method  of  descent.  In  the  zero-th  order  approximation,  equation  (3) 
describes  pyramids.  For  a  pyramid  with  a  sharp  leading  edge,  calculations 
show  that  the  larger  the  span  (21),  the  higher  is  the  aerodynamic  quality. 

At  moderate  !\ ,  the  quality  K  of  such  pyramids  and  the  slope  P  of  the  lower 
surface  toward  the  x  axis  only  weakly  depend  on  i\  (Fig.  1  and  2).  At 
moderate  the  friction  T  is  one  third  of  the  total  drag. 
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constant  over  the  entire  range  of  R  in  contrast  to  the  pyramids  with  sharp 
leading  edges.  The  leading  edge  bluntness  does  not  significantly  affect 
the  aerodynamic  quality  of  a  lifting  body.  The  values  of  parameters  m,  / ,  k 
calculated  from  (3)  in  the  first  and  second  approximation  indicate  that  the 
aerodynamic  quality  of  optimum  bodies  differs  by  less  than  0.  1%  from  that 
of  the  corresponding  pyramids.  The  data  are  in  qualitative  and  quantitative 
agreement  with  the  earlier  Western  data  approximating  cone-shaped  bodies 
with  sharp  leading  edges. 

It  is  concluded  that  the  geometry  of  the  body  type  studied, 
having  a  given  volume  and  a  maximum  aerodynamic  quality,  is  nearly 
pyramidal  and  in  good  agreement  with  the  earlier  advanced  hypothesis. 


Temkin,  L.  A.  Approximate  solution  to 
the  problem  of  the  ground  zone  in  a 
rarefied  gas.  MZhiG,  no.  1,  1972,  139-143. 

Two-dimensional  equilibrium  flow  in  a  monatomic  rarefied 
gas  around  a  semi-finite  rectangular  bar  is  analyzed,  on  the  basis  of  the 
Krook  model  for  collision  processes  in  gases.  The  Krook  equation  is 
given  in  the  form 


where  f  =  f(r,u)  is  the  coordinate  and  velocity  distribution  function  of  gas 
molecules,  fQ(r,u)  is  the  local  Maxwellian  distribution  function,  and  V 
is  the  integral  kinetic  operator.  A  set  of  equations 


"(r)  =  JJJ  *7o(r,  u)du, 


T(r)  = 


771 

3kn(r) 


u  <r) = ' r/* (r’ u) 

—  •• 

UJi  u  —  U(r)]*F/0(r,u)du 


(2) 


was  derived  from  (1)  for  gas  density  n,  macroscopic  velocity  U,  and 
temperature  T.  The  set  of  equations  (2),  where  k  is  the  Boltzmann  constant 
and  m  is  the  molecular  mass,  was  solved  by  iteration  in  a  first-order 
approximation  both  with  and  without  allowance  for  conditions  at  infinity. 

Using  the  M-20  computer,  the  macroscopic  flow  parameters  n,  U,  and  T 
were  calculated  and  plotted  against  coordinates  x  and  y,  where  the  x-axis 
is  the  flow  direction  at  infinity.  The  plots  calculated  for  the  Knud  sen 
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number  K*  =  i  and  Mach  number  Moo  =  2,  and  for  K  «  =  1  and  M  oo  =  0.  75, 
dipplayed  a  sharp  rise  in  T  in  the  ground  zone  of  the  bar,  near  the  corners,' 
followed  by  a  smooth  decrease  to  the  T  level  at  infinity.  T  on  the  axis  of 
symmetry  is  maximum  at  the  points  x  «  -  1  and  0,  6  at  M  „  =  2  and  0.  75, 
respectively.  At  both  Moo  values,  n  and  pressure  p=nT  increase  monotonically 
to  the  levels  at  infinity.  At  Moo  =2,  the  stream  line  deviates  from  the 
body  near  the  ground  zone  because  of  decelerating  gas  flow.  Results  are 
in  qualitative  agreement  with  an  earlier  solution  to  the  ground  zone  problem 
obtained  within  the  framework  of  the  Navier-Stolces  equations. 


Shapiro,  Ye.  G.  Shock  layer  radiation  during 
hypersonic  air  flow  around  a  spherical  segment. 

MZhiG,  no.  1,  1972,  101-106. 

A  theoretical  analysis  is  presented  of  the  effect  of  shock 
layer  radiation  on  gas  dynamic  parameters  of  hypersonic  air  flow  (V  n  = 
10-16  km/sec)  around  a  spherical  segment  of  a  base  with  radius  R.  The 
shock  layer  flow  (Fig.  1)  is  described  by  a  set  of  integro-differential 


Fig.  1.  Shock  wave  profiles  and  flow 
line  system  for  a  segment  of  R  =  1  m  at 
V  eo  =15  km/sec:  the  flow  line  numerals 
I  thru  V  correspond  to  iJ  =  0.  1, 0.  3,  0.  5, 
0.7  and  0.9,  respectively. 
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equations  in  the  spherical  coordinate  system  r,  d  ,  whose  origin  is  shifted 
at  a  distance  /  =  0.  65  R  from  the  segment  center.  The  term  Q  or  the 
divergence  of  the  integral  radiation  flux  in  the  equations  was  apparently 
known.  The  equations  were  then  solved  for  gas  dynamic  parameters 
using  an  iterative  method  and  an  independent  variable 

1=  (r  —  rb)(r.  —  rb)~l  (1) 

where  rs  is  the  shock  wave  profile  and  r^  is  the  body  profile.  All  gas 
dynamic  parameters  and  rs  were  approximated  by  Lagrangian  polynomials 
with  Hugoniot  relations  and  impermeability  as  boundary  conditions  in  the 
shock  wave  and  on  the  body  surface,  respectively. 

Calculations  were  made  for  segments  of  R  =  0-4  m  over  the 
cited  range  of  V oo  .  In  most  cases,  temperature  Too  ahead  of  the  shock 
wave  was  assumed  to  be  250°,  and  in  all  calculations  Poo  was  assumed  to 
be  10"4  atm.  As  with  flow  around  a  sphere,  the  radiation  effect  on  p  and 
V  was  found  to  be  practically  nonexistent  and  extremely  weak  with  respect 
to  the  distribution  of  the  velocity  component  v.  Also  analogous  to  the  flow 
around  a  sphere,  the  flow  difference  with  and  without  radiation  was  greatest 
at  the  body  surface.  The  temperature  distribution  parameter  was  the  one 
most  affected  by  shock  layer  radiation  (Fig.  2).  The  higher  the  flight  speed. 


Fig.  2.  Relative  temperature  distribution 
at  Voo  =  9.  85  km/sec,  R  =  1  m  with  (solid  lines) 
and  without  (dashed  lines)  allowance  for  radiation 
flow  along  the  lines  i?2  =  0.  747,  i?4  =  1.141. 
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the  more  noticeable  was  the  radiation  effect  on  T  distribution.  Radiation 
cooling  was  also  intensified  by  an  increase  in  R  (Fig.  3)  .The  plots  in 
Fig.  2  versus  those  in  Fig.  3  illustrate  the  effect  of  the  flow  velocity  around 
the  segment.  A  decrease  in  T  of  the  shock  layer  results  in  a  decrease  in 


Fig.  3.  Relative  temperature  distribution 
at  Voo  =12  km/sec  and  R  =  0.05-4. 

layer  thickness.  It  was  shown  that  the  radiation  energy  flux  to  the  stagnation 
point  increases  with  an  increase  in  R,  primarily  due  to  spectral  line 
radiation,  in  agreement  with  spectroscopic  theory.  The  universal  dependence 
of  the  ratio  of  radiation  flux  to  a  body  and  at  the  stagnation  point  is  also 
valid  for  a  spherical  segment.  An  earlier  conclusion  that  the  nose  of  a 
streamlined  body  is  one  area  where  radiation  heating  exceeds  convective 
heating  applies  also  to  a  spherical  segment. 

The  effect  of  radiation  near  the  rounding  point  on  gas  dynamic 
parameters  distribution  in  a  shock  layer  was  also  analyzed  on  the  basis  of 
the  energy  conservation  equation.  This  equation  indicates  that  as  Q  — ►  0,  the 
total  enthalpy  H  along  different  flow  lines  is  even  affected  by  small  radiation 
fluxes  (Fig.  4).  Fig.  4  illustrates  the  phenomenon  of  radiative  freezing 
which  occurs  when  a  particle  is  carried  out  of  the  region  affected  by  radiation 
to  a  colder  region  before  the  particle  energy  is  significantly  affected.  This 
phenomenon  determines  changes  in  distribution  of  gas  dynamic  parameters 
near  the  symmetry  axis  and,  in  particular,  causes  a  sharp  drop  of  H.  In 
regions  far  from  the  symmetry  axis,  particles  propagating  along  the  flow 
lines  do  not  lose  energy  because  of  the  high  velocity  and  low  Q  of  radiation 
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Fig,  4.  Distribution  of  total  enthalpy 
H  along  different  flow  lines.  £  is  the 
distance  traveled  by  a  particle  from  the 
shock  wave. 

flux.  The  radiation  near  the  rounding  point  therefore  does  not  significantly 
affect  distribution  of  aerodynamic  parameters  of  the  shock  layer.  Calcu¬ 
lations  showed  that  allowance  for  the  geometry  of  the  radiative  volume  is 
not  required  in  the  vicinity  of  the  rounding  point. 


Kokin,  G.  A.  and  Ye.  V.  Lysenko.  Pitot 
tube  pressure  in  supersonic  rarefied  gas 
flow.  MZhiG,  no.  1,  1972,  195-199. 

Readings  of  four  Pitot  tubes  in  a  low  density  wind  tunnel  at 
M  oo  =  2.5,  2.6,  3.2,  and  5.4  are  compared  to  Pirani  and  McLeod  pressure 
gauge  data.  The  experiments  were  designed  to  evaluate  a  correction  to 
Pitot  tube  readings  in  a  supersonic  rarefied  gas  flow.  The  d/D  parameter 
of  the  experimental  Pitot  tubes  (Fig.  1)  was  0.  1-0.2  and  0.  5-0.9  for  1,  2  and 
3  tube  types,  and  the  !/d  parameter  was  1  and  50-100  for  the  1  and  2,  3 
types.  The  mean  free  path  X  oo  of  molecules  in  a  free  stream  flow  was  0.  5, 
1.2-0.  9,  3.  3-2.4,  and  0.  5-0.  4  mm,  respectively,  at  the  cited  M  oo  values. 

The  Knudsen  number  (Kn  oo  =  X  oo  /D)  dependence  on  the  relative  Pitot  tube 
readings  at  the  cited  M  oo  is  also  plotted  in  Fig.  1,  where  Po*  is  the  measured 
Pitot  tube  pressure  and  po2  is  the  Pitot  pressure  behind  the  compression  shock. 
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Fig.  1.  Geometry  of  Pitot  tubes  and  experi¬ 
mental  plots  of  Po*/P02  versus  Kn  00  . 

Dotted  lines  -  theoretical  plots  for  r/R  =  0.5 
or  ~  0.863.  Arrows  indicate  theoretical 
values  corrected  for  thermal  conductivity  of 
the  tube. 


To  evaluate  the  correction  to  the  Pitot  tube  readings,  the 
correction  to  the  Pitot  tube  readings,  the  asymptotic  pressure  ratio  P0jr7Po2 
in  free  molecular  flow  around  the  tube  was  formulated  using  known  relations, 
including  that  of  Rayleigh.  The  derived  formula  served  to  calculate  and 
plot  P0f/p02  ratios  for  thermally  insulated  tubes  of  different  f/d  parameters 
and  for  Moo  =  1-9.  The  relation 
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was  derived  for  a  thermally  insulated,  long  cylindrical  tube  in  a  nearly 
free  molecular  streamline  flow,  where  S  x  =  M  oo  /V«,  ,  the  ratio  between 
the  stream  mass  velocity  and  the  most  probable  thermal  velocity  of  molecules; 
and  r  and  R  are  the  inner  and  outer  radii  of  the  tube.  This  relation  was 
derived  in  approximation  of  a  single  collision  between  the  incident  molecules 
and  those  emitted  by  the  tube.  The  p0*/p02  values  for  Moo  =2.6  and  3.2 
calculated  from  (1)  are  in  good  agreement  with  corresponding  experimental 
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values  (Fig.  1).  An  additional  correction  factor  C  for  thermal  conductivity 
was  introduced  into  (1)  to  obtain  theoretical  p0*/p02  values  at  M«  =5.4 
in  agreement  with  the  experimental  data.  The  experimental  and  theoretical 
data  were  used  to  plot  a  universal  Kn  dependence  of  the  parameter  £=  (p0*-p02) 
(PoF-Po2)~  »  independent  of  M  «  .  The  plot  can  then  be  used  to  determine 
M  oo  from  measured  Pitot  tube  pressure. 


Zakharchenko,  V.  F. ,  M.  V.  Tsvetkova, 
and  Ye.  E.  Borovskiy.  Supersonic  flow 
around  a  porous  cone.  PM,  no.  2,  1972, 

132-134. 

Supersonic  air  flow  around  a  circular  pointed  cone  with  a 
permeable  surface  was  investigated  in  connection  with  the  control  of 
incident  heat  flux  and  drag  on  streamlined  bodies.  Shadow  photographs 
were  obtained  of  the  unperturbed  flow  around  a  0.  197  rad  half  angle  cone 
at  0-0.24  values  of  the  injection  parameter  (  p  V)},  where  p\  and  Vj  are 
the  density  and  velocity  of  the  injected  gas.  The  Mach  number  M  oo  of  the 
incident  unperturbed  flow  was  2.9,  the  Reynolds  number  Re  oo  /m  was 
2.  6  x  10  ,  and  the  ratio  Tqoo  /Tqi  of  stagnation  temperatures  of  the  external 
and  injected  gas  was  1.  Photographs  were  used  to  determine  the  angle  ©s 
between  the  compression  shock  front  and  the  direction  of  the  unperturbed 
flow.  It  was  assumed  that  the  free  stream  boundary  is  also  conic  at  uniform 
gas  injection  through  the  conic  surface  and  that  the  pressure  along  the 
boundary  is  uniform.  The  boundary  can  then  be  assimilated  to  an  impervious 
surface  of  an  imaginary  cone  as  well  as  its  position,  i.  e.  the  angle  A© 
between  the  boundary  and  the  cone  surface.  Pressure  p  along  the  boundary 
can  therefore  be  calculated  as  in  a  problem  of  supersonic  flow  around  a 
cone.  Calculations  indicate  that  both  0g  and  A0  increase  continuously  with 
increases  in  injection  (  p  V)i,  the  rate  of  increase  depending  on  the  (  p  Y)j 
value.  The  pressure  coefficient  pig  on  the  boundary  was  calculated  using 
available  data  on  supersonic  flow  around  cones.  Plots  of  the  relative  pressure 
coefficient  pi  -  "p/p  along  the  cone  surface  at  different  (  p  V)i  values  (Fig.  1) 
show  a  qualitative  pattern  of  variations  of  the  pressure  coefficient  Pi  in  the 
presence  of  air  injection  in  comparison  with  the  flow  (p)  around  an  impervious 
cone,  i.  e.  in  the  absence  of  air  injection.  The  experimental  p^  -  p/p  value 
is  shown  (Fig,  2)  to  depend  significantly  on  the  coordinate  x/L  and  (  p  V)^ . 


Fig.  1.  Relative  pressure  coefficient 
versus  injection  parameter. 


Fig.  2.  Relative  pressure  coefficient 
versus  coordinate  x/L  at  fixed  points 
of  the  cone  generatrix. 


The  data  reveal  a  significant  dependence  of  incident  flow 
parameters  on  the  injection  pressure  and  make  it  possible  to  calculate  the 
drag  on  a  permeable  surface  cons.  Total  drag  calculations  should  include 
an  allowance  for  injection  drag  which  can  amount  to  10%  of  the  total  drag 
at  a  low  injection  pressure. 
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Bykov,  B.  P. ,  and  S.  I.  Kozlov.  Ion  and 
electron  kinetics  in  a  disturbed  ionosphere 
at  altitudes  of  100^-200  km.  Geomagpotizm 
i  aeronomiya,  no.  2,  1972,  340-342. 

Results  are  described  of  an  investigation  on  kinetics  of  ions 
and  electrons  at  altitudes  100  <  h  <  200  km  in  a  disturbed  atmosphere.  The; 
disturbing  action  is  assumed  to  significantly  increase  thb  degree  of  ionization 
a  of  air,  so  that  N^+,  and  O’*”  are  formed  approximately  in  direct 

relation  to  the  N2,  O2  and  O  contents  at  these  altitudes.  Investigation  ie 
limited  to  values  a  «  1,  so  it  is  considered  that  the  atmosphere  is  not  warmed 
up  due  to  a  sharp  rise  of  a}i.  e.  electrons,  neutrons  and  ions  of  gas  have  the 
high  level  of  temperature  distribution  existing  in  normal  ionosphere. 


F  orm. 
no. 

Reaction 

Velocity  const  *  cm^  per  sec 

1 

O’  +  x.  NO’  -I  N 

at«2.10-u 

2 

Of  -L  0.  os  0 

a.  -.=  4  ■  lO-ii 

3 

XS  •+■  0  -  0+  -1-  N.J 

33  <  l-10-i' 

-1 

i\v  +  0  NO*  +  N 

y,  -■■■  2 , 3 ■  1 0~w 

5 

N.'*  +  Oi-»  o,h  +  x. 

ttii~  MO-1’ 

6 

0.'  +  Ni-»  NO*  -[•  NO 

eta  <  1  •  10-'s 

7 

OS  +  e  ->  0  -!-  0 

r  (Oj+)  —  2,2-  It)-7  ( V’./OOO)-1 

8 

N.+  -\-  e  ->  N  -{•  N  ■ 

:<*  (NY)  3-  10”T  »0)'a  ■ 

0 

NO*  -|-  e  N  4  0 

a*  (NO ')  --1,4- 10-;  (7V:i00)-«,  200'  K  < 

<  Te  <  500 3  K 

a*(\0’)  **  0,'J- 10-7  /1000)-*.'»,  Te  >  500°  K 

Table  1. 


The  general  scheme  of  elementary  processes  used  in  this  work 
and  their  velocity  constants  are  given  in  Table  1.  On  the  basis  of  these 
reactions  kinetic  differential  equations  are  obtained  which  were  solved  by 
computer.  Calculations  are  done  for  a  wide  range  of  altitudes  and  values  of 
Neo,  at  day  and  night  conditions.  Results  of  calculations  for  daylight 
conditions  at  03  =  10"H  and  =  10"^  cm^/sec  are  plotted.  The  main 
characteristic  features  of  ion  and  electron  kinetics  in  a  disturbed  ionosphere 
are  seen  to  be  as  follows:  1)  The  role  of  N2+  ions  is  significant  only  at  the 
initial  moment  after  disturbance,  not  exceeding  in  most  cases  —  1  sec; 
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2)  At  lower  altitudes,  as  in  a  normal  atmosphere,  recombination  takes 
place  mainly  according  to  a  square  law.  However,  with  an  increase  of 
Neo  after  the  disappearance  of  N^  over  a  certain  time,  velocity  variation 
may  follow  a  linear  law.  The  duration  of  this  time  significantly  increases 
at  higher  altitudes;  3)  The  time  behavior  of  [NO+J  differ s  little  from 
predictions.  Even  at  upper  altitude  regions,  the  role  of  this  ion  is  small, 
except  when  Neo <  10^  cm*^;  4)  To  a  first  approximation,  at  all  altitudes 
the  main  molecular  ion  is  C>2  ;  5)  The  principal  difference  in  behavior  of 
ions  and  electrons  in  night  conditions  from  day  consists  in  the  fact  that 
the  dissociative  recombination  level  increases,  so  that  the  role  of  N2*  and 
NO+  ions  substantially  decreases;  and  6)  The  decay  time  of  electron  concen¬ 
tration  from  Neo  to  normal  values  of  an  undisturbed  ionosphere  has  only  a 
slight  dependence  on  the  value  of  Neo. 

These  conclusions  also  hold  up  under  varying  values  of  the 
unknown  velocity  constants  o.  3  and  All  calculations  were  done  assuming 

al  -  a(,  to  be  independent  of  temperature. 


Alimov,  V.  A.  ,  L.  M.  Yerukhimov,  and  T.  S. 

Pyrkova.  Theory  of  the  Fapypa/|  phenomenon 
in  the  ionosphere.  Geomagnetizm  i  aeronomiya, 
no.  5,  1971,  790-797. 

An  analysis  is  given  of  possible  factors  which  cause  stretching 
of  the  return  r-f  pulse  from  an  ionospheric  probe,  or  the  FSpreacj  effect. 

The  authors  use  the  diffraction  theory  of  pulsed  signal  spread  in  a  nonuniform 
ionosphere,  together  with  data  on  nonuniformity  distribution  vs.  altitude 
in  the  F-layer,  and  diffusivity  data  on  returned  signals.  Both  ground-based 
and  satellite  transmitters  are  considered,  and  the  variation  in  spread  effect 
at  different  latitudes  is  considered.  Using  an  averaged  model  of  the  F-layer 
structure,  the  authors  obtain  expressions  for  electron  density  variation 
AN/N  in  the  equatorial,  middle  and  high  latitudes.  Following  this  the  character 
istic  of  returned  pulsewidth  T  is  obtained  for  several  possible  AN/N  models, 
assuming  either  an  exponential  or  parabolic  variation  within  a  layer,  as 
tabulated  in  Table  1.  Graphical  comparisons  are  given  of  theoretical  and 
calculated  behavior  of  the  depth  of  the  diffusivity  zone  Az  =  c  T  as  a  function  of 
signal  frequency.  Fig.  1  gives  this  relation  for  the  extremes  of  equatorial 
and  polar  latitudes,  based  on  ionogram  data  from  the  Alouette  satellite.  These 
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=  const 

t~6 r* 

t  ~  <0 
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t  as  const 

T~<t>* 

(IN?  (2)  =  B/A  —  z 

—  . 

tw  const 

(AAT//V)1  (z)  =  Aieaz 

T~(0 

— 

(*NlN)*(i)  =  Ci+Ble'“  . 

T  ~  V  Cj  +  Bilto* 

— 

Table  1. 


Fig.  1.  Az(f)  for  satellite  signals, 
a  -  equatorial,  b  -  polar  latitudes 

results  show  a  distinct  variance  from  analogous  data  using  a  ground-based 
probe  system,  as  seen  in  Fig.  2,  indicating  the  discrete  nature  of  inhomo- 


Fig.  2.  A  z(f)  for  ground-based  signals, 
a  -  equatorial,  b  -  polar  latitudes 
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geneities  in  differing  portions  of  the  F-layer.  The  findings  appear  to 
establish  that  the  FSprea(3  mechanism  is  generally  caused  by  scattering 
of  the  signal  pulse  on  local  electron  density  variations  in  areas  on  the  order 
of  1--5  km  across.  The  variation  in  spread  behavior  is  then  governed  by 
the  noncorrelation  in  spectral  response  of  the  pulse,  as  it  is  reflected  from 
different  levels  of  the  F-layer. 
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Vashchenko,  V.  I.,  Yu.  N.  Matyushin,  A.  K. 

Parfenov,  Yu.  A.  Lebedev,  and  A.  Ya.  Apin. 

Heat  liberation  in  a  low-velocity  detonation  regime. 

FGiV,  no.  1,  1971,  121-126. 

The  heat  of  explosion  and  the  composition  of  the  products  of 
powdered  tetryl  explosive  conversion  were  determined  experimentally  for 
two  discrete  detonation  regimes.  Low-velocity  (D  =  2.  1  km/sec)  detonation 
was  emphasized  to  verify  an  earlier  proposed  mechanism.  Cylindrical  tetryl 
charges  (  f)  =  0.  95  g/crm ,  d  =  20  mm,  H  -  100  mm)  in  0-25  mm  thick  water 
envelopes  were  detonated  in  a  calorimeter  by  TNT/NaCl  (50/50)  or  RDX 
primers  initiating  low-velocity  detonations,  or  an  RDX  primer  initiating  a 
normal  velocity  (D  =  4.5  km/sec)  detonation.  The  measured  heats  of  explosive 
conversion  Q  outlined  in  Fig.  1  show  the  differences  between  the  two  detonation 
regimes.  The  gaseous  products  of  the  low-velocity  detonation  contained  a 
significant  amount  of  nitrogen  oxide.  The  solid  residue  of  unreacted  tetryl  was 


Fig.  1,  Heat  of  explosion  of  tetryl  vs. 
thickness  of  water  envelope.  1  -  normal 
detonation  velocity;  2  -  low-velocity 
detonation  '".ithout  allowance  for  unrelated 
tetryl;  3  -  low- velocity  detonation  with 
allowance  for  unreacted  tetryl. 

20%  of  the  initial  charge  in  a  12-20  mm  thick  envelope.  These  experimental 
data  confirmed  the  earlier  proposed  point  mechanism  for  low-velocity  explosive 
combustion,  in  which  explosive  particles  ignite  at  separate  points  on  the 
surface  and  burn  up  in  parallel  layers.  Detonation  velocity  D,  calculated  by 
the  method  of  successive  approximations,  was  3.7  km/sec  which  is  higher  than 
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the  experimental  D(2.  1  km/ sec).  This  discrepancy  is  apparently  due  to 
the  much  higher  relative  amount  of  the  unreacted  tetryl  residue  (~70%) 
in  the  reaction  zone,  which  determines  D.  The  measured  explosion  pulse 
of  a  low- velocity  detonation  was  lower  by  a  factor  of  1.  7  than  the  pulse  of 
a  normal  detonation  due  to  a  significantly  lower  detonation  wave  pressure. 

A  low-velocity  detonation  regime  is  extremely  unfavorable  in  industrial 
applications  of  solid  explosives,  because  of  the  presence  of  toxic  gases,  the 
increased  danger  of  ignition  of  CH4-O2  mixtures,  and  the  decreased  effective  - 
ness  of  the  explosion. 


Merzhanov,  A.  G.  ,  and  E.  A.  Shtessel'.  Thermal 
explosion  in  liquid  reactive  systems  in  the  presence 
of  heat  convection.  FGiV,  no.  1,  1971,  68-76. 


Thermal  explosion  experiments  were  carried  out  using  1:2  hexa- 
methylene  di-isocyanate  (HDI)  -  alcohol  mixtures  in  flat  and  cylindrical 
containers,  under  conditions  favoring  heat  transfer  by  natural  convection.  The 
main  purpose  was  to  correlate  the  thermal  explosion  limit  8^  with  the  Rayleigh 
criterion  of  natural  convection  where  g  is  gravity,  <P  is  the  coefficient  of  volume 
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-  «  *  R  Tt  ri 
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expansion,  T*  is  the  critical  spontaneous  ignition  temperature,  r  is  the 
characteristic  dimension,  v  is  the  kinematic  viscosity,  a  is  the  conductivity 
temperature,  and  E  is  the  activation  energy  of  chemical  reaction.  The  corre¬ 
lation  to  be  determined  was  formulated  as 


'•*  -  s  •  /  ;R;i)  (2) 

and  the  problem  was  reduced  to  finding  f(Ra).  The  kinetic  parameters  in  (1) 
and  in  the  expression  for  8^  and  v  of  the  reactive  mixtures  were  determined 
experimentally.  Thermophysical  parameters  for  calculations  of  Ra  and  8^ 
were  taken  from  the  literature.  T*  was  determined  by  the  abscissa  of  the 
inflection  point  of  the  experimental  plot  of  ^kTmax  versus  T0  of  the  heat  carrier. 
The  Frank-Kamenetskiy  parameter  8Q,  without  allowance  for  convection,  was 
calculated  using  a  known  formula.  The  tabulated  T*  experimental  data  exceeded 
by~30°C  the  T*  theoretical  data  owing  to  the  effect  of  natural  convection. 

The  experimental  plot  of  f(Ra)  and  its  interpolation  formula  were  used  to 
evaluate  the  coefficients  of  internal  aeff  and  external  a  heat  transfer,  and  the 
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temperature  distribution  within  the  reaction  zone  in  the  presence  of  free 
convection.  The  effect  of  convection  on  8#  dependence  versus  Ra  was 
strongest  at  high  a  and  disappeared  at  low  or,  corresponding  to  Biot  numbers 
ranging  from  2.4  to  6.  The  theoretical  analysis  of  thermal  explosion  in 
liquids  reacting  without  gas  evolution  in  the  presence  of  convection  revealed 
a  satisfactory  agreement  between  theoretical  and  experimental  f(Ra)  data. 


Kukhtevich,  V.  I,  ,  I.  V.  Goryachev,  and  L.  A. 

Trykov.  Zashchita  ot  pronikayushchey  radiatsii 
yadernogo  vzryva  (Protection  against  the  penetrating 
radiation  of  a  nuclear  explosion).  Moscow,  Atomizdat, 

1970,  190  p. 

An  engineering  design  method  is  presented  for  the  protection 
of  ground  structures  against  the  penetrating  radiation  of  a  nuclear  explosions. 

The  book  is  based  on  Soviet  and  Western  open  source  materials  published  in 
1958-1969;  Soviet  sources  make  up  about  50%  of  the  123  references. 

Physical  aspects  of  nuclear  explosion  are  treated  only  briefly, 
since  they  are  covered  well  enough  in  several  recent  Soviet  and  Western  mono¬ 
graphs.  Emphasis  is  placed  on  neutrons  and  gamma-rays  produced  at  the 
atmosphere  -  ground  boundary  and  their  penetration  through  a  typical  protective 
shield.  Numerous  published  data  on  characteristics  of  ionizing  radiation  are 
assembled  for  use  in  the  design  of  protection  against  nuclear  radiation. 

Formulas  and  numerical  data  are  given  for  calculation  of  the 
dose  rate  of  primary  gamma  radiation  propagating  from  a  point  source  and 
from  a  spherical  radioactive  cloud  through  an  infinite  atmospheric  medium. 

The  effect  of  a  shock  wave  created  by  an  explosion  on  the  propagation  of 
primary  gamma-rays  is  evaluated  in  the  cases  of  aerial  and  ground  explosions. 
Formulas  are  also  given  for  the  dose  rates  of  primary  and  secondary  gamma 
-radiation  with  allowance  for  the  shock  wave  effect,  as  well  as  the  time, 
explosive  force,  distance  from  the  explosion  center,  and  ascent  velocity  of 
the  radioactive  cloud.  Numerical  values  for  an  integral  dose  of  gamma 
radiation  at  a  given  distance  from  the  explosion  center  are  described  by  formulas 
for  low  and  high  explosive  forces,  and  plots  of  gamma- ray  accumulation  rate 
versus  time  are  shown.  Theoretical  data  are  included  on:  spatial  distribution 
of  slow  neutrons  released  by  a  nuclear  explosion  and  formed  in  air  by  moder¬ 
ation  and  diffusion  of  delayed  neutrons;  integral  neutron  albedo  from  the 
ground  surface  and  coefficients  of  the  angular  distribution  function  of  scattered 
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neutrons;  and  characteristics  of  neutrons  reflected  from  the  ground  by  a 
nuclear  explosion  in  the  atmosphere.  Calculated  data  are  also  presented 
on  integral  gamma -radiation  albedo  in  a  nuclear  explosion. 


Other  tables  give  theoretical  and  experimental  daift  on  transmission 
of  nonscattered  neutrons  and  total  neutron  dose  through  an  iron  shield  of 
varying  thicknesses  and  at  different  angles  of  incidence;  transmission  of  a 
total  dose  of  neutrons  through  water,  polyethylene,  and  iron  shields;  trans¬ 
mission  of  nonscattered  neutrons  from  a  nuclear  reactor  through  an  iron, 
concrete,  or  polyethylene  shield;  and  angular  distribution  of  the  dose  rate 
of  scattered  neutrons  beyond  an  iron  or  polyethylene  shield.  Experimental 
and  theoretical  data  are  plotted  and  tabulated  on  angular  distribution  (differ¬ 
ential  albedo)  and  integral  albedo  of  neutrons  back- scattered  from  water, 
ground,  and  iron  at  different  neutron  beam  angles  of  incidence. 

Theoretical  mass  coefficients  of  attenuation  of  nonscattered 
gamma- radiation  by  different  materials  are  tabulated  and  discussed.  Theo¬ 
retical  transmission  coefficients  of  gamma- radiation  scattered  by  iron  shields 
of  varying  thicknesses  and  at  different  angles  and  characteristic  angles  of 
scattering  for  specific  energies  and  materials  are  tabulated.  Formulas  used 
to  calculate  the  tabulated  data  are  given.  Empirical  formulas  are  cited  for 
calculating  integral  energy  albedo  of  gamma  quanta  back- scattered  from 
protective  shields.  The  approximate  formula  for  bremsstrahlung  energy  of 
beta-electrons  is  given  and  data  on  the  relative  contribution  of  spectral 
components  of  bremsstrahlung  were  tabulated.  Secondary  gamma- radiation 
yield  beyond  the  shield  is  analyzed,  and  the  numerical  values  of  the  constants 
for  calculation  of  the  yield  are  tabulated  for  water,  iron,  and  concrete. 

Calculated  yield  functions  for  different  shield  thicknesses  are  plotted. 

Two  methods  of  designing  the  protection  of  ground  structures 
are  outlined:  (1)  an  exact  method  involving  the  use  of  computers  and,  (2)  a 
simplified  method  using  approximate  formulas.  Both  methods  require  a 
knowledge  of  the  barrier  and  geometry  factors  of  radiation  attenuation.  The 
formulas  for  calculation  of  these  factors  are  given,  and  the  sequence  of 
mathematical  operations  associated  with  the  two  methods  is  described.  The 
precise  method  consists  of  six  consecutive  steps. 

The  effects  are  evaluated  of  the  multiple  scattering  of  radiation 
on  the  inner  walls  of  a  structure. 

Appendices  contain  8  tables  of  barrier  factors  of  transmission 
of  neutrons  and  gamma-radiation  through  various  materials,  and  two  tables 
for  yields  of  trapped  gamma-radiation  for  various  materials. 
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Korobeynikov,  V.  P.  ,  V.  A.  Levin,  and 
V.  V.  Markov.  Explosion  of  a  combustible 
gaseous  mixture.  IN:  Nauch.  trudy  In-t  mekh. 

Mosk.  un-ta  no.  11,  1971,  83-89.  (RZhMekh, 

1/72,  #1B21 1) 

A  theoretical  study  was  made  of  the  explosion  of  a  stationary 
inviscid  and  thermally  nonconductive  gas  mixture,  and  of  its  combustion 
with  allowance  for  ignition  delay  and  a  subsequent  reversible  reaction 
proceeding  simultaneously  in  both  directions.  The  initial  gas  flow  is  described 
by  strong  point-explosion  formulas.  A  strongly  supercompressed  detonation 
wave  decomposes  into  compression  shock  and  ignition  fronts.  Front  separation 
increases  in  time  depending  on  the  symmetry  pattern,  explosion  energy, 
energy  of  activation,  and  the  adiabatic  exponent.  In  the  second  phase,  energy 
liberated  by  a  chemical  reaction  is  taken  into  account  in  the  numerical 
calculation  of  the  flow.  It  is  shown  that  detonation  combustion  is  not  restored 
in  this  case. 


Rodionov,  V.  N.  ,  and  A.  P.  Sukhotin. 

Parameters  of  elastic  waves  generated  by  a 
spherical  explosion  in  metals  of  different 
strength.  FGiV,  no.  1,  1971,  142-146. 

An  experimental  study  was  made  of  the  effect  of  yield  point  o 
in  uniaxial  elongation  of  A-00  aluminum  (1),  D-16  annealed  (2),  and  D-16 
hardened  duralumin  (3)  on  the  parameters  of  an  elastic  compression  wave 
generated  by  a  0.  2  g  spherical  charge  of  TEN.  A  diagram  of  the  experiment 
is  shown  in  Fig.  1.  The  metal  samples  had  similar  mechanical  characteristics, 


P - - - izoi:r,ot.<kt - — . 


Fig.  1.  Experimental  setup:  1  -  explosive 
charge;  2  -  metal  block;  3  -  pickup  plate; 

4  -  guard  ring;  5  -  insulators;  6  -  cable. 
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with  the  exception  of  differences  in  strength.  The  measured  Owas  450, 
1600,  and  3(100  kg/cm  for  samples  (1),  (2)  and  (3).  Velocity  rise  time 
t  positive  phase  duration  x  and  maximum  mass  velocity  Um  were 
measured  on  oscilloscope  traces  of  free  metal  surface  displacement.  The 
experimental  plots  Tj  (r),  1 2  (r),  (Fig.  2),  and  lgUm(lgr)  indicate  that 
elastic  waves  are  generated  by  the  explosion,  and  the  observed  differences 
in  temporal  wave  parameters  are  the  effect  of  strength  only.  This  effect 


Fig.  2.  Velocity  rise  time  (a)  and  positive 
phase  duration  (b)  versus  distance  r  from 
the  explosion  center. 

was  also  evident  in  variations  between  the  measured  volume  of  the  explosion 
cavities  of  the  three  metal  samples  (1.3,  0.75,  and  0.42  cu  cm,  respectively). 
The  linearity  of  the  calculated  plots  of  the  relative  maximum  mass  velocity 
versus  distance  from  the  explosion  center  confirms  that  the  extent  of  plastic 
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deformation  is  determined  by  the  static  o  of  the  metal.  The  motion 
pattern  and  state  of  the  metal  in  the  elastic  region  was  reconstructed 
from  the  experimental  data  and  the  equation  of  plasticity.  The  effect 
of  o  on  the  motion  of  the  boundary  line  R*  at  which  the  difference  of 
principal  stress  is  equal  too  of  the  metal  is  shown  in  Fig.  3.  An  increase 


Fig.  3.  The  effect  of  0  of  the  2  and  3 
metals  on  boundary  line  R# 

in  a  affects  both  the  amplitude  and  temporal  parameters  of  R#  motion. 


Krupnikov,  K.  K.  ,  V.  F.  Kuropatenko,  A.  T. 

Sapozhnikov,  B.  N.  Simanov,  and  V.  A.  Simonenko. 

Calculation  of  explosions  in  a  medium  with  poly¬ 
morphic  phase  transition.  DAN  SSSR,  v.  202, 
no.  2,  1972,  300-301. 

An  explosion  in  a  quartz-like  medium  was  calculated  on  a 
computer  using  a  RAND(raschetadiabaticheskikh  nestatsionarnykh  dvizheniy) 
program  for  nonstationary  adiabatic  motion.  The  program  takes  into  account 
the  polymorphic  phase  transition  of  quartz  into  "stishovit".  A  two-liquid 
model  was  adopted  for  the  medium  with  a  gas  bubble  in  its  center;  high  pressure 
expansion  of  the  bubble  simulated  a  point  explosion.  In  the  region  of  continuous 
flow,  pressure  P  and  internal  energy  E  of  each  phase  were  calculated  from 
the  equations 

dE  +  PdV  =  0,  (1) 

p  =  jpx(p)  +rp(F-F,)  (2) 


where  P«  =  —p ($" -  1).  £,  4  8  -  f>/p0  =  0.  605». 
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nl  =  4,  112  =  2.7,  I\  =  r2  =  0.55,  E0  =  const.  P,  V,  and  E  were  corre¬ 
lated  by  a  P  =  P(T)  equilibrium  curve  and  the  equation  of  state  of  the  mixture, 

E-E.=  [TdP(T)  ldT-P(T)](V-V.) 

where  V*  and  E#  are,  respectively,  specific  volume  and  internal  energy 
of  one  of  the  phases  along  the  equilibrium  curve.  In  many  cases,  expansion 
shock  and  bifurcated  compression  waves  appear  at  the  phase  mixture  boundary. 
Allowance  was  made  in  the  RAND  program  for  these  waves.  The  calculated 
shock  wave  parameters  were  plotted  against  a  dimensionless  variable 
r  =  r/rg,  where  rg  =  [Q/(  p  0C^)J  x/3  (Fig.  1).  The  plots  show  that  P  and  U 


Fig.  1.  Dimensionless  pressure  P  and 
velocity  U  versus  time:  1  -  before  bifurcation 

of  the  first  shock  wave;  2,  3,  and  4  -  after 
bifurcation. 

behind  the  first  shock  wavefront  increase  significantly  before  dropping 
stepwise  at  the  front  of  the  expansion  shock  wave. 


Ocheretin,  V.  N.  Study  of  shock  wave  formation 
from  exploding  wires  in  water.  EOM,  no.  5, 

1  ‘>71,  39-41. 


A.  description  is  giv^n  of  an  experiment  on  the  formation  of 
shock  waves  in  water  from  an  exploding  wire.  The  experimental  sketch  is 
given  in  Fig.  1.  Tests  were  conducted  with  steel  and  copper  wires  of  0.  15  mm 
dia  and  90  mm  length,  each  with  an  initial  voltage  of  50  kv  and  a  discharge 
capacity  of  1/Af.  Photographs  were  taken  at  the  rate  of  600,000  frames 
per  second;  two  photographs  are  given  illustrating  steel  and  copper  wire 
explosions.  It  was  noted  that  the  channel  formed  from  the  explosion  of  a 
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Fig.  1.  Experimental  diagram.  TR  - 
high  voltage  step-up  transformer;  V  - 
rectifier;  IVPl,  IVP2  -  spark  gaps; 

B  -  test  vessel;  <t>  -  aluminum  foil; 

BP  -  exploding  wire. 

steel  wire  had  a  uniform  structure,  while  that  formed  by  copper  wire  was 
nonuniform.  Energy  levels  in  the  channels  were  measured  at  150  kj/cm3 
and  115  kj/cm3  for  steel  and  copper  wire  respectively.  From  the  test 
results,  the  author  recommends  use  of  steel  wires,  rather  than  copper  or 
brass,  for  its  superior  hydrodynamic  parameters. 
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Aristov,  V.  V. ,  K.  S.  Karplyuk,  and  V.  P. 

Pavlenko.  Three- wave  interaction  coupling 
effect  on  development  of  explosive  instability. 

UFZh,  no.  2,  1972,  307-314. 

A  theoretical  analysis  is  given  of  the  nonlinear  interaction  of 
positive  and  negative  energy  waves,  for  the  case  of  mutual  coupling  action 
among  the  waves.  The  study  derives  from  the  established  fact  that  interaction 
of  three  waves  with  a  negative  energy  wave  results  in  a  so-called  explosive 
instability,  where  the  amplitudes  of  all  waves  increase  without  bound  during 
some  limited  interval.  Specifically  for  the  four  wave  case  defined  by 
wave  vector  kj  and  frequency  the  problem  assumes  the  resonance 

conditions 

k  k*  =  kj,  k  +  k2  =■  kJt 

“ft  —  “ft.  «=  “ft,,  “k .+  “ft,  =  CD*,. 


Two  limiting  cases  are  examined:  the  first  in  which  the  phases  of  the 
interacting  waves  are  taken  to  be  correlated  over  the  interaction  interval, 
and  the  second  in  which  a  chaotic  phane  relationship  exists.  In  the  latter 
case  it  is  shown  that  for  the  four-wave  case  with  one  negative  energy  wave, 
all  interacting  with  a  weakly -turbulent  plasma,  explosive  instability  will 
always  occur. 


Vaynshteyn,  B.  I  V.  I.  Zenin,  and  I.  L. 

Babict  enko.  On  the  effect  of  dynamic  compression 
of  explosive  charges  on  detonability.  IN: 

Detonatsiya  vzryvchatykh  veshchestv  i  bezopasnost' 
vzryvnykh  rabot.  Moskva,  Izd-vo  Nedra,  1967,  77-82. 

A  study  is  reported  on  the  effect  of  explosive  charge  compaction 
on  its-  detonability.  Factors  treated  include  crushing  of  explosive  constituents 
and  uneven  compaction  within  the  charge  volume.  It  is  found  that  crushing 
of  components  of  type  PZhV-20  ammonite,  particularly  trotyl,  by  dynamic 
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compression,  causes  no  significant  loss  in  detonability,  nor  does  the 
uneven  density  resulting  from  the  compression.  It  is  assumed  that  the 
reduction  in  detonability  of  dynamically  compressed  explosives  results 
from  uneven  density  of  the  explosive  in  the  direction  of  the  diameter  of  the 
charge,  associated  with  significant  distortion  of  the  shape  of  the  detonation 
wave. 


Dremin,  A.  N.  ,  S.  A.  Koldunov,  and 
K.  K.  Shvedov.  Shock  wave-initiated  detonation 
of  low-density  explosive  charges.  FGiV, 
no.  1,  1971,  103*111. 

Detonation  experiments  using  low-density  (0.  78-lg/cm3) 
charges  of  TNT  and  tetryl  are  described.  Particle  size  was  0.  1-1  mm 
at  8-25  kbar  initial  shock  wave  pressure.  All  experimental  wave  velocity 
profiles  during  detonation  transition  exhibited  a  steady  decay,  in  contrast 
to  compacted  and  cast  explosive  charges.  Energy  liberation  is  consequently 
most  intensive  immediately  behind  the  shock  front.  The  depth  profiles  of 
the  shock  front  parameters  and  the  extent  of  the  transitior  region  depend 
on  the  initial  shock  pressure,  as  well  as  particle  size  a-  i  density. 


Lidorenko,  N.  S.  and  A.  V.  Sidyakin.  On 
the  possibility  of  thermal  exp  osion  in  fine 
metal  powders.  DAN,  v.  202,  no.  3, 

1972,  566-569. 

The  phenomenon  of  sudden  high  temperature  jump  during 
the  heating  of  fine  metal  powders  and  its  subsequent  almost  instantaneous 
caking  is  explained  using  the  characteristics  of  caking  kinetics.  The  inter¬ 
dependence  of  the  '  ate  of  caking  and  the  temperature  increase  can  lead  to 
conditions  analogous  to  thosenecessary  for  a  heat  explosion.  Frenkel's 
(ZhETF,  16,  29,  (1946))  model  is  used,  where  the  process  of  powder  caking 
is  represented  as  the  viscous  flow  of  powder  particles  under  the  effect  of 

surface  tension  forces.  An  advantage  of  this  model,  besides  simplicity,  is 
in  its  explanation  of  the  process  of  the  transformation  of  the  surface  energy  of 
powder  particles  into  heat.  An  equation  describing  caking  kinetics  was 


obtained  by  Frenkel  by  equating  the  surface  tension  forces  to  the  viscous 
friction  acting  on  each  particle.  This  model  is  applicable  only  ^or  the 
extremely  fine  powder,  with  particle  sizes  of  the  order  of  100  A  and  low 
caking  temperatures,  when  the  viscous  flow  effect  is  predominant.  Starting 
with  the  heat  conductivity  equation  and  equations  derived  for  the  Frenkel 
viscous  flow  model,  the  following  well-known  heat  explosion  equation  is 
obtained: 


dT  0  _r/T  .  .  m 

-5T  =  — '  '  +*4f 


(1) 


where  T  =  powder  temperature;  ©  =  (Z/4)  (3  o  /ePQ  a);  Z  =  powder 
coordination  number;  0  =  specific  surface  energy;  c  =  heat  capacity; 
a  =  particle  radius;  tq  =  (2rj0a/3  <y  );  0  =  characteristic  constant;  E  * 

characteristic  constant,  (E2£>  T);  X  =  heat  conductivity  of  powder.  The 
solutions  to  (1)  agree  with  the  experimental  data  depicted  in  Figure  1. 


Fig.  1.  Characteristic  caking  temperature 
T  as  a  function  of  time.  Solid  line  -  T  >  Tc, 
Dotted  line  -  T  Tc. 


Independent  of  the  model  used  to  describe  the  caking  phenomenon,  the  heat 
explosion  equation  gives  a  good  estimate  of  the  surface  energy  of  the 
substance.  This  relationship  has  been  experimentally  confirmed.  If  the 
powder  is  placed  in  a  container  having  a  characteristic  radius  r,  its  behavior 
is  determined  by  the  so-called  explosion  parameter  8.  (D.  A.  Frank-Kam- 

enetskiy.  Diffusion  and  heat  transfer  in  chemical  kinetics,  M.  ,  1967,  ch.  7) 
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where  T0  =  temperature  of  the  container  walls  and  r  =  characteristic  radius 
of  the  container.  For  the  values  of  8  ^  8C  (8  critical  depends  on  the  type 
of  the  container),  a  certain  stationary  temperature  distribution  prevails 
in  the  powder.  For  8  2.  8C,  the  stationary  temperature  distribution  is  no 
longer  posjible  and  a  thermal  explosion  takes  place.  The  rate  of  caking  is 
deduced  from  the  above  equation,  as  well  as  Tc  as  a  function  of  the  system's 
size.  The  authors  suggest  that  experimental  verification  of  this  equation  would 
be  satisfactory  for  other  models,  and  point  out  that  if  it  were  so  verified  the 
value  of  self-diffusion  activation  energy  would  be  definable. 
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Popov,  V.  N.  and  V.  I.  Chernyshov.  Equation 
of  state  of  a  helium-nitrogen  mixture.  IAN 
Energ,  no.  3,  1971,  162-165. 

A  procedure  for  statistical  mechanical  treatment  o.f  experimental 
compressibility  data  is  described,  to  determine  virial  coefficients  for  He, 

N2,  and  He-N2  mixtures  in  the  273  -  773°K  range  of  temperatures  and  at 
pressures  to  400  bar.  The  virial  coefficients  B,  C,  D,  E  are  given  in  poly¬ 
nomial  form  as  functions  of  T.  The  coefficients  of  the  polynomials  are 
tabulated.  The  virial  equations  of  state  for  He,  N2,  and  He-N2  mixtures 
were  developed  from  the  described  treatment  of  the  experimental  data.  The 
equation  of  state  for  a  He-N2  mixture  where  the  index  m  refers  to  the  mixture, 

'm  -  1  +  "rrf>  +  <t»Ts  +  ;fn!‘3  +  W-  ( 1 ) 

describes  the  experimental  results  within  +  0.20%.  Compressibility  factors 
and  density  p  values  calculated  from  (1)  were  in  satisfactory  agreement 
with  earlier  experimental  values  over  the  experimental  temperature  and 
pressure  ranges.  This  fact  indicates  that  the  caloric  data  calculated  from 
(1)  are  sufficiently  accurate  and  reliable  for  design  of  technological  processes 
and  solution  of  structural  design  problems. 
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3. 


Geosciences 


A.  Abstracts 


Bui  in,  N.  K.  New  data  on  the  crustal  structure 
of  the  Baltic  shield.  IN:  AN  SSSR.  Doklady, 
v.  198,  no.  3,  1971,  657-660. 

The  deep  structure  of  the  Earth1!?  crust  in  the  Baltic  shield 
has  been  studied  by  deep  seismic  sounding  (DSS)  as  well  as  by  seismological 
and  gravimetric  methods.  DSS  data  were  obtained  along  four  profiles  in 
the  northern  part  of  the  shield,  while  seismological  data  were  obtained 
along  the  350-km-long  Vyborg  -  Sortavala  -  Spaskaya  Guba  profile,  as 
well  as  at  four  permanent  seismograph  stations  (see  Figure  1).  The 
article  describes  the  results  of  seismological  research  and  gives  an 
estimate  of  the  crustal  thickness  for  the  entire  Baltic  shield. 


(Continued  with  figure  on  next  page) 
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Fig.  1.  Location  of  Profiles  and 
Seismographic  Stations  in  the  Eastern 
Part  of  the  Baltic  Shield. 

a  -  DSS  profiles  (1958-1965);  b  -  seismo 
graphical  profile  (1965-1967);  c  -  perma 
nent  or  mobile  seismographic  stations; 

(1  -  4  -  stations  of  the  USSR  Academy 
of  Sciences,  5-9  mobile  stations 
(number  in  denominator  is  the  average 
crustal  thickness,  based  cn  earthquake 
converted  waves,  in  km), 


The  section  along  the  Vyborg  -  Spaskaya  Guba  profile 
(shown  in  Fig.  2)  was  determined  from  both  converted  PS  and  com- 
pressional  waves  generated  by  distant  earthquakes  and  engineering 
explosions. 


Fig.  2.  Section  of  the  Southern  Part  of 
the  Baltic  Shield  (simplified),'  Constructed 
using  Converted  PS  Waves  of  Earthquakes 
and  Explosions  (1965-1967).  1  -  recording 

points  reduced  to  the  profile  line;  2  -  zones 
of  abrupt  changes  of  the  characteristics  of 
the  section  and  conversion  interface  relief, 
possibly  related  to  deep-seated  faults; 

3-7  -  averaged  position  of  the  seismic 
horizons  (conversion  interfaces): 

3  -  A  horizon  in  "granite"  layer,  4  -  Conrad 
discontinuity,  5  -  horizon  in  "basaltic"  layer, 
6  -  Mohorovicic  discontinuity;  7  -  other 
(unidentified)  conversion  interfaces  in  the 
crust  and  upper  mantle. 
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As  can  be  seen  in  figure  2,  the  crust  and  upper  mantle  are 
distinctly  divided  into  blocks  down^to  the  depths  of  60  km.  A  refracting 
interface  with  a  refractor  velocity  vr  =6.8  km/sec  at  a  depth  of  about 
10  km  and  a  refracting- reflecting  interface  with  vr  =  8.  5  km/sec  at  a 
depth  of  47-51  km  have  been  identified  by  compressional  waves.  The 
latter  has  been  interpreted  as  the  Mohorovicic  discontinuity.  The  main 
seismic  discontinuities  identified  through  earthquake  converted  waves 
are  the  M  horizon  in  the  39-50  range  and  the  A  horizon  in  the  8.  5  -  18  km 
range  (corresponds  to  the  refracting  interface  with  vr  =  6.  8  km/sec). 

The  average  parameters  of  the  crustal  layers  are:  1)  crustal 
thickness  is  44  -  5  km,  2)  depth  to  the  A  horizon  is  13  -  5  km,  while 
the  thickness  of  the  layer  bounded  by  the  A  and  M  horizons  is 
30.  5-1.5  km,  3)  the  ratio  of  the  thicknesses  of  the  upper  (above 
the  A  horizon)  and  lower  partn  of  the  crust  is  0.43.  Thus,  the  depth 
to  the  Mohorovicic  discontinuity  in  the  southern  part  of  the  shield 
exceeds  the  previously  accepted  depth  by  5-15  km  or  more  (based  on 
gravity  data). 

The  average  crustal  parameters  obtained  previously  by  DSS 
along  a  100-km-long  segment  of  the  Pechenga  -  Lovno  profile  are 
reviewed.  From  this  review,  it  wa|  determined  that:  1)  the  depth  to 
the  Mohorovicic  discontinuity  is  42  -  3  km;  2)  the  depth  to  the  horizon 
with  v  =6.8  km/sec  is  13  km  and  the  thickness  of  the  "layer"  with 
velocity  vp  =  6.  8  -  7.  5  km/sec  (lower  part  of  a  "granulite  -  basite" 
layer, after  Litvinenko,  1968)  is  29  km;  and  3)  the  ratio  of  the  thicknesses 
of  the  upper  (above  the  hurizon  with  vr  =  6.8  km/ sec)  and  lower  parts 
of  the  crust  is  0.  45. 

On  the  basis  of  seismological  studies  in  the  southern  part  of 
the  shield,  very  reliable  results  from  DSS,  and  data  obtained  in 
Finland,  the  crustal  thickness  of  the  Baltic  shield  has  been  estimated 
to  be  38  -  50  km,  with  an  average  value  of  42  -  44  km.  This  estimate 
of  the  crustal  thickness  of  the  Baltic  shield  is  in  a  eement  with  DSS 
data  on  the  Ukraine  and  Voronezh  crystalline  massifs  (Kosminskaya, 

I.  P.  ,  et  al.  Geofiz.  byull.  Mezhvedomstv.  geofiz.  komit.  no.  20, 

1969).  Based  on  the  present  data  derived  for  the  Baltic  shield,  it 
is  recommended  that  previously  derived  data  (Deminitskaya,  R.  M.  , 

Kora  i  mantiya  Zemli,  1967),  be  reinterpreted,  since  they  yield  a 
10-15  km  lesser  thickness  value  for  the  Baltic  shield  as  compared 
with  the  Ukrainian  shield  and  Voronezh  massif. 
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Andreyev,  A.  A.  Some  general  problems  of  the 
crustal  structure  of  Sakhalin,  based  on  deep 
seismic  sounding  and  gravity  data.  Geologiya 
i  geofizika,  no.  2,  1971,  61-68. 

A  combined  interpretation  of  deep-seismic  sounding  (DSS)  and 
gravity  data,  measured  on  the  island  of  Sakhalin,  was  conducted  to 
obtain  additional  information  on  the  deep  structure  of  the  crust  in  the 
region.  In  this  article,  sections  along  profiles  VI  -  VI  and  III  -  III 
(see  Fig.  1)  are  discussed,  with  emphasis  on  the  geological  inter¬ 
pretation  of  velocity  (or  density)  discontinuities  within  the  sedimentary 
cover. 
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Fig  1.  Location  of  Profiles. 

1  -  DSS  profiles  and  their  numbers. 

2  -  geological-geophysical  sections. 
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The  section  along  profile  VI  -  VI  was  derived  from  DSS  data 
from  profile  19  and  gravimetric  data,  as  shown  in  Figure  2.  The 
velocity  discontinuities  were  interpolated  so  that  their  morphology 
satisfied  the  observed  Ag  field. 
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Fig.  2.  Section  VI  -  VI. 


1  -  velocities  (densities)  3  -  faults: 

a)  from  geophysical  data  and 

b)  from  geological-geophysical  data. 

Seismic  discontinuities  with  vr  =  4.  8  and  vr  =  5.  0  -  5.  1  km/sec, 
observed  in  the  West  Sakhalin  mountains,  Susunay  depression,  and 
Korsakov  plateau  regions, are  interpreted  as  the  surface  of  the 
Upper  Cretaceous  and  Mesozoic  -  Paleozoic  formations,  respectively. 
However,  in  the  Tartar  Straits  region,  which  is  characterized  by  a 
density  distribution  different  from  adjacent  regions,  two  seismic  dis¬ 
continuities  with  v„  =  5.0  and  6.2  -6.4  km/sec  were  assumed  to 
correspond  to  the  above  interfaces.  Seismic  discontinuities  with 
vr  =  7.  4  and  v  =  8.  0  -  8.  1  km/sec  observed  over  the  entire  profile 
were  interpreted  as  the  top  of  basaltic  layer  and  the  bottom  of  the  crust, 
respectively.  The  crust  is  divided  into  blocks  in  the  upper  portion 
(above  the  basaltic  layer),  with  the  blocks  being  characterized  by 
different  velocity  and  density  distributions.  Gravity  effect  computed  for 
the  above  interfaces  showed  that  no  single  interface  constitutes  the  pre¬ 
dominant  part  of  the  calculated  Ag  anaomalies. 
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The  section  along  profile  III  -  III  was  constructed  using 
seismic  (DSS  profiles  18,  20,  30  and  a  RWCM  profile)  and  gravity 
data  (profile  III  -  III).  Two  versions  of  the  sections  which  differ  only 
in  the  structure  of  the  sedimentary  cover  are  shown  in  Figures  3a  and  3b. 
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Fig.  3.  Two  Versions  of  the  Section 
Along  Profile  III  -  III. 


The  crust  is  divided  into  three  blocks  corresponding  to 
1)  the  West  Sakhalin  anticlinorium,  2)  the  Poronay  depression  and 
western  part  of  the  East  Sakhalin  anticlinorium,  and  3)  the  eastern  part 
of  the  East  Sakhalin  anticlinorium.  Velocity  discontinuities  with 
vr  =8.0,  6.  7  -  7.  0,  and  6.  0  -  6.  5  km/sec  are  interpreted  as  the 
Mohorovicic  discontinuity,  the  top  of  basaltic  layer,  and  an  interface 
underlying  some  volume  of  Mesozoic-Palezoic  rock,  respectively.  It 
is  concluded  that  the  results  obtained  confirm  the  existence  of  faults 
extending  to  a  depth  of  8  -  10  km,  and  that  the  crustal  thickness  beneath 
the  axial  part  of  the  island  of  Sakhalin  amounts  to  33  -  35  km. 
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Antonenko,  A.  N.  ,  and  G.  K.  Dubrovin. 
Structural  characteristics  of  the  upper  mantle 
of  Mugodzhary.  IN:  AN  KazSSR.  Institut 
geologicheskikh  nauk.  Trudy,  no.  30,  1971, 
22-28. 


A  comprehensive  analysis  of  geophysical  data  obtained  in  the 
central  Ural  and  Mugodzhary  regions  is  made,  and  certain  assumptions 
are  derived  on  the  structure  of  the  upper  mantle  beneath  the  Mugodzhary 
folded  system.  The  results  of  seismological,  magnetotell  uric,  geothermal, 
and  gravimetric  studies  are  briefly  reviewed.  The  following  mani¬ 
festations  led  to  the  hypothesis  of  the  existence  of  a  structural  block  in 
the  upper  mantle  of  the  Mugodzhary  folded  system:  a)  a  conversion 
interface  at  a  depth  of  150-155  km;  b)  azimuthal  anomalies  of  -  4° 
(observed)  and  -  11°  (actual)  for  distant  earthquakes;  and  c)  the 

presence  of  a  high-conductivity  layer  with  its  surface  at  a  depth  of 
160-200  km. 

The  upper  boundary  of  the  block  is  assumed  to  coincide  with 
the  Mohorovicic  discontinuity  and  its  lower,  with  a  conversion  interface 
di  150-155  km;  its  thickness  is  100  km  with  its  width  corresponding  to 
the  Mugodzhary  system.  Its  existence  is  confirmed  by  gravimetric  and 
geothermal  data  (low  temperature  gradient),  as  well  as  seismological 
data  (foci  of  local  earthquakes  in  the  upper  mantle).  The  block  is 
characterized  by  high  seismic-wave  velocity,  excessive  density,  relative 
coolness,  high  specific  resistivity,  and  significant  viscosity.  It  is 
surrounded  by  a  low-velocity  layer  with  its  depth  of  occurrence  at 
80-100  km.  In  the  adjacent  Tyumen  -  Kustanay  downwarp,  the 
eastern  part  of  the  Russian  platform,  and  the  northern  Ciscaspian 
depression  regions,  the  proximity  of  this  layer  to  the  diurnal  surface 
causes  regional  negative  Ag  anomalies  and  high  geothermal  gradients. 


Suvorov,  V.  D.  ,  S.  V.  Krylov,  A.  L,.  Rudnitskiy, 
and  A.  L.  Krylova.  Deep  seismic  studies  of  the 
Earth's  crust  on  the  south  of  the  West  Siberian 
Platform.  Geologiya  i  geofizika,  no.  5,  1971, 

11 0"- 11 7. 

A  study  of  the  deep  crustal  structure  of  the  West  Siberian  platform 
has  been  conducted  using  point  deep-seismic  sounding  along  the  600  km 
long  Ishim  -  Barabinsk  profile.  The  techniques  of  generating  and  recording 
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of  seismic  waves,  as  well  as  the  dynamic  and  kinematic  characteristics 
of  the  recorded  wave  fields,  are  discussed.  As  a  result  of  the  study, 
data  have  been  obtained  on  crustal  layering,  the  block  structure  of  the 
consolidated  crust,  and  the  relief  of  the  Mohorovicic  discontinuity. 

The  velocity  section  along  the  Ishim  -  Barabinsk  profile  is  shown  in 
Figires  la  and  lb. 


U 


Fig.  1.  Velocity  Section  of  Ishim  -  Barabinsk  Profile. 

a  -  Velocity  section  of  the  upper  part  of  the  crust  and 
grapn  showing  refractor  velocity  at  the  surface  of  the 
basement,  b  -  Velocity  section  of  the  crust:  1  -  seismic 
interfaces,  <t>  -  surface  of  Prejurassic  basement,  I  -  II 
-  intercrustal  interfaces,  M  -  Mohorovicic  discontinuity; 

2  -  interfaces  from  less  reliable  data;  3  &  4  -  depths  from 
refracted  and  reflected  waves ;  5  -  boreholes;  6  -  probable 
zones  of  deep-seated  faults;  7  -  refractor  velocity,  km/sec; 
8  &  9  -  II  type  record  of  waves  from  the  Mohorovicic 
discontinuity. 
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The  basement  surface  is  a  sharp  discontinuity  with  vr  =  5.  3  -  6.  2  km/ sec 
at  a  depth  of  2  -  4.  5  km.  It  separates  the  Mesozoic- Cenozoic  sedimentary 
cover  (layer  velocity  2.  6  km/sec)  and  the  consolidated  crust.  The 
bottom  of  the  crust  is  at  a  depth  of  35  -  38  km,  with  a  refractor  velocity 
of  8.  0  km/sec  (determined  only  for  the  eastern  segment  of  the  orofile). 

The  mean  layer  velocity  of  the  consolidated  crust  is  6.  4  km/sec.  In  the 
upper  part  of  the  consolidated  crust,  two  interfaces  have  been  identified: 
intermittent  refracting  surface  I  at  a  depth  of  6  -  10  km  (v  =  6.  2  -  6.  4  km/ sec) 
and  reflection  surface  II  at  a  depth  of  22  -  26  1cm,  which  was  not  observed 
with  sufficient  confidence  along  the  entire  profile. 

Three  large  blocks  with  differing  deep  structures  and  velocity 
distribution  have  been  identified  and  designated  as  follows: 

1.  The  western  block  is  separated  by  a  zone  of  deep-seated 
faulting  on  the  east.  Its  crustal  thickness  amounts  to  34  -  38  km.  The 
basement  surface  is  characterized  by  a  significant  variation  of  refractor 
velocity,  i.  e.  ,  5.  4  -  6.  1  km/sec.  Refracting  surface  I  is  at  a  depth  of 

6  -  10  km,  with  vr  =  6.  3  -  6.  4  km/sec.  Reflecting  surface  II  at 
24  -  f  6  km  is  not  identified  with  sufficient  confidence. 

2.  The  central  block  which  is  characterized  by  a  more  homo¬ 
geneous  consolidated  crust  and  basement  surface,  with  a  stable  refractor 
velocity  of  6.  25  km/sec.  The  depth  to  the  Mohorovicic  discontinuity  is 
35  -  36  km,  and  it  is  a  very  sharply  defined  seismic  discontinuity,  as 

in  the  western  block. 

3.  The  eastern  block  is  separated  by  the  Omsk  zone  of  deep- 
seated  faults  on  the  west.  The  upper  part  of  the  consolidate  i  crust  is 
characterized  by  very  sharp  heterogeneity.  The  basement  surface 
plunges  from  2.  5  to  4.  0  km  in  the  Barabinsk  region.  Refracting  surface  I 
is  at  a  depth  of  7  -  8  km,  with  v  =  6.  2  -  6.  7  km/sec.  Reflecting 
surface  II  plunges  smoothly  toward  the  east  from  22  to  26  km.  The 
Mohoricic  discontinuity  in  this  block  is  not  as  sharply  defined  a 
discontuity  as  in  the  other  two  blocks. 
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Puchkov,  S.  V.  ,  and  D.  Garagozov.  Calculation  of 
the  effects  of  ground-water  level  on  ground -motion 
intensity  during  earthquakes  .  IN:  AN  TurkSSR. 

Izvestiya.  Seriya  fiziko-tekhnicheskikh,  khimicheskikh 
i  geologicheskikh  nauk,  no.  4,  1971,  55-63. 

The  effect  of  the  ground  water  level  in  unconsolidated  sandy- 
clayey  soil  on  the  intensity  of  ground  motion  during  earthquakes  is  con¬ 
sidered.  A  method  for  calculating  motion  velocity  in  compressional  and 
shear  waves  in  the  ground-water  level  is  developed  as  are  formulas  for 
evaluating  the  increase  in  earthquake  intensity  due  to  the  change  of  the 
ground-water  level.  It  is  determined  that  motion  velocity  in  seismic 
waves  varies  with  the  depth  of  the  ground-water  level  according  to  a 
hyperbolic  law  with  an  exponent  of  1/6.  The  increment  of  intensity  of 
an  earthquake  due  to  a  change  of  the  depth  of  ground  water  level  from 
Z2  to  ^  can  be  expressed  as  AB  =  1.  1  log  £  ,  where  z2  >  z^. 


An  experimental  hyperbolic  curve  representing  the  dependence 
of  the  increment  of  the  earthquake  intensity  on  the  ground  wp'.er  level 
has  been  obtained  for  the  Ashkhabad  district.  It  was  found  that  ground 
water  at  depths  exceeding  10-12  m  does  not  affect  ground-motion  inten¬ 
sity  at  the  surface.  However,  ground  water  at  shallower  depths  has  a 
considerable  effect  on  the  intensity  of  ground  motion. 


Bezrodnyy,  Ye.  M.  Seismic  shear-wave 
reflections  from  the  Conrad  discontinuity,  based  on 
seismological  data.  IN:  Uzbekskiy  geologicheskiy 
zhurnal,  1971,  1,  20-23. 

An  analysis  of  kinematic  and  dynamic  characteristics  of  shear 
waves  from  the  Tashkent  earthquake  aftershocks  was  made  for  the  purpose 
of  demonstrating  the  possibility  of  recording  near -normal  ( subcritical) 
reflections  of  SH  waves  from  crustal  interfaces.  On  the  basis  of  the 
analysis  of  observational  data  and  the  velocity  distributions  for  the 
region,  four  pronounced  phases  of  shear  waves  were  identified.  The&e 
were  interpreted  as:  S'  and  S"  -  waves  singly  and  doubly  reflected 
between  the  free  surface  and  the  Paleozoic  basement;  and  S'j^  -  waves 
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reflected  from  a  crustal  interface  at  a  depth  of  12-15  km.  Taking  into 
account  the  kinematic  characteristics  of  S'^  waves  and  previously 
obtained  data  on  the  Conrad  discontinuity  (H  «  10  km),  the  conclusion  was 
drawn  that  there  exist  possibilities  for  recording  near-normal  reflec¬ 
tions  from  the  Conrad  discontinuity.  The  "shadow"  zone  observed  for 
S'k  waves  was  explained  as  due  to  a  zone  of  crushed  rock  within  the 
focus  of  the  Tashkent  earthquake. 


Kurbanov,  M.  ,  and  V.  I.  Lykov.  The  relationship 
of  recent  tectonic  movements  and  the  seismicity  of 
southern  Turkmeniya  to  characteristics  of  crustal 
structure.  IN:  AN  TurkSSR.  Izvestiya.  Seriya 
fiziko-lekhnicheskikh,  khimicheskikh  i  geologicheskikh 
na.uk,  no.  1,  1972,  32-37. 

A  comprehensive  analysis  of  seismological,  geodetic,  and 
geological  data  was  conducted  to  determine  the  relationship  between 
recent  tectonic  movements  and  the  seismicity  of  southern  Turkmeniya 
to  its  tectonic  structure.  ReleveTng  showed  continuing  formation  of 
structure  of  the  first  kind  in  the  Ciskopetdag  downwarp  region,  where 
negative  structures  are  characterized  by  the  negative  sign  of  recent 
motion.  The  faults  identified  as  being  transverse  to  the  trend  of  the 
downwarp  coincide  with  zones  of  changes  of  sign  of  recent  motions. 

The  zones  of  highest  velocity  gradient  of  recent  motions  coincide  with 
the  zones  of  deep-seated  faults  occuring  in  recent  times.  The  hypo- 
centers  of  earthquakes  are  concentrated  along  these  faults,  from  the 
uppermost  part  of  the  crust  down  to  the  upper  mantle.  The  intersection 
of  two  systems  of  regional  faults  characterized  by  opposite  horizontal 
motions  leads  to  the  formation  of  zones  of  compression  and  tension.  A 
circular  seismicaLy  active  zone  of  the  transcrustal  type  divides  the 
South  Caspian  intrageosynclinal  megadepression  (with  a  suboceanic 
crust)  from  regions  with  continental  crust.  Recent  vertical  movements 
of  the  Earth's  surface  in  southern  Turkmeniya  are  the  result  of  tectonic 
development  of  a  vast  region  of  outhwestern  Asia. 
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Zunnunov,  F.  Kh.  Interpretation  of  reflected  waves 
in  deep  seismic  sounding.  IN:  AN  UzSSR.  Doklady, 
no.  10,  1970,  50-52. 

The  problem  of  selecting  the  effective  velocity  for  the  inter¬ 
pretation  of  time  distance  curves  of  reflected  waves  in  DSS  is  discussed. 

The  effective  velocity  in  the  layer  above  the  Mohorovicic  dis¬ 
continuity  was  evaluated  from  the  time  distance  curves  of  normal  and 
wide-angle  reflections  using  the  method  of  velocity  model  selection 
based  on  the  equality  of  distances  between  actual  and  apparent  shot 
points.  Time  distance  curves  were  computer  converted  into  lines 
toi(x)  using  velocities  of  6000,  6200,  6400,  and  6600  m/sec.  It 
2  was  found  that  the  best  result  was  obtained  with  a  velocity  of 

6200  m/sec.  The  time  distance  curve  converted  with  that  velocity 
was  then  used  for  plotting  the  seismic  section. 


Yepifanov,  A.  A.  The  crustal  thickness  in  eastern  and 
central  Ciscaucasia.  IN:  AN  GruzSSR.  Soobshcheniya, 
v.  65,  no.  2,  1972,  313-315. 

The  mapping  of  the  Mohorovicic  discontinuity  beneath  the 
eastern  and  central  Ciscausasian  region  has  been  accomplished  based 
on  multi-parameter  correlation  of  seismic  and  gravimetric  data. 
Multiparameter  correlative  equations  =  f  Agres  h  (where 
Hm  -  depth  of  the  Moho  discontinuity,  Ag  res  ^  '-  residual 
gravity  anomaly)  obtained  for  two  DSS  profiles  were  extrapolated 
throughout  the  entire  region.  The  crustal  thickness  over  a  large 
part  of  eastern  and  central  Ciscaucasia  does  not  vary  significantly, 
and  the  minimum  thickness  of  40  km  occurs  in  a  belt  confined  to  the 
Manych  downwarp  system.  It  increases  slightly  toward  the  north 
reaching  42-43  km,  while  it  increases  considerably  toward  the  south. 
The  maximum  thickness  of  62-63  km  correspcnds  to  the  central  zone 
of  the  Greater  Caucasus.  Another  depression  in  the  Mohorovicic  dis¬ 
continuity  is  confined  to  the  eastern  zone  of  the  Greater  Caucasus 
where  the  crustal  thickness  is  53  km. 


Bulin,  N.  K.  ,  Ye.  A.  Fronyayeva,  V.  I.  Bubnova, 

V.  N.  Tolmanov,  and  Ye.  I.  Erglis.  The  deep 
structure  of  the  southwestern  Altay  region  from 
seismological  data.  Sovetskaya  geologiya,  no.  4, 

1969,  97-109. 

An  analysis  was  made  of  PS,  P  and  S  waves  recorded  by  mobile, 
permanent,  and  temporary  seismographic  stations  from  (a)  distant,  nearby, 
and  local  earthquakes  and  (b)  engineering  explosions.  The  purpose  ox  the 
analysis  was  to  study  the  deep  structure  of  the  crust  and  upper  mantle  in 
the  southwestern  Altay  region.  In  the  May-October  1966  period,  observa¬ 
tions  were  conducted  with  mobile  seismographic  stations  at  51  points 
along  the  Charsk  -  Leninogorsk  (I  -  I)  and  Kokpetky  -  Zyryanovsk  - 
Stolboukha  (II  -  II)  profiles,  as  well  as  at  several  regional  seismographic 
stations  (see  Fig.  1). 
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Fig.  1.  Chart  Showing  Seismographic  Stations,  Shot  Points, 

and  Profiles. 

1  -  Points  of  observation  by  mobile  seismographic  stations; 

2  -  permanent  and  temporary  seismographic  stations;  and 

3  -  shot  points  of  engineering  explosions. 

Sections  (see  Fig.  2a  and  2b)  along  the  profiles  were  derived  from 
seismological  data,  taking  into  account  existing  data  on  gravity  and  magnetic 
anomalies  and  the  density  and  magnetization  of  rock  in  southwestern  Altay 
and  adjacent  regions. 
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1  -  Observation  points;  2  -  depth  of  conversion  interface,  based  on 
reliable  observation  data  on  PS  waves  from  distant  earthquakes; 

3  -  the  same,  based  on  less  reliable  data;  4  -  the  same,  based  on 
nearby  earthquakes;  5  -  the  same,  based  on  engineering  explosions; 
6  -  Predevonian  basement  surface;  7  -  reference  horizon  A; 


8  -  horizon  B  (Conrad  surface);  9  -  Mohorovicic  discontinuity 

10  -  strong  conversion  interfaces  in  the  upper  mantle; 

11  deep-seated  fault  zone;  12  -  lower  edges  of  magnetized 
bodies  (a)  less  than  30  km  from  the  profile  and  (b)  30  —  60  km 
from  the  profile;  13  -  earthquake  hypocenters  K  =  10;  14  -  K  = 

7  -  9;  15  -  faults,  based  on  geological  data;  16  -  reflecting 
horizons,  based  on  seismic  prospecting  data;  17  -  faults, 
based  on  seismic  prospecting  data.  Curves  above  sections 
represent;  maximum  Vp  distribution;  maximum  density  (  <r) 
distribution;  distribution  of  the  relative  amplitudes  of  the 
horizontal  components  of  compress ionai  (q  )  and  converted 
(qps)  and  waves  for  the  Mohorovicic  discontinuity  (W  is 
the  vertical  compression  wave  component);  distribution  of 
magnetic  (  ATa)  and  gravity  (Ag)  anomalies;  and  mean  rock 
magnetization  (J). 


The  following  conversion  interfaces  are  identified: 

1.  A  conversion  interface  at  a  depth  of  3.  2  -  9.  0  km  (corres¬ 
ponding  vy  =  5.7  -  6.1  km/sec),  interpreted  as  the  surface  of  the 
Predevonian  basement  or  as  the  surface  of  the  consolidated  crust. 

Its  depth  gradually  increases  in  a  southeasterly  direction. 

2.  Horizon  A  at  a  depth  of  13  -  19  km,  interpreted  as  an  interface 
within  the  "granitic-metamorphic"  layer. 

3.  Horizon  B  at  a  depth  of  23  -  29  km  (v  =  6.  7  -  6.  8  km/sec), 
interpreted  as  the  surface  of  the  "basaltic"  layer  or  Conrad 
discontinuity. 

4.  Horizon  M  at  a  depth  of  42  -  50  km  (v  =  8.  1  -  8.  3  km/sec), 
interpreted  to  be  the  Mohorovicic  discontinuity. 

5.  Horizons  M^,  M^.  and  M^  at  depths  of  72  -  93,  83  -  108,  and 
103  -  114  km,  respectively.  These  horizons  plunge  in  a  northeasterly 
direction. 

The  Mohorovicic  discontinuity  was  mapped,  with  two  zones  of 
minimum  depths  (43  -  45  km)  observed  which  are  confined  to  the 
Charsk  -  Kalba  -Narym  and  Rudno-Altay  structural  zones.  Maximum 
depths  (48  -  50  km)  were  observed  in  a  narrow  band  which  corresponds 
to  the  Irtysh  folded  zone. 
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The  maj  ority  of  the  vertical  seismic  discontinuities  identified 
were  confidently  interpreted  as  deep-seated  fault  zones.  All  of  these, 
except  for  the  Charsk  zone,  are  characterized  by  an  almost  vertical 
dip.  The  Charsk  and  Kalba-Narym  faults  extend  to  a  depth  of 
130  km,  while  others  are  considerably  less  deep. 

The  relief  of  the  Mohorovicic  discontinuity  and  crustal  thickness 
compare  favorably  with  results  derived  from  gravity  and  seismic 
data;  however,  a  significant  number  of  previous  determinations  of 
crustal  thickness  are  5-10  km  less  than  those  derived  in  the 
present  article. 


Kondorskaya,  N  .  V.  ,  and  A.  R.  Rakhimov. 

Results  on  the  computer  determination  of  the 
focal  position  of  the  1948  Ashkhabad  earthquake, 
using  observations  of  the  WWNSS.  IN:  AN  TurkSSR. 
Izvestiya.  Seriya  fiziko-tekhnicheskikh, 
khimicheskikh  i  geologicheskikh  nauk,  no.  3, 

1971,  53-62. 

Computer  calculations  of  the  coordinates  of  the  1948  Ashkhabad 
earthquake  focus  have  been  carried  out  using  the  observational  data 
of  the  world-wide  network  of  134  seismograph  stations,  within  a 
focal  distance  range  of  7  -  130  .  Epicentral  coordinates  were  calcu¬ 
lated  at  fixed  focal  depths  (h  =  0,  10,  15,  33)  from  various  sets 
of  data:  seismograms  recorded  at  134  stations  withA=7  -  130° 
and  seismograms  recorded  by  70-80,  48-50,  30-32,  20-22  and 
10-12  distant  stations  with  A  >20°.  The  results  obtained  are  shown 
in  the  form  of  tables  and  polar  diagrams.  The  accuracy  of  deter¬ 
mination  of  epicentral  coordinates  is  discussed  and  found  to  be 
sufficiently  high  with  a  0°.  1  semi  major  axis  of  "error"  ellipse. 

It  was  shown  that  epicentral  coordinates  do  not  vary  with  fixed 
focal  depth,  when  only  a  network  of  distant  stations  with  A  >20°  is 
used  for  the  calculation,  and  that  the  network  of  48-50  stations 
uniformly  distributed  around  the  epicenter  at  a  distance  A  >20° 
provided  limited  accuracy.  A  focal  depth  of  about  25  km  was 
determined  with  rather  low  accuracy,  since  the  nearest  observation 
station  was  at  A  =  7  .  The  coordinates  of  the  epicenter  were  com¬ 
pared  with  the  results  obtained  by  macroseismic  study  of  the  earth¬ 
quake.  The  epicenter  of  the  Ashkhabad  earthquake,  as  determined  from 
instrumental  data,  was  found  to  be  shifted  by  20  km  toward  the  north¬ 
east  with  respect  to  the  epicenter  determined  from  macroseismic 
data. 
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Kurbanov,  M. ,  K.  Dzhamanov,  and  T.  Tokareva. 
Tilt  observations  in  the  Ashkhabad  area.  IN: 

AN  TurkSSR.  Izvestiya.  Seriya  fiziko- 
tekhnicheskikh,  khimicheskikh  i  geologicheskikh 
nauk,  no.  5,  1971,  124-127. 

Data  from  tilt  measurements  in  the  seismically  active  area 
of  Ashkhabad  and  the  parameters  of  earth  tides  computed  from 
these  data  are  examined.  Observations  during  December  1968  - 
September  1970  were  conducted  at  the  Ashkhabad  station  by  two  sets 
of  tiltmeters  along  the  N-S  and  E-W  directions,  and  the  positioning 
and  installation  of  the  instruments  is  described.  A  tilt  of  the 
earth's  surface  of  4  arc  seconds  (northward)  along  the  N-S  leg 
and  1.25  arc  seconds  (eastward)  along  the  E-W  leg  were 
observed  over  the  above  period.  The  average  values  of  tidal 
parameters  were  computed  by  two  methods,  with  values  of 

7M2  =  °*  and  7M2  =  °*  obtained  by  fcbe  Matveyev  and 
Pertsev  methods,  respectively.  The  results  differ  significantly 
from  results  previously  obtained  for  the  same  area,  as  well  as 
from  the  results  for  adjacent  central  Asian  regions. 
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4.  Particle  Beams 


A.  Abstracts 


Abramovich,  V.  U.  ,  and  V.  I,  Shevchenko. 

Nonlinear  theory  of  dissipative  instability 
of  a  relativistic  beam  in  plasma.  UFZh, 
no.  2,  1972,  329-332. 

A  brief  treatment  is  given  of  conditions  for  dissipative 
instability,  occurring  upon  interaction  of  a  relativistic  particle  beam 
with  a  dense  plasma.  The  analysis  assumes  a  high  collision  frequency, 
and  is  performed  using  a  quasilinear  approximation  for  beam  relaxation. 

A  one-dimensional  case  is  first  considered  with  the  plasma  assumed  to 
be  in  a  strong  magnetic  field,  for  which  the  differential  equations  of 
relaxation  are  given.  A  three-dimensional  model  is  then  treated,  but 
without  the  presence  of  the  external  magnetic  fit' Id.  It  is  shown  that  the 
beam  energy  required  to  generate  oscillatory  instability  is  of  the  same 
order  of  magnitude  for  both  the  assumed  models,  in  any  case  a  criterion 
for  dissipative  instability  is  that  &  v  <  6/  k,  where  A.  v  is  thermal  velocity 
spread  in  the  beam,  5  is  increment  of  beam  instability,  and  k  »  wave  number 
of  the  unstable  harmonic.  It  follows  from  this  that  a  sufficient  spread  of 
electron  velocities  in  the  beam  would  prevent  the  cited  instability  from 
occurring. 


Ayvazov,  V.  Ya.  ,  and  B.  O.  Bertush. 

Forming  films  of  stoichiometric  composition 
during  electron- irradiated  evaporation  of 
AlpQ^  in  vacuum.  NM,  no.  2,  1972,  259-262. 

A  study  is  reported  on  factors  affecting  the  composition  and 
physical  properties  of  the  condensate  during  electron- irradiated  evaporation 
of  corundum  in  vacuum,  and  also  on  the  determination  of  optimum  conditions 
for  formation  of  AI2O3  films  of  stoichiometric  composition.  Dieleotric 
constant  e,  loss  tangent  8,  film  thickness  Y,  and  coefficient  of  optical 
absorption  e  were  experimentally  measured;  and  their  characteristic  curves 
as  a  function  of  the  parameter  Ng/Nm  are  drawn  (Fig.  1),  where  Ng  =  flow 
of  active  (liberated)  gas  and  Nm  =  molecular  flow  of  the  evaporating  substance. 
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Optimum  conditions  for  formation  of  AI2O3  film  close  to  stoichiometric 
composition  is  found  in  the  range  15<Ng/Nm<80,  where  e  =  8.  6  and 

9  g/cm2.  From  the  test  results,  the  following  conclusions  are  drawn: 
1)  Heating  of  AI2O3  by  a  focused  electron  beam  in  vacuum  gives  rise  to 
a  significant  dissociation  of  the  evaporating  substance;  2)  Condensate 
properties  are  governed  by  the  rate  of  precipitation,  temperature  of  the 
substrate  and  the  degree  cf  vacuum;  and  3)  The  composition  and  dielectric 
properties  of  the  film  are  determined  by  the  relation  Ng/Nm. 


-104- 


Kovalev,  V.  P. ,  V.  P.  Kharin,  V.  V.  Gordeyev 
and  V.  I.  Isayev.  Isotropic  neutron  source 
based  on  the  LUE-25  linear  electron  accelerator. 

Atomnaya  energiya,  v.  32,  no.  2,  1972,  173-175. 

Results  are  described  of  an  experiment  to  obtain  an  isotropic 
neutron  source  by  proper  selection  of  target  shape  and  material  in  a  LUE-25 
linear  electron  accelerator.  Angular  distributions  and  discharge  of  photo¬ 
neutrons  from  a  lead  target  bombarded  by  23  Mev  electrons  were  studied. 
Threshold  reactions  P^(n,  p)Si^,  Al^^(n,  p)Mg^7  and  Al^(n,  a  )Na^4 
having  effective  thresholds  of  2.7,  4.5  and  8.  1  Mev  respectively,  were  used 
for  detection  of  neutrons.  A  peculiarity  of  the  LUE  neutron  source  is  the 
presence  of  a  high  bremsstrahlung  background  noise.  Angular  distribution 
of  bremsstrahlung  from  the  target  was  also  studied  and  the  noise  level  was 
measured.  A  nearly  isotropic  angular  distribution  (within  limits  of  15%) 
was  observed  when  using  a  spherical  lead  target  of  60  mm  diameter  with  a 
20  mm  deep  well.  The  photoneutron  discharge  for  such  a  configuration  was 
registered  at  4  x  10*0  neutron//*  a  x  sec.  The  minimum  value  of  the  brems¬ 
strahlung  for  the  above  target  was  measured  at  an  angle  of  145°  to  the  electron 
beam  direction,  and  was  equal  to  1.  7  r//4.a  x  min.  The  peak  value  was  at 
zero  angle  and  at  90°,  at  which  the  noise  level  was  3.  3  r//*a  x  min.  Fig.  1 
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Fig,  1. 

shows  the  neutron  yield  as  functions  of  target  geometry. 
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Kolesnikov,  P.  M. ,  N.  S.  Kolesnikova  and  I.  B. 

Gabris.  Acceleration  of  plasma  in  a  coaxial 
device  with  optimum  inductive  energy  storage. 

I-FZh,  v.  21,  no.  6,  1971,  1117. 

A  theoretical  study  is  given  on  the  operation  of  a  coaxial 
plasma  accelerator  from  an  inductive  energy  storage  with  a  ferromagnetic 
core.  The  case  of  minimum  energy  loss  in  the  inductive  storage  is  studied, 
when  the  current  source  follows  the  relation  I  =  A  (Sh  at  +  b  Ch  at),  where 
I  =  current  of  the  charging  cycle;  a,b  =  constants  and  t  =  time.  Equations 
for  transient  processes  in  the  feed  circuit  and  in  plasma  acceleration  are 
solved  by  taking  into  consideration  the  mass-transfer  processes  expressed 
in  dimensionless  variables.  Results  of  the  quantitative  investigation  show 
that  inductive  energy  storage  with  an  optimum  law  of  current  variation  in  the 
charging  circuit  significantly  changes  the  character  of  transient  processes 
in  the  storage  accelerator  system  and  improves  characteristics  of  the 
accelerated  plasma.  Optimum  selection  of  parameters  can  lead  to  an  increase 
by  several  times  in  plasma  velocity,  compared  to  a  similar  accelerator 
operating  under  a  non-optimum  law  of  current  variation. 


Morozov,  A.  I.  Plasma  accelerators. 

Conference  in  Moscow.  VAN,  no.  7,  1971,  95-97. 

A  summary  is  given  of  some  200  papers  presented  at  the 
First  All-Union  Conference  on  Plasma  Accelerators,  held  February  8  -  12, 

1971  in  Moscow  under  sponsorship  of  the  Scientific  Council  on  Plasma  Physics  , 
AS  USSR,  and  the  Ministry  of  Higher  and  Secondary  Specialized  Education. 

All  plasma-dynamic  devices  discussed  are  included  under  the  general  heading 
of  the  plasma  accelerator(PA).  After  the  opening  speech  by  L.  A.  Artsimovich, 
A.  I.  Morozov  discussed  general  physical  problems  of  the  PA  which  he 
defined  as  a  system  for  ion  acceleration  without  disturbing  quasineutrality. 

The  two  basic  physical  problems  of  development  and  adjustment  of  the  PA 
arise  in  connection  with  the  near-wall  effects  and  plasma  instabilities.  An 
up-to-date  state-of-the-art  of  PA  research  and  development  was  summarized 
in  review  papers  presented  by  A.  A.  Porotnikov,  A.  A.  Kalmykov,  L.  V. 
Leskov,  M.  I.  Pergament,  A.  V.  Zharimov,  Yu.  V.  Yesipchuk,  and  A.  A. 
Plyutto.  The  review  makes  it  evident  that  a  further  development  of  the  PA, 
specifically  to  a  high-energy  (10-100  keV)  particle  accelerator  badly  needed 
for  fusion  and  astrophysical  R  and  D,  would  t?ke  several  years  of  collective 
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effort.  A  review  by  S.  D.  Grishin  and  the  papers  of  I.  N.  Golovin, 

G.  L.  Grozdovskiy,  I.  M.  Podgornyy,  L.  S.  Polak,  and  N.  N.  Rykalin 
dealt  with  the  present  and  potential  applications  of  PA  in  nuclear  fusion 
research,  aerodynamics,  experimental  astrophysics,  plasma  chemistry, 
and  metallurgy,  respectively.  Discussion  of  the  papers  led.  to  the  conclusion 
that  current  research  on  applications  of  high-energy  (^10  eV)  plasma  is 
not  yet  fully  developed.  The  author  concludes  that  development  of  the  PA 
would  considerably  advance  the  solution  of  the  nuclear  fusion  problem. 


Ternov,  I.  M.  ,  V.  G.  Bagrov,  Yu,  I.  Klimenko, 

O.  S.  Pavlova,  and  V.  R.  Khalilov.  Stimulated 
emission  from  relativistic  electrons  moving  in  a 
low  intensity  plane  wave.  VMU,  no.  6,  1971, 

691-695. 

An  approximation  of  weak  waves  analysis  is  presented  on  the 
problem  of  stimulated  emission  from  relativistic  electrons  moving  in  a 
monochromatic  plane  electromagnetic  wave.  Equations  are  obtained  for 
the  probability  and  the  intensity  of  radiation,  and  spin  effects  are  analyzed. 
Wavefunctions  describing  electron  motion  in  an  elliptically  polarized  plane 
monochromatic  electromagnetic  wave  are  given  as  exact  solutions  to  Dirac 
equations.  The  radiation  emitted  by  a  moving  electron  is  caused  by  the 
perturbation  effect  of  an  incident(quautized)  electromagnetic  wave.  The 
assumption  1  ^2»  which  can  be  experimentally  realized  with  powerful 

laser  beams  and  for  which  a  single  photon  emission  (or  absorption)  process 
is  the  most  probable,  is  used.  General  expressions  for  the  probability  of 
stimulated  single  photon  radiation  (dDSorption)  emitted  by  an  electron  in  a 
unit  time  interval  are  discussed.  The  given  equations  are  cited  as  being 
useful  for  numerical  computations  of  the  stimulated  emission  under  appropriate 
experimental  conditions.  The  total  probability  of  stimulated  emission  from 
an  unpolarized  electron  is  computed  as  well  as  the  total  emission  intensity. 

For  the  purposes  of  qualitative  analysis,  the  above  expressions  are  simplified 
using  additional  restrictions  on  radiative  transitions.  The  probability  of 
transitions  with  spin  reorientation  and  their  dependence  on  the  initial  (or  final) 
spin  orientation  are  considered.  The  spin  reversal  probability  determines  the 
degree  of  electron  beam  polarization.  Maximal  polarization  effect  is  obtained 
when  the  electron  spin  vector  is  parallel  to  the  direction  of  incident  electro¬ 
magnetic  wave  propagation.  Analysis  of  a  simplified  expression  for  the 
probability  of  transition  with  spin  reversal  shows  that  maximum  polarization 
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and  probability  are  obtained  with  a  circularly  polarized  incident  wave. 
Conditions  are  given  for  the  complete  polarization  of  electrons.  It  is 
estimated  that  considerable  polarization  is  possible  for  electrons  with 
energies  up  to  50  mev;  for  higher  electron  energies  the  lifetime  of  the 
initial  electron  state  and  the  transition  probability  decrease  rapidly.  In 
principle  it  is  therefore  possible  to  obtain  polarized  electrons  using  powerful 
laser  beams. 


Doroshenko,  A.  N. ,  and  I.  M.  Rayevskiy. 

Short-circuiting  of  pulse  solenoid  by  semi¬ 
conductor  diodes.  PTE,  no,  1,  1972,  118. 

Results  are  described  of  the  application  of  power  diodes  in 
series  for  shunting  of  a  7.  8  ka  pulsed  current  solenoid  at  a  storage  battery 
energy  of  100  kilojoules.  Four  power  diodes  were  connected  in  series,  as 
shown  in  Fig.  1,  each  operating  at  an  inverse  voltage  of  900  v  and  an  operating 
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Fig.  1.  Experimental  diagram. 


current  of  500  a.  The  sclenoid,  with  an  inductance  L  =  0.  7  mh,  was  fed 
from  a  condenser  battery  with  a  capacitance  C  =  15y«f  through  an  ignitron. 

The  diodes  were  connected  in  the  reverse  direction  with  respect  to  the  source, 
and  therefore  do  not  affect  circuit  operation  until  the  capacitor  voltage  drops 
to  zero.  As  soon  as  a  positive  potential  ?»ppears  on  the  anode  of  the  diode,  the 
Q  of  the  discharge  circuit  drops  sharply,  and  the  solenoid  current  shorts 
through  the  diodes.  This  circuit  thus  allows  the  short  circuited  solenoid  to 
pull  a  current  up  to  7.  8  ka  at  a  voltage  on  the  condenser  baxtery  of  3  kv. 
Comparative  oscillograms  of  current  and  voltage  pulses  with  and  without  the 
diodes  are  included,  showing  the  stretching  effect  of  the  diodes  on  the  current 
pulse. 
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Bakal,  S,  Z.  Effect  of  air  pressure  on 

the  erosion  action  of  anode  and  cathode  flares. 

EOM,  no.  5,  1971,  20-21. 

One  of  the  ways  to  decrease  the  erosion  of  electrodes  during  a 
drop  in  air  pressure  is  to  lower  the  erosion  action  of  a  flare.  In  this  work, 
the  relationship  between  the  erosion  effect  of  flare  and  the  air  pressure 
was  investigated.  Figure  1  gives  the  experimental  diagram.  The  cylindrical 


Fig.  1.  Experimental  sketch. 

1  -  ceramic  tube;  2  -  cylindrical 
electrodes;  3  -  vacuum  chamber; 

4  -  test  specimen;  5  -  high  voltage 
circuit  breaker.  A  -  discharge  gap; 
B  -  distance  between  flare- shaped 
electrode  and  specimen. 


electrodes  are  made  of  steel  with  a  1.  5  mm  diameter.  Erosion  was  deter¬ 
mined  on  a  cylindrical  steel  and  copper  specimen  of  6  mm  diameter  and 
10-12  mm  height.  Experiments  were  conducted  at  pressures  of  760,  100, 

10  and  10-  torr.  At  each  value  of  pressure,  six  series  of  five- pulse 
combinations  were  applied,  and  erosion  was  measured  after  the  fifth  pulse. 
After  every  discharge  series,  the  steel  electrodes  were  moved  by  0.  3-0.  5  mm, 
and  flush  ground  using  the  plane  face  of  a  ceramic  tube.  Results  of  the 


Fig.  2.  Relationship  of  erosion  action 
of  (1)  anode;  and  (2)  cathode  to  air  pressure. 
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experiments  are  given  in  Fig.  2.  The  erosion  effect  of  an  anode  flare 
at  atmospheric  pressure  is  seen  to  be  J.  5  -  2  times  greater  than  that  of 
a  cathode  flare.  When  the  air  pressure  was  lowered  from  760  to  10  "1  torr, 
the  erosion  effect  of  the  anode  and  cathode  flares  decreased  by  10  and  2 
times,  respectively. 


Zakatov,  L.  P. ,  A.  S.  Kingsen,  A.  G. 

Plakhov,  and  V.  D.  Ryutov.  Spectral 
investigation  of  electromagnetic  self- radiation 
in  a  plasma-beam  system.  ZhETF,  v.  61, 
no.  3,  1971,  1009-1015. 

The  object  of  this  work  was  to  measure  the  frequency  spectrum 
of  e-m  radiation  from  plasma  interacting  with  a  powerful  electron  beam. 


Fig.  1.  Experimental  device. 

I  -  electron  gun;  2  -  plasma  injector; 

3  -  main  magnetic  field  coil;  4  -  screw  typw 
coil;  5  -  vacuum  chamber;  6  -  horn 
antenna;  7  -  waveguide;  8  -  InSb  crystal; 

9  -  detector  head;  10  -  preamplifier; 

II  -  solenoid  for  linear  magnetic  field  in  crystal. 
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A  sketch  of  the  experimental  arrangement  and  the  ultrahigh  frequency 

spectrum  analyzer  is  shown  in  Figure  1.  The  magnetic  field  intensity 

was  1  koe  at  the  center  of  the  chamber,  and  5.25  koe  at  the  end.  The 

electron  beam  and  hydrogen  plasma  were  injected  in  the  trap  along  the 

axes  of  the  system  so  as  to  collide.  Plasma  concentration  was  10*2-  1013  crr>-? 

and  electron  temperature  Te  was  5  ev.  Beam  parameters  were:  energy 

®  ~  30kev,  current  10a,  pulse  duration  t  =  250^cs,  and  particle  concentration 

nb^  2  x  10*^  cm“ 3.  A  new  type  of  spectrum  analyzer  is  introduced  here  which 

is  based  on  the  existence  of  slow-damping  electromagnetic  waves  in  solid 

plasmas.  Oscillograms  obtained  using  the  analyzer  are  shown  in  Fig.  2, 

and  analytical  results  are  illustrated  in  Fig.  3.  It  is  seen  that  the  intensity 


Fig.  2.  Spectrum  analyzer  oscillograms 
at  various  plasma  concentration, 
a  -  n~7  x  10*2  cm“3;  b  -  n  ~ 4  x  10  ^2cm-  3; 
c - n~l x  10*2cm-3;  d  -  linear  magnetic 
field  H  . 

peak  lies  close  to  the  plasma  frequency  cupf,  and  the  radiation  power  of  a 
1  cm3  plasma  equals  Pwp,f2:2  x  10-2  watt/cm3  (Fig.  3a):  and  P2  a»  f  a 
5  x  10’3  watt/cm3  at  a  concentration  of  n=  4  x  1012  cm-3  (Fig.  3b).  The 
energv  density  of  the  Langmuir  noise  in  the  plasma,  estimated  on  the  basis 
of  radiation  with  a  frequency  of  2  c opf,  is  1.  4  ergs/cm3  -  2  x  1012  ev/cm3 
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Fig,  3.  Radiation  power  of  1  cm^ 
plasma  at  a  given  frequency  interval 
(fk»  fk+i)  *or  different  plasma  concentrations, 
a  -  n  ~ 7  x  10^2  cm"3;  b  -  n  ~  4  x  lO*2  cm" 3 ; 
c  -  n  ~  1  x  10^2  cm" 3 . 

at  n  s 4  x  1012  cm"3.  The  presence  of  radiation  with  a  frequency  wp(f 
indicates  that  low  frequency  oscillations  are  also  produced  in  plasma  along 
with  Langmuir  radiation. 


Lopatin,  V.  V.  ,  and  V.  Ya.  Ushakov. 

Equipment  for  investigating  optical  and 
electrical  characteristics  of  nanosecond 
breakdown  in  condensed  media.  PTE, 
no.  1,  1972,  144-146. 

A  description  is  given  of  a  device  consisting  of  a  nanosecond 
pulse  generator  of  amplitudes  to  400  kv  with  regulated  duration,  a  recording 
circuit  of  pulsed  predischarge  current  of  an  amplitude^  5  ma,  an  electron- 
optical  gate  and  an  optical  amplifier  (Fig.  1).  The  glycerin  condenser,  with 
a  capacitance  of  2 yuf,  is  charged  through  a  choke  coil  by  an  Arkrd'yev-Marks 
generator.  During  the  operation  of  the  arc  commutator  in  N2  at  a  pressure 
of  15  atm  a  voltage  wave  is  propagated  with  a  rise  time  of  3  nsec.  The 
experimental  discharge  chamber  is  located  in  a  section  of  the  Lj  line  and  is 
adjusted  in  accordance  with  the  impedance.  A  constant  discharge  up  to  96  nsec 
may  be  obtained;  pulse  length  is  controlled  by  the  shearing  discharger,  in 
which  discharge  takes  place  on  the  teflon  surface. 
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Fig.  1.  Experimental  device  for  observing 
nanosecond  pulses.  1  -  microsecond  pulse 
generator;  2  -  charge  choke  coil;  3  -  glycerin 
condenser;  4  -  arc  commutator  in  compressed 
N2;  5  -  shearing  discharger;  6,  8  -  capacitive 
voltage  divider;  7  -  oscillograph,  12-7;  9  -  T-junction; 

10  -  compensating  capacitor;  11  -  experimental 
discharge  chamber;  12  -  current  transformer; 

13  -  amplifier,  UZ-5A;  14  -  lens;  15  -  electron 
-optical  gate;  16  -  light  amplifier;  17  -  camera. 

Experiments  were  performed  with  two  types  of  dischargers  - 
a  nonregulated  dual  electrode  and  a  trigatron  type.  In  the  first  case,  the 
electrodes  were  placed  in  the  teflon  cylinder,  filled  with  an  insulating  liquid 
and  the  discharge  was  generated.  The  time  spread  of  discharger  operation 
does  not  exceed  50%  and  decreases  with  decreasing  pulse  duration.  An 
intensifier  arc  in  the  trielectrode  discharger  was  formed  in  the  1  mm  ring 
gap,  between  a  grounded  copper  cylinder  placed  in  the  teflon  cylinder  and  a 
di  sc  type  ignition  electrode,  fed  by  L4.  The  gap  between  main  electrodes 
is  regulated  during  varying  voltages  to  maintain  a  steady  over-voltage  coefficient. 
The  time  spread  of  discharger  operation  at  pulse  duration  T  j  =  40  nsec  was 
+  8  nsec,  and  at  t  j  =  12  nsec  =  +  1.  Onsec.  One-fourth  of  the  cylinder  volume 
was  filled  with  air  in  both  types  of  dischargers,  to  improve  the  shear  stability 
and  protect  the  teflon  surface  from  the  discharge.  The  liquid  dischargers, 
despite  their  simplicity  in  construction,  are  only  slightly  inferior  to  the 
more  complex  discharger  described  by  G.  S.  Korshunov  and  others  (ZhTF, 
no.  8,  1969,  1430). 

Compensation  for  capacitive  current  during  recording  of  the 
predischarge  current  was  done  using  a  bridge  circuit  with  Rogovskiy  trans¬ 
formers.  Transformer  parameters  were  selected  experimentally  using  a 
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model  comprising  a  G5-12  generator,  a  UZ-5A  amplifier  and  a  12-7 
oscillograph.  A. current  of  5  Ma  was  recorded  by  the  amplifier. 

Optical  recording  of  a  predischarge  phenomena  was  obtained 
according  to  the  method  given  by  Yu.  Ye.  Nesterikhin,  and  R.  I.  Soloukhin 
(Methods  of  high-speed  measurements  in  gasdynamics,  Nauka,  1967).  For 
optical  image  projection,  a  Yupiter-3  type  short  focusing  objective  lens  was 

used. 


Nagaybekov,  R.  B.  Ionization  processes  and 
ion  supercharge  in  the  cathode  spot  of  a  vacuum 
arc  discharge.  ZhTF,  no.  11,  1971,  2350-2352. 

Calculations  hav«  shown  that  for  the  existence  of  a  fast-moving 
arc  discharge  the  ion  current  density  should  be  very  high,  on  the  order  of 
105  -  106  a/cm2.  This  density  is  possible  only  in  the  case  of  intensive 
evaporation  of  the  cathode  and  a  high  effectiveness  of  ionization  by  electron 
collision.  A  mathematical  expression  is  therefore  derived  for  the  ionization 
probability  of  atoms  in  the  cathode  spot  and  curves  are  drawn  in  relation 
to  their  concentration  for  current  levels  of  10  -  10  a/cm  for  a  copper 
cathode  (Fig.  1).  It  is  seen  that  the  ionization  probability  is  close  to  unity  up 


Fig.  1.  Relationship  of  ionization 
probability  of  atoms  in  a  cathode  (copper) 
spot  to  concentration.  1  -  103  a/cm2; 

2  -  10^  a/cm2. 

to  an  atom  density  na  =7.5  x  101?  cm"3,  and  starts  dropping  at  na  >  7.  5  x  10*9 
cm- 3.  Results  of  the  calculation  of  atom  flow,  evaporation  rate,  cathode 
spot  temperature  and  atom  concentration  for  a  copper  cathode  at  various  E  and 
distances  from  the  cathode  are  given  in  Table  1.  Due  to  the  result  of 
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E,  j.’cu 

*«. 

CM-1  •  COK.-1 

a, 

t/cm*  •  ceK. 

T.  'K 

n(,  cif* 

10< 

7.5.10*2 

8 

2850 

7.75.10” 

10S 

7.5.102J 

80 

3400 

7.1.101* 

10«  . 

7,5  •  1024 

800 

4350 

6.25-101* 

Table  1. 

overcharging  of  ions  at  distances  of  5  x  10”^  -  10“^  cm  from  the  cathode, 
a  layer  of  ions  is  seen  to  exist,  creating  a  field  of  5  x  10^  v/cm. 


Bagirov,  M.  A.  ,  V.  P.  Malin,  and  Yu,  N, 

Gazaryan.  Breakdown  of  polymer  dielectrics 
under  the  action  of  electric  discharges 
resulting  from  thermo- oxidation  destruction 
of  their  macromolecules.  IAN  Az,  no.  3, 

1971,  9-14. 

Based  on  previous  results  of  Cooper  and  Prober  (Polymer 
Sci,  44,  no.  144,  I960),  the  possibilities  were  studied  of  erosion  due  to  the 
thermal  destruction  of  polymeric  bonds  at  contact  points  of  a  specific  arc 
discharge  channel  with  a  polymeric  surface.  Mathematical  expressions  are 
developed  and  calculations  are  made  for  the  temperature  field  of  an  arc 
channel  of  an  electric  discharge  in  a  gaseous  gap  between  dielectrics,  one  of 
which  is  a  polymer.  Results  are  plotted  in  graphs,  and  the  following  conclu¬ 
sions  are  made:  1)  At  contact  points  of  the  polymer  surface  with  the  arc, 
the  temperature  is  insufficient  for  thermal  decomposition  of  polymer  bonds 
in  an  oxygen-free  medium,  but  is  sufficient  for  thermo-oxidizing  decompo¬ 
sition;  2)  The  polymer  film  erosion  rate  increases  with  a  rise  in  the  tempera¬ 
ture  of  the  atmospheric  medium,  in  which  ionization  aging  takes  place. 


-115- 


Namitokov,  K.  K. ,  D.  P.  Solopikhin  and 
I.  Ya.  Surovtsev.  Radiographic  study  of 
structural  changes  in  electrodes  from  the 
effect  of  discrete  electrical  pulse  discharges. 

FiKhOM,  no.  6,  1971,  11-16. 

Structural  changes  (mainly  in  the  crystallization  process)  were 
studied  inside  and  near  a  unit  hole  on  electrodes  maae  of:  a)  annealed  and 
deformed  nickel,  and  b)  a  mechanically  compressed  mixture  of  copper  and 
nickel  powders.  The  electrodes  were  subjected  to  the  effect  of  a  single  high 
current  pulsed  discharge.  A  nickel  specimen  was  heated  in  vacuum  to  900°C 
for  two  hours  and  deformed  by  uniaxial  compression  so  that  its  thickness 
decreased  by  64%.  With  this  specimen  as  cathode  and  an  undeformed  nickel 
element  as  anode,  a  single  discharge  was  generated  in  air  from  a  battery 
with  C  =  800/^f,  voltage  U  =  240v  and  resistance  R  =  0. 1  ohm.  An  x-ray 
was  obtained  of  the  affected  area  after  the  discharge.  For  the  second 
specimen,  copper  powder,  pre-reduced  by  hydrogen  at  800°C  for  one  hour, 
was  mixed  with  nickel  and  pressed  at  5  T/cm2.  X-ray  records  obtained  for 
both  of  these  specimens  show  the  diverse  effect  of  the  temperature  gradient 
on  the  structural  change  of  electrode  materials.  A  special  x-ray  camera 
developed  for  the  experiment  is  described  and  illustrated. 


-116- 


*i**.»**m. 


Bek- Bulatov,  I.  Kh.  ,  and  R.  B.  Nagaybekov. 

On  the  state  of  the  cathode  spot  region  during 
arc  discharges  in  vacuum.  ZhTF,  no.  11, 

1971,  2383-2384. 

Characteristics  of  cathode  spots  as  functions  of  high  tempera¬ 
ture  and  voltage  gradients  are  examined.  On  the  basis  of  the  thermal 
balance  relaxation  time  of  atoms  and  the  time  between  two  successive 
collisions,  it  is  shown  that  the  basic  cathode  characteristics  (work  function 
of  electrons,  heat  of  evaporation,  temperature,  etc)  are  also  applicable 
to  a  cathode  spot,  taking  into  account  the  high  temperature  effect  and  the 
field  voltage  (~  107  v/cm).  The  temperature  T  and  its  rate  of  increase  in 
a  cathode  spot  were  determined  for  copper.  At  t  =  10-8  ~10"3  8ec  and  Q  =  107 
watt/ cm^,  re8ults  were:  T  =  4000  ~  4400°K  and  =  108~  109  deg/sec. 

The  findings  disagree  with  those  of  L.  B.  Leb  on  cathode  surface  properties 
(Fundamental  processes  of  electric  discharges  in  gas.  Gostekhizdat,  1950). 
Conditions  for  which  Leb's  conclusions  possibly  hold  true  would  require  an 
ion  current  density  of  10®  -  10^  a/cm^.  This  does  not  appear  to  be  a  reasonable 
figure,  since  to  maintain  such  an  ion  current  density  would  require  an  ion 
concentration  at  the  spot  close  to  that  of  a  solid  (~  1022  -  1023  cm-3). 


Konobeyevskiy,  S.  T.  ,  Zh.  S.  Takibayev, 

A.  G.  Abdullin  and  V.  F.  Grishchenko. 

Investigation  of  parameter  changes  in  a 
diamond  lattice  irradiated  at  various  tempera¬ 
tures  by  high-speed  electrons.  DAN  SSSR, 
v.  202,  no.  3,  1972,  560-561. 

Natural  diamond  powder  was  irradiated  in  a  U-10  linear  heavy 
current  accelerator  with  a  water-cooled  attachment  at  temperatures  of  80°C 
and  236  C.  The  specimen  was  examined  using  a  URS-60  x-ray  device  with 
a  RKU-114  camera  before  and  after  irradiation  by  2.3  Mev  electrons.  Results 
show  that  the  cyclic  change  of  a  diamond  lattice  is  a  function  of  electron 
radiation  dosage,  which  is  approximately  expressible  as  :  I  =  0.  11. 102* 

el/cm2,  as  shown  in  Fig.  1.  This  is  about  70  times  less  than  the  effect  of 
neutron  irradiation  (10/102^  n/cm2).  The  coefficient  is  valid  for  irradiation 
at  T  =  236°C  with  radiation  doses  of  4.  66xl020;  7.02xl020  and  1.05xl021  el/cm2, 
reflecting  a  nearly  linear  growth  rate  by  dosage.  A  similar  linear  relationship 
was  also  found  for  irradiation  at  80°C  (Fig.  1)  which  shows  that  a  change  in 
irradiation  temperature  does  not  change  the  stability  of  resulting  defects. 
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Fig.  1.  Change  of  parameters  of  diamond 
lattice  a  from  electron  radiation  dose  I. 

1  -  at  80°C,  2  -  at  236°C. 

It  was  also  noted  that  the  lattice  defects  induced  by  electron  bombardment 
were  permanent,  evidently  of  a  graphite  type  of  center,  which  remained 
even  after  anneal  at  1000°C. 


Lavrovskiy,  V.  A. ,  Ye.  G.  Shustin,  and 
I.  F.  Kharchenko.  Experimental  observation 
of  oscillations  in  electron  distribution 
functions  during  beam- plasma  interaction. 

ZhETF  P,  v.  15,  no.  2,  1972,  84-87. 

A  plasma-beam  discharge  was  obtained  in  hydrogen  in  a 
linear  magnetic  field  of  0.  2  ml  with  an  initial  beam  energy  of  1  kev  at 
currents  of  20-40  ma,  A  sawtooth  pulse  was  applied  as  a  beam  retarding 
voltage.  Pulse  amplitude  was  1.  5  kv  at  a  rise  rate  of  6.  109  v/sec.  The 
measuring  scheme  worked  in  a  single  stage  triggering  mode  during  stationary 
beam  injection,  with  a  measuring  circuitry  bandwidth  upper  limit  of  not 
les^  than  200  MHz.  High  frequency  oscillations  were  noted  in  a  range  of  1000 
MHz  while  LF  oscillations  were  observed  at  up  to  600  KHz.  At  sufficiently 
low  pressures  (p  =  5xl0"6  torr),  the  delay  curve  is  in  the  form  of  a  stepped 
pulse  and  its  differentiation  yields  a  bell-shaped  electron  distribution  function 
with  an  initial  dispersion  of  40  ev.  By  increasing  the  pressure  to  p  =  (1~3)10-3 
torr,  the  form  of  the  curve  changes  and  excitation  of  both  HF  and  LF  oscillations 
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occurs.  A  nonmonotonic  change  in  slope  of  the  delay  curve  with  a  change 
in  retarding  voltage  (Fig.  1)  ii,  a  significant  factor,  attesting  to  the  multi¬ 
peak  distribution  function  of  electrons.  Specific  electron  groups  have  an 
energy  dispersion  in  the  range  of  20-100  ev,  and  some  of  thes  s  groups  occur 
simultaneously.  With  a  further  increase  in  pressure,  the  multipeak  distri¬ 
bution  curve  characteristic  is  maintained,  but  a  general  spreading  of  the 
curve  in  the  direction  of  lower  velocities  is  noticed,  along  with  the  appearance 
of  particles  with  velocities  higher  than  Vo.  At  a  higher  recording  time 
constant  of  (T  =  10  msec),  oscillation-free  smoothing  of  the  distribution 
function  takes  place  with  the  formation  of  a  plateau  in  the  excitation  regimes 
of  the  most  intensive  H  F  oscillations. 
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Fig.  1.  Electron  velocity  distribution 
vs.  pressure.  Pressures  are:  a  - 
5  x  10"6  torr;  b  -  1,3  x  10"3  torr; 
c  -  2  x  10-3  torr‘  and  d  -  2.  7  x  10“3  torr. 


Vorob'yev,  A.  A.,  O.  B.  Yevdokimov,  and 
N.  P.  Tubalov.  Effect  of  charge  reversal 
in  a  dielectric  char ;(ed  by  an  electron  beam. 

FTT,  no.  12,  1971*  3691-3692. 

Tests  arc  described  which  demonstrate  that  the  volumetric 
charge  of  an  electron-bombarded  dielectric  may  reverse  after  sufficient 
elapsed  time.  Plexiglas  disks  of  4,  5  mm  thickness  were  covered  by  5-micron 
aluminum  foil  and  irradiated  by  electrons  with  0.  8-1.  2  Mev  energy  at  a 
current  density  of  0.5  /ua/cm^  for  15-20  sec.  To  determine  space  charge, 
the  effect  of  an  electric  field  on  the  passage  of  beta-particles  through  the 
dielectric  was  utilized.  The  charged  specimen  was  periodically  irradiated 
by  beta-particles  from  a  strontinm-iridium  source  and  the  number  of 
particles  before  and  after  radiation  was  measured.  The  relationship  of 
the  change  with  time  of  beta-particles  as  a  function  of  electric  field  is  shown 
in  Fig.  1.  At  energies  of  0.8  and  1.2  Mev,  the  quantity  AN/Nq  is  seen  to 
change  sign,  where  N0  and  N  are  the  total  particle  counts  before  and  after 
radiation.  This  corresponds  to  sign  change  of  effective  space  charge  in  the 


Fig.  1.  Relative  change  of  beta-particle 
count  with  time.  Charge  energies  shown 
on  curves. 

dielectric.  This  effect  is  attributed  to  the  electret  state  of  plexiglas  since 
the  effect  is  not  observed  in  polystyrene,  which,  in  effect,  does  not  display 
an  electret  state. 
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Musin,  A.  K. ,  and  M.  A.  Tyulina,  Deioni¬ 
zation  of  plaama  dispersing  in  vacuum.  ZhTF, 
no.  12,  1971,  2539-2544. 

Deionization  of  an  expanding  plasma  cloud  during  arcing  between 
refractory  electrodes  in  vacuum  was  investigated.  Experiments  were 
conducted  with  tungsten  electrodes  in  an  arc-quenching  chamber  as  shown 
in  Fig.  1.  Arc  quenching  was  simulated  by  a  thyratron.  The  theoretical 


Fig.  1.  Experimental  vacuum- 
arc  quenching  chamber  with  probes. 

Distance  from  the  chamber  axis  In 
mm:  1  -  23;  2  -  15;  3  -  40;  4  -  26. 

distribution  of  charged  particles  and  experimental  results  are  compared. 

A  change  in  charged  particle  concentration  with  time  was  determined  from 
the  rate  of  probe  current  drop.  At  an  arc  channel  of  radius  r~10"*  cm  and 
plasma  dispersion  rate  of  v<j  ~  10^  cm/sec,  the  theoretically  calculated 
concentration  drop  on  a  probe  placed  2  cm  from  the  arc  axis  begins  ~  2.  10“^  sec 
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after  arc  quenching,  whereas  the  characteristic  time  of  the  change  of 
concentration  at  a  given  point  is  ~10"7  sec.  These  values  show  gcod 
agreement  w?'h  experiments!  values. 


Ivanov,  N.  V. ,  G.  P.  Maksimov,  A.  I, 

Markin,  and  V.  P.  Smirnov,  Inductive 
titanium  hydride  film  injector.  ZhTF, 
no.  12,  1971,  2631-2633. 

An  inductive  titanium  hydride  film  injector  iB  described,  in 
which  a  plasma  bunch  is  formed  by  an  inductive  discharge  on  the  surface  of 
a  ceramic  disk  with  an  applied  layer  of  titanium.  The  disk  diameter  was 
11.  4  cm.  The  coil  had  a  diameter  of  13  cm  and  a  width  equal  to  10  cm. 
Voltages  to  50  kv  were  applied  through  a  discharger  with  a  capacitance  of 
0.  1  or  0.7/af.  The  value  of  the  electric  field  on  the  disk  periphery  reached 
650  v/cm.  Intensity  of  the  constant  magnetic  field  varied  from  250  to  1000  oe. 
Plasma  bunch  motion  and  parameters  were  studied  in  a  180  cm  magnetic 
circuit.  The  injector  contributed  to  the  attainment  of  a  lateral  cros Bisection 
homogeneous  plasma  bunch  of  10  cm  diameter  with  a  concentration  of 
1012  ~  1014  cm“3,  and  a  total  of  ionized  particles  up  to  10*8.  Results  of 
the  experimental  study  of  plasma  parameters  in  relation  to  discharge 
characteristics  are  shown  in  Figs.  1  and  2. 


Fig.  1.  Concentration  and  lateral 
electron  pressure  of  plasma  vs. 
discharge  voltage.  C  =  0.  l^af. 

Ho  =  500  oe. 
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Fig.  2.  Azimuthal  and  radial 
distribution  of  plasma  concentration. 
1  -  10;  2  -  40^tsec. 


Golovaniv skiy ,  K.  S. ,  and  A.  D.  Guseynova. 

Study  of  a  plasma  beam,  generated  by  a 
Penning  type  discharger  with  a  strong  radial 
electric  field.  ZhTF,  no.  9,  1971,  1876-1880. 

It  is  experimentally  shown  that  a  plasma,  generated  by  a 
Penning  type  source,  can  possess  an  anisotropic  electron  temperature.  The 
appearance  of  anisotropy  is  connected  with  electron  acceleration  in  a 
direction  perpendicular  to  a  magnetic  field  under  the  influence  of  a  strong 
electric  field  of  excited  discharge.  The  Penning  type  plasma  source  is  shown 
in  Fig.  1;  Fig.  2  outlines  the  overall  experimental  setup.  The  source 
discharge  chamber  is  a  copper  cylinder  of  140  mm  length  and  44  mm  diameter. 
A  stainless  steel  anode  of  38  mm  length  and  37  mm  diameter  is  placed  in 
the  middle  of  the  chamber.  The  anticathode  is  a  stainless  steel  disc  of 
42  mm  diameter  and  3  mm  thickness,  with  a  10  mm  diameter  aperture  at 
its  center.  The  cathode  is  a  tungsten  wire  of  0.  5  mm  diameter,  which  is 
heated  by  a  current  of  about  10  a.  The  cathode  and  anticathode  are  held  at 
zero  potential  while  a  positive  potential  of  up  to  500  v  is  applied  to  the  anode. 
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Discharge  current  was  varied  over  the  200 ~ 400  ma  range.  The  discharge 
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Fig.  1.  Plasma  source. 

1  -  cathode;  2  -  water  coil; 

3  -  anode;  4  -  screen;  5  -  anti¬ 
cathode;  6  -  plasma  beam;  7  -  glass. 


Fig.  2.  Experimental  setup. 

1  -  valve;  2  -  plasma  source; 

3  -  probe;  4  -  frequency- modulated 
signal  (  \  =  12  mm);  5  -  recording  of 
resonance  curve;  6  -  pumping. 

chamber  was  cooled  by  water.  The  magnetic  field  in  the  vicinity  of  the 
source  was  varied  between  100  to  600  oe;  argon  pressure  in  the  source  was 
kept  at  2. 10“  torr.  A  relationship  of  the  energy  of  fast  and  slow  electrons 
as  a  function  of  source  anode  voltage  Ua  was  obtained  (Fig.  3).  Results 
of  measurements  of  lateral  and  linear  electron  temperatures  Te  and  Te>, 
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Fig.  3.  Relationship  of  ion  energy 
and  space  potential  versus  source 
anode  voltage.  1  -  fast  ions;  2  -  slow 
ions;  3  -  space  potential. 


Fig.  4.  Relationship  of  lateral  and 
linear  electron  temperature  to  anode 
voltage. 

in  relation  to  Ua  are  plotted  in  Fig.  4,  The  maximum  anisotropy  factor 
T  =  Tej_  /Te  ||  was  5.  It  is  concluded:  1)  A  Penning  source  with  a  strong 
radial  electric  field  can  generate  a  plasma  beam  with  a  concentration  to 
10!0  cm” ^  without  additional  ion  extraction  from  the  discharge;  2)  Beam 
ions  are  divided  into  two  groups  according  to  energy:  fast  ions,  with  a  non- 
compensated  space  charge,  which  accelerate  with  a  cathode  potential  drop 
and  an  energy  level  equal  to  50  -  70%  of  the  anode  voltage;  and  slow  ions, 
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with  energies  of  20  -  40  ev,  which  fall  on  the  beam  due  to  the  ambipolar 
propagation  of  gas-discharge  plasma  along  the  magnetic  field. 


Gavrilenko,  T.  B. ,  and  G.  O.  Karapetyan. 

Resistance  of  activated  glasses  to  a  current 
of  accelerated  electrons.  ZhPS,  v.  15,  no.  6, 

1971,  1110-1112. 

Tests  on  cathode  resistance  to  electron  bombardment  are 
described.  The  article  deals  with  several  cathodolumine scent  glasses, 
and  treats  the  composition  and  excitation  state  factors  determining  the 
resistance  characteristic.  Experiments  were  conducted  with  borax,  silicate 
and  borosilicate  glass,  activated  by  trivalent  terbium  and  cerium.  Chara¬ 
cteristic  curves  for  aluminized  specimens  of  various  glass  under  similar 
radiation  conditions  (V  =  20  kv,  j  =  5.  10“7  a/cm2)  are  given  in  Fig.  1. 


Fig.  1.  Relative  decrease  of  brightness 
of  the  cathode-luminescence  glow  of  glass 
specimen  during  5  hours  of  continuous  exci¬ 
tation  (V  =  20  kv,  j  =  5. 10"7  a/cmZ), 

Dashed  curve  is  for  18-2  glass  at  V  =  20  kv, 
j  =  ia$/cm2.  Curve  numbers  correspond 
to  glass  numbers  in  Table  1. 

The  molecular  composition  of  the  glasses  is  listed  in  Table  1,  in  which 
activator  concentration  is  given  in  wt  %  oxide  above  100.  It  was  noted 
that  the  oxide  content  of  the  glass  is  the  determining  factor  in  resistance 
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to  accelerated  electron  current.  An  optimum  concentration  of  boric 
oxide  exists  in  a  series  of  borosilicate  glasses,  which  if  increased  sharply 


Table  1.  Composition  of  glass  (Mol  %) 
and  disintegration  parameter  c,  given 
by*’  I  =  I0/(l  +  CN),  where  I  -  initial 
intensity,  IQ  •  intensity  of  aged  screen, 

N  -  total  no.  of  electrons  falling  on  the 
specimen  during  radiation. 

lowers  the  stability  of  the  specimen  (Fig.  2).  Measurements  revealed  that 
glass  containing  about  6-7  mol  %  of  boric  oxide  was  the  most  stable.  Oxides 


Fig,  2.  Relative  decrease  of  brightness 
of  borosilicate  glass  with  0.  5%  Ce203 
and  varying  boric  oxide  content.  Curve 
numbers  are  the  same  as  in  Table  1. 

of  alkali-earth  elements  in  glass,  and  their  atomic  weights,  do  not  affect 
resistance  to  electron  irradiation,  but  additions  in  alkaline  oxide  content  do 
lower  the  resistance.  The  effect  of  temperature  on  the  darkening  of  the 
operating  layer  of  glass  was  also  investigated.  Dark  spots  appearing  in  glass 
specimens  after  electron  bombardment  vanished  when  the  specimens  were 
heated  to  450°C. 
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Kalyatskiy,  I.  I. ,  G.  S.  Korshunov,  and 


G,  A.  Kiselev.  Variation  in  spark 


resistance  in  water  from  the  action  of  a 
high  pulsed  voltage.  EOM.  no.  6.  1971.  32-37. 


Relationships  are  experimentally  determined  for  high  voltage 
resistance  in  a  spark  gap  in  water,  in  terms  of  the  specific  resistance  of 
the  water,  the  electric  field  configuration  and  the  spark  gap  width  6.  Fig.  1 


Fig.  1.  Experimental  diagram. 
Li,L2»L3  -  cable;  Gen  -  generator 
(Arkad'yev-Marks),  Pj  -  peaking  dis¬ 
charger  in  transformer  oil;  Pd  -  air 
gap.  §  -  spark  gap  in  water 


is  a  diagram  of  the  principal  experimental  circuit  and  Fig.  2  illustrates 
the  most  complex  element  of  the  experimental  arrangement,  the  coaxial 
discharge  chamber.  The  experimental  results  are  plotted  (3  curves  and 
1  oscillogram  are  given),  from  which  the  authors  make  the  following 


Fig.  2.  Experimental  discharge  chamber. 

1  -  cable  KPV- 1-300;  2  -  peaking  discharger; 
3  -  outputs  to  voltage  divider;  4  -  spark  gap 
in  water;  5  -  transformer  oil. 
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conclusions: 


1)  The  findings  have  application  in  the  calculation  and 
construction  of  high  voltage  electrode  systems  operating 
in  water,  and  in  the  devebpment  of  high  voltage  pulse¬ 
shaping  circuits,  loaded  by  spark  gaps  in  water, 

2)  A  hemisphere-hemisphere  electrode  system  possesses 
the  highest  resistance  during  the  application  of  w  high 
voltage  pulse  in  using  a  uniform  electric  field.  This  has 
significance  in  the  construction  of  a  peaking  nanosecond 
type  discharger  operating  in  water. 

3)  A  sharp  decrease  of  high  voltage  resistance  was  noted 
in  a  non-uniform  electric  field  (pointed-base  electrode 
system).  This  should  be  considered  when  designing  an 
electrode  system  for  operation  in  water. 

4)  To  reduce  load  pulse  deformation  and  to  obtain  a  high 
resistance  spark  gap  in  water,  high  voltage  rectangular 
pulses  with  nanosecond  fronts  are  recommended  for  various 
electrophysical  methods  of  treatment  and  destruction  of 
solids. 


Artyukh,  V.  G. ,  L.  G.  Lisenko,  and  S. 
Smirnov.  Circuit  for  rapid  switching  of 
heavy  currents  in  an  inductive  storage, 
no.  1,  1972,  119-120. 


A. 

PTE, 


A  circuit  is  described  which  allows  commutation  of  currents 
at  a  given  moment  of  time  (Fig,  1).  Capacitor  batteries  1  and  2  are  pre¬ 
charged  by  a  charging  device,  not  shown  in  the  figure.  On  triggering  the 
discharger  Pj  in  the  discharge  circuit  (battery  1,  discharger  Pj,  inductive 
storage  L  and  wire  P^  in  series),  a  sinusoidal  commutation  current  flows 
through  the  circuit.  Loads  Ljj  and  Rjj  are  shunted  by  P^,  so  that  practically 
no  current  flows  through.  On  triggering  discharger  P2  in  its  discharge 
circuit  (battery  2,  discharger  P2,  wires  Pfc  and  Fg  in  series),  a  regulating 
current  of  a  value  higher  than  the  commutation  current  starts  to  flow.  The 
discharger  P2  is  triggered  at  the  moment  of  the  first  maximum  of  switched 
current.  Due  to  the  regulating  current,  the  wires  Pfc  and  Pg  explode  and 
the  switched  current  flows  through  loads  Ljj,  Hjf  and  Rsh* 
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Fig.  1.  Experimental  scheme. 

Pj,  P2  -  regulating  discharger; 

L  -  inductive  storage,  130/xh;  Pg  - 
additional  exploding  wire;  P^.  -  commu¬ 
tation  wire;  Rjj  -  resistive  load  0. 1  ohm; 

Ljj  -  inductive  load,  2/<.h;  Ksh  -  shunt; 

Battery  1  -  3 3 60yU  f,  5  kv;  Battery  2  - 
140/if,  5  kv. 

To  increase  the  pressure  in  the  gap  formed  after  the  explosion 
of  wires,  the  wires  are  placed  in  dis  “'argo  chambers  of  the  type  shown 
in  Fig.  2.  A  description  of  the  dischargers  and  trigger  devices  used  in 
this  experiment  are  given  in  a  paper  by  I.  I.  Aksenov  et  al  (PTE,  no.  3, 
1970,  166). 


Fig.  2,  Discharge  chamber  for  wires. 

1  -  insulator;  2  -  steel  tube;  3  -  outlets; 

4  -  wire;  5  -  discharge  channel;  6  -  channel 
widening. 

The  maximum  value  of  switched  current  obtained  in  the 
experiment  was  approximately  25  ka  for  wires  of  14  mm  diameter.  The 
time  spread  of  the  initial  moment  of  switching  was  mainly  determined  by  a 
momentary  instability  of  the  synchronizing  device,  and  did  not  exceed  a 
few  tenths  of  a  microsecond.  Tests  showed  that  the  device  permits  simple, 
reliable  switching  of  currents  up  to  10^-- 10^  a,  using  laboratory  type 
inductive  storage. 
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Krivitskiy,  Ye.  V. ,  V.  V.  Shamko  and 
V.  L.  Apostoli.  Analysis  of  energetic 
parameters  of  an  underwater  spark  discharge 
path.  EOM,  no.  5,  1971,  48-50. 

Approximate  mathematical  expressions  are  derived  for  para¬ 
meters  of  an  underwater  arc  discharge  channel.  The  expressions  define 
time  behavior  of  current  power  and  resistance  in  terms  of  the  remaining 
circuit  parameters.  From  these  relations  an  expression  is  derived  for 
maximum  power  Nmax  in  an  underwater  path  of  a  given  length.  This 
expression  has  been  experimentally  verified  and  results  tabulated  (Table  1). 


Path 

No. 

CIO*, 

* 

U0  ■  10-«, 

V 

L '  10*. 

h 

i .  in-* 

m 

Wm.*  •  »_,.  W 

UiXpt. 

l 

9.0 

4,7 

1,16 

10 

11,8 

10,7 

•t 

9,0 

4,7 

1.16 

4 

6,3 

6.35 

3 

9.0 

4,7 

12,0 

9 

2,65 

2.78 

4 

2,83 

4,55 

7.5 

12 

2,9  . 

2,83 

5 

2,83 

4.65 

7.5 

6 

2,0 

1.65 

0* 

1,0 

7.0 

4,3 

5 

2,6 

7* 

1,0 

6.2 

3,8 

8 

3,85 

3.5 

8 

1,0 

4.0 

3,2 

6 

1,98 

1,78 

9 

1,0 

4,0 

3,2 

4 

1,51 

1,47 

10* 

1,0 

3,8 

3,2 

6 

1.9 

1.6 

11 

1,0 

4.4 

18,9 

4 

0,41 

0,3 

12 

0,1 

4,7 

4,01 

6 

0,76 

0,53 

13 

0.1 

4.8 

4.04 

4 

0,64 

0.55 

U 

0,1 

4.7 

4,04 

2 

0,4 

0,33 

Table  1.  Maximum  power  transmitted 
in  paths  of  length  f . 

*  Data  of  V.  V.  Arsent'yev  (PMTF,  no.  5, 
1965)  and  G.  S.  Bondarenko  (IVUZ  Elek- 
tromekhanika,  no.  5,  1968). 


Table  1  shows  that  the  difference  between  the  experimental  and  calculated 
results  does  not  exceed  25%,  which  is  acceptable  for  practical  applications. 

It  is  thus  shown  that,  given  only  the  initial  parameters  of  the  discharge  circuit 
and  the  distance  /  between  electrodes,  it  is  possible  to  calculate  magnitude 
and  time  to  maximum  power,  eliminating  the  need  for  experimental  work. 

No  other  channel  parameters  are  given  in  the  cited  study. 
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5.  Miscellaneous  Interest 


A.  Abstracts 


Mandzhikov,  V.  F. ,  A.  P.  Darmanyan, 

V.  A.  Barachevskiy,  and  Yu.  N.  Gerulaytis. 

Photochromism  in  organic  compounds  from 
the  action  of  laser  radiation.  OiS,  v.  32, 
no.  2,  1972,  412-413. 

Photochromic  intramolecular  transformation  of  a  spiropyran 
molecule  was  observed  in  the  experiments  with  spiropyran  solutions  in 
benzene,  ethanol,  and  polymethylmethacrylate,  which  were  excited  by  a 
ruby  laser.  Output  energy  and  pulse  length  were  0.  4j  and  3  x  10“8  sec 
in  the  experiments  with  a  Q-switched  laser,  4j  and  3  x  10"4  sec  in  those 
with  a  free-running  laser,  respectively.  A  colored  filament  appeared 
in  a  colorless  spiropyran  solution,  ranging  from  8-10  mm  to  40  mm 
long,  after  a  laser  beam  was  focussed  on  the  center  of  the  cell.  Coloration 
vanished  spontaneously  in  ~10  min  and  in  a  few  minutes  in  liquid  and  in 
polymer  solutions,  respectively.  The  colored  filament  was  not  generated  by 
ruby  laser  radiation  in  spiropyran  solutions  which  did  not  absorb  in  the 
300-400/<  range,  nor  by  neodymium  glass  laser  in  solutions  of  photochromic 
spiropyrans.  The  experimental  plot  of  D  versus  I2  (Fig.  1)  and  the  absence 
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Fig,  1.  Optical  density  D  of  the  8'-allyl  - 
6'-nitro  -  1,3,3  -  trimethylspiro  -  (2'H-l- 
benzopyran  -  2,  2'-indoline)  solution  in  benzene 
versus  square  of  intensity  of  exciting  ruby 
radiation. 

of  electronic  absorption  bands  at  14,400  cm"1  frequency  in  the  solutions 
studied,  indicate  a  two^photon  mechanism  of  the  spontaneously  reversible 
coloration  of  these  solutions. 


Dunin,  S.  Z.  Theory  of  spin  waves  in 
metals.  FMiM,  no.  1,  1972,  28-35. 


In  a  contribution  to  the  spin  wave  theory  of  a  Fermi  fluid  with 
allowance  for  anisotropy  of  the  Fermi  surface,  an  electron  fluid  with  a 
random  Fermi  surface  was  analyzed.  It  was  assumed  that  the  spin  component 
of  the  interelectronic  correlation  function  is  expressed  by 

V  (p,  p')  =  +  «V  (p)  aMVfv  (P')  ( 1 ) 

where  ey,„  =  ct  v/*  is  a  Landau  correlation  coefficient  of  conduction  electrons. 
In  approximation  of  (1),  the  dispersion  equation 


1  +  f +  “1 1*  <w-  k>  +  *v  («.  k)  Ovj^m  (®.  k)]  =0  (2 ) 

Ll  +  P,  WT  J 

was  derived  for  excitations  of  an  electron  fluid  in  a  steady  and  uniform 
magnetic  field.  In  (2),  fiQ  is  another  Landau  correlation  coefficient  and  t  is 
the  pulse  relaxation  duration.  The  equation  defines  the  spectrum  of  spin  wave 
frequencies  week.  An  analysis  of  Eq.  (2)  showed  that  the  w  spectrum  of  long 
waves  including terms  depends  significantly  on  the  Landau  anisotropic 
coefficients.  A  simple  analysis  of  (2)  revealed  that  w  of  short  waves  propagating 
normally  to  the  magnetic  field  is  limited  by  (nS2  extr  +.  Qq)»  where  flextr 
are  the  extremums  of 


°  -  2*t*  ar 

on  the  Fermi  surface  e  (p),  and 


(3) 


S20  =  2YB0/h  (4) 

In  (3)  and  (4),  B0  is  the  magnetic  induction,  S  (  e  ,  pz)  is  the  area  of  intersection 
of  the  e  (p)  =  surface  with  the  pB0  =  const  plane,  Y  is  the  effective  momen¬ 
tum  of  a  conduction  electron  and  h  is  the  Planck  constant  divided  by  2v  .  At 
n  =  0,  the  frequency  spectrum  of  spin  waves  monotonically  approaches  k  + 
value.  At  n  4  0,  w(  Q  ,k),  where  0  is  the  angle  of  the  magnetic  field  vector 
with  Z  crystallographic  axis,  oscillates  around  its  mean  (nS2extr  (  Q  )  +  fl0) 
v/hich  can  be  determined  experimentally.  In  principle,  the  Landau  coefficients 
of  the  Fermi  liquid  can  be  determined  from  the  dependence  of  the  spectra 
obtained  in  the  short  wavelength  range  on  the  reciprocal  orientation  of  magnetic 
field  and  crystallographic  axes. 
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Kiryushchenko,  A.  I.  and  M.  A.  Lebedev. 

Development  of  a  low-voltage  arc  discharge 
in  cesium  vapors.  I.  Development  of  a 
discharge  with  formation  of  a  spherical  plasma 
bunch.  ZhTF,  no.  6,  1971,  1170-1173. 

In  connection  with  the  problem  of  direct  conversion  of  thermal 
into  electrical  energy,  electrical  discharges  in  cesium  vapors  were  investigated. 
The  conditions  were  examined  of  luminous  spherical  plasma  bunch  formation 
in  an  inter-electrode  gap.  Discharge  volt-ampere  characteristics  were 
simultaneously  measured  by  an  electrically-heated  movable  tungsten  wire, 
and  emission  was  photographed.  The  volt-ampere  characteristic  and  photo¬ 
graphs  of  a  discharge  in  a  11.2  mm  gap  between  the  18  mm  diameter  electrodes 
at  1.  4  x  10“2  torr  cesium  vapor  pressure  and  cathode  temperature  Tc  =  790°K 
are  shown  (Fig.  1).  Spherical  plasma  bunch  formation  and  spread  over  the 
discharge  gap  coincide  with  a  1-1.  5  v  jump  in  voltage  and  point  6  of  the  v-a 
characteristic,  respectively.  Spherical  bunch  duration  is  limited  by  the  cesium 


and  emission  picture  of  a  low- voltage 
discharge  in  different  development  phases. 

vapor  pressure  (5  x  10-3  to  5  x  10“2  torr),  the  gap  length  (8sdsl8  mm), 
and  the  electrode  diameter  (8SDS20  mm).  The  position  of  the  bunch  within 
the  gap  depends  on  Tc  and  the  external  magnetic  field.  The  v-a  characteristic 
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of  a  Langmuir  probe  shows  that  electron  energy  distribution  within  the 
spherical  bunch  does  not  obey  Maxwell's  law.  A  sharp  emission  boundary 
is  observed  at  4.  5  x  10“^  torr  cesium  vapor  pressure  coincidentally  with 
a  change  in  electron  concentration  and  potential  along  the  spherical  plasma 
bunch  radius. 


Shul'gin,  B.  V. ,  A.  L.  Shalyapin,  V.  V. 

Ilyukhin,  I,  Rogovich,  N.  A,  Nosyrev, 

F.  F.  Gavrilov,  and  N.  V.  Belov.  Na-Oxy- 
zirconosilicate  (NaaZrSiOO.  Crystalline  structure 
and  optical  spectra.  DAN  SSSR,  v.  202, 
no.  5,  1972,  1068-1070. 

A  spectrographic  study  of  Na-oxyzirconosilicate  is  reported, 
indicating  its  potential  as  a  laser  material.  A  complete  analysis  of  spectro¬ 
scopic  characteristics  of  the  Na2ZrSi05-Eu3+  matrix,  dependable  character¬ 
ization  of  the  luminescence  spectrum  of  Eu3+,  calculation  of  intensity  of 
spectral  lines  and  the  magnitude  of  Stark  splitting  of  energy  levels  of  Eu3+ 
ion,  are  accordingly  given.  The  authors  discuss  the  crystalline  structure 
of  the  matrix  and  possible  isomorphic  replacement  of  Zr  by  Eu3+  ions  which 
point  to  the  existence  of  two  types  of  Eu  centers;  this  was  experimentally 
confirmed  by  observing  lines  in  Na2ZrSi05-Eu3+  luminescence  spectrum 
characteristic  of  the  two  types  (Euj  and  Euh)  of  centers.  The  method  for 
identifying  spectral  lines  corresponding  to  Euj  and  Eun  is  discussed. 
Inasmuch  as  the  intensities  of  spectral  lines  and  the  magnitudes  of  Stark 
splitting  are  approximately  the  same  for  both  cenfers,  experimental  and 
theoretical  results  are  given  only  for  one  of  these.  Figure  1  shows  the  lumi¬ 
nescence  spectrum  of  Eu  +  in  Na2ZrSi05»  Excitation  spectra  for  80°K  and 


V-*/. 


Fig.  1. 
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295°K  were  also  studied.  Results  are  presented  in  Figure  2,  where 
the  excitation  maxima  for  T  =  80°K  are  marked  by  the  heavy  black  lines. 


Fig.  2." 


rj\*i2aa 


0 


Energy  levels  for  Eu3+  ion  in  Na2ZrSiOg  were  computed  on  the  basis  of 
luminescence  and  excitation  spectra  and  are  shown  in  the  right  hand  side 
of  Figure  2. 


Experimental  and  calculated  magnitudes  of  Stark  splitting  of 
'F1  and  7F2  energy  levels  of  Eu3+  ion  in  Na2ZrSi05  are  given,  as  well  as 
the  relative  intensities  of  emission  transitions  of  Eu3+  in  Na2ZrSiOs  crystal 
from  the  5Dq  energy  level.  Authors  also  point  out  that  the  population  of  7F2 
energy  level  by  Eu3+  ion  in  Na2ZrSiC>5  at  room  temperature  is  considerably 
smaller  than  in  ZrSiO^  or  YNbO,^  matrices  which  is  an  important  factor  in 
considering  Na2ZrSiC>5  for  a  basic  laser  material.  Judd's  (B.  B.  Judd,  Phys. 
Rev.,  127,  750  (1962))  and  Ofelt's  (G.  S.  Ofelt,  J.  Chem.  Phys.,  37,  511 
(1962))  techniques  were  used  to  calculate  the  crystalline  field  parameter* 
in  Na^ZrSiOg  and  the  relative  intensity  of  lines  in  the  emission  spectrum 
of  Eu3+. 
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SUBJECT:  RECENT  SOVIET  ACTIVITIES  IN  UNDERSEA  RESEARCH 
VEHICLE  DEVELOPMENT 


SOURCES:  (1)  Bezdomnyy,  N.  The  "Odissey"  puts  to  sea.  Vokrug 
sveta,  no.  3,  1972,  50-53. 

(2)  Kovadlo,  M.  Afalina  studies  the  sea  bottom.  Bakinskiy 
rabochiy,  22  March  1972,  p.  4,  col.  1. 

(3)  New  equipment  for  science  and  technology.  Sovetskaya 
Moldaviya,  13  August  1971,  p.  4,  cols.  7-8. 


One  of  the  first  concrete  articles  (1)  on  the  full -configuration 
use  of  the  most  sophisticated  Soviet  undersea  research  vehicle  (URV), 
the  Sever-2,  appeared  recently  and  is  accompanied  by  several  photo¬ 
graphs.  The  article  gives  a  popularized  account  of  the  Sever-2's 
acceptance  test  (July  1971,  in  the  Black  Sea)  with  almost  half  its  emphasis 
on  a  general  description  of  the  R/V  "Odissey",  the  mother  ship  and  vessel 
used  to  transport  the  Sever-2.  The  Odissey  is  a  converted  BMRT 
(large  freezer /trawler)  which  has  been  fitted  with  a  large  "hangar"  on 
the  port  side  near  the  stern.  This  hangar  houses  the  Sever-2  in  a  con¬ 
trolled  environment,  and  large  doors  give  access  to  the  port  side.  In 
launching  the  Sever-2,  the  travelling  crane  operator,  sitting  in  a  glass 
enclosure  over  the  hangar,  controls  the  entire  launch  procedure,  includ¬ 
ing  control  of  ballast  tanks  selectively  flooded  to  provide  list  to  port 
during  launching.  Photographs  in  the  article  show  a  top  view  and  bow 
view  of  Sever-2  being  launched  (or  recovered),  with  two  extendable 
booms.  The  photographs  indicate  that  freeboard  to  the  hangar  deck 
is  quite  small  (perhaps  8--10  ft.).  This  would  indicate  that  the  doors 
are  watertight  to  protect  the  controlled-environment  hangar  against  sea 
water  entry.  Sever-21  s  actual  design  is  quite  faithful  to  previously 
published  photos  of  a  model,  with  the  exception  of  four  housings  on  the 
deck.  Judging  from  the  photographs,  these  comprise  the  port  and  starboard 
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vertical-propulsion  systems.  The  after  housing  consists  of  what  appears 
to  be  either  a  watertight  electric  motor  or  gear  box,  a  four-blade 
propeller,  and  a  Kort  nozzle  mounted  on  a  shaft.  This  shaft  is  connected 
to  the  forward  housing  and  allows  the  vertical  propulsion  systems  (port 
and  starboard)  to  be  swung  outward  90  into  their  operating  positions 
or  "tucked  in"  over  the  hull  during  launch  or  recovery.  It  is  not  yet 
known  whether  this  motion  can  be  controlled  from  within  Sever-2  or 
whether  it  is  performed  manually  on  deck  by  the  launch/recovery  team. 

As  Sever-2  submerged,  it  was  tracked  by  the  Odissey's 
echo  sounder  and  graphic  recorder.  Aboard  Sever-2,  an  open-mike 
tape  system  was  used  to  record  all  commands  and  other  sounds 
throughout  the  vehicle  and  hydroacoustic  communication  was  used 
between  Odissey  and  Sever-2. 

The  acceptance  program  was  kept  simple  and  amounted  to 
submerging  vertically  to  a  specific  depth  (probably  2000  m  -  Sever-21  s 
rated  operating  depth),  measuring  water  parameters  from  the  surface 
to  the  above  depth,  and  checking  the  operation  of  the  instruments  and 
the  integrity  of  the  hull  penetrations.  It  is  interesting  to  note  that  at  a 
depth  of  somewhere  below  1000  m,  when  Sever-2  reached  equilibrium 
in  the  denser  water,  it  was  decided  not  to  risk  additional  ballast-tank 
flooding  by  pump.  At  this  point,  the  vertical  propulsion  system  was 
used  to  drive  Sever-2  downward  to  the  final  test  depth  (it  is  mentioned 
that  the  last  50m  took  a  half  an  hour).  Recorders  plotted  current 
speed,  water  temperature,  salinity,  density,  transparency,  and  other 
parameters.  After  reaching  the  desired  depth,  the  Odissey  was 
informed  that  Sever-2  showed  no  leaks  or  damage  and  was  ready  to 
ascend.  During  ascent,  the  vertical-propulsion  motors  were  used 
again.  At  1000m,  ballast-discharge  pumping  began  and  continued  to 
about  300m  at  which  point  Sever -2  began  horizontal  movement.  After 
a  brief  excursion,  the  vertical -propulsion  motors  were  used  to  bring 
Sever-2  to  the  surface. 

Two  principal  figures  associated  over  the  years  with  the  devel¬ 
opment  and  designing  of  the  Sever-2  and  other  Soviet  UR V's  were  involved 
in  the  acceptance  tests.  The  first,  A.  N.  Dmitriyev,  served  as  the 
Chairman  of  the  State  Commission  for  the  Sever-2  tests,  while  the  second, 
M.  N.  Diomidov,  was  the  vehicle  pilot. 

Although  earlier  indications  were  that  Sever-2  was  to  carry  a 
crew  of  three,  the  vehicle  now  reportedly  accomodates  four  crewmen. 
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A  very  brief  news  item  (2)  appearing  in  a  March  newspaper 
indicates  that  a  small,  relatively  unknown,  two-man  TJRY,  the  Afalina, 
is  presently  being  used  for  bottom  mapping  ("topographic  surveying"  in 
the  original  article)  and  for  studying  bottom  flora  and  fauna.  The 
Afalina  was  designed  and  built  by  students  at  the  Leningrad  Shipbuilding 
Institute  and  was  awarded  a  prize  at  the  permanently  operating  Exhibition 
of  the  Achievements  of  the  National  Economy  in  Moscow.  The  only 
previous  information  on  Afalina  appeared  in  an  August  1971  newspaper 
item  (3),  in  which  it  was  implied  that  the  vehicle  had  only  been  designed 
and  produced  in  a  small-scale  model  (shown  in  a  photo)  which  was  on 
show  at  the  above  exhibit.  The  information  available  on  Afalina  is  too 
scant  to  determine  when  construction  actually  began,  but  if  it  is  indeed 
operational,  this  should  constitute  some  sort  of  a  record  for  the  Soviets 
in  getting  a  vehicle  off  the  drawing  board  and  into  the  water.  The  fact 
that  it  was  a  "student"  project,  possibly  not  subject  to  conventional 
government  red  tape,  may  have  contributed  to  the  apparently  short 
development  time.  In  this  regard,  it  ie  interesting  to  note  that  Sever-2, 
discussed  above,  was  through  its  design  phase  (first  model)  in  1963  and 
had  a  projected  delivery  date  of  1965. 
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